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Electromagnetics Introduction 



K. M. Siegel* 



Professor Keller was originally chosen to lead the 
single-body scattering and diffraction effort, but 
because he spent the latter part of the period in 
Brazil he was not able to make a contribution to this 
report. The final reports on rough surfaces by 
Anient and on coherence theory by Zucker are not 
yet available at this writing, but should be available 
at least as supplementary material at the General 
Assembly in Japan. 

1. The Progress Report on Antennas, 1960 to 1962, 
has proven to be an impressive document. It points 
out new results obtained on log periodic antennas, 
on conical and multi-element spirals, and on the 
utilization of small antennas as elements of large 
arrays. It describes a significant expansion of effort 
in electrical scanning and array techniques, pointing 
out the relative decline in favor of ferrite phase 
shifters and frequency variation techniques for the 
scanning of arrays. Those devices and techniques 
have been replaced by the utilization of multiple 
beams and digital arrays which offer more versatility 
with little or no increase in complexity. Significant 
advances have been made in the utilization of non- 
uniform arrays which have fewer elements than 
uniform arrays without sacrificing the gain of the 
main lobe. It is quite clear that an increasing 
amount of future effort will go into nonuniform 
arrays, especially for the purpose of reducing side 
lobes. It is also clear that low-noise antenna 
systems have made significant advances. 

This time period has been marked by the first 
significant publications in the United States litera- 
ture on synthetic arrays, and the use of these arrays, 
especially in view of their application to ground 
mapping, definitely represents a major advance in 
antenna theory. The large buildup of information 
on this subject during this period, however, is 
deceptive since much of the work was done many 
years ago but is now reaching the open literature for 
the first time. Of course, the principle of synthetic 
arrays has been familiar in radio astronomy for some 
time, but the new application of this principle to 
ground mapping by airborne radar has made it 
possible to obtain absolute azimuth resolution with- 
out sacrificing gain. It is clear that absolute 
azimuth resolution will allow radar techniques to 
compete for the first time with more conventional 
techniques associated with certain portions of the 
optical spectrum. 

2. Statistical problems in electromagnetics have 
received an increasing amount of theoretical at- 
tention, especially in the field of partial coherence 
and the construction of possible distribution func- 
tions to model rough surfaces mathematically. 

The utilization of statistical methods to interpret 
scattering data from the moon has not as yet effected 
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a prediction of the lunar return for any case not 
previously measured. Of course, with enough ar- 
bitrary parameters present in the statistical theory 
one can almost predict the results of all previous ex- 
periments. One of the first tests of the utility of 
statistical theories of radar scattering should come 
from experimental work on the scattering of electro- 
magnetic waves by planets. It can be expected that 
the next three-year period will see a multiplication 
of the effort in this area by a big factor in comparison 
to that of the last three-year period. Critical anal- 
yses have been made which increase our under- 
standing of the effect of turbulence and other ran- 
dom fluctuating media on electromagnetic fields. 

New approximations have been introduced in the 
studies of propagation in random media and diffrac- 
tion from rough surfaces. The use of impedance 
approximations for such problems has been found 
justified and improved results have been obtained. 
However, many of the approximate solutions in- 
volving impedance boundary conditions with differ- 
ent numerical techniques have in fact yielded 
apparently different answers, and it is hoped that dur- 
ing the next three years a better understanding will 
come into existence of which approximate numerical 
methods are the correct ones under various imped- 
ance boundary condition situations. It is, for 
example, clear that when one type of approximation 
is used a roughening of a surface always increases the 
radar cross section of the surface. Yet when a differ- 
ent approximation is used this is not the case for an 
instantaneous solution, and only if one takes an 
average over many different instantaneous solutions 
does it become the case. It is for this reason, among 
others, that confusion exists not only in connection 
with the utility of particular methods but in the 
arguments presented by different authors on results 
obtained presumable for the same models. 

The work of partial coherence has been aimed 
primarily at the region of optical wavelengths al- 
though it is beginning to find application to other 
portions of the spectrum, especially in radio astron- 
omy. It is clear that if appropriate data processing 
techniques are used after measurements are made 
both of the amplitude and the phase of signals, a 
significant increase in information on such quantities 
as the size of objects being measured and the distance 
between random sources will result. The theory of 
partial coherence should find further application in 
optical analog computing devices, soon to be instru- 
mented to supplement, and, in some cases, replace 
digital computer techniques when large quantities of 
information must be processed. 

3. The work on radiation has pretty well followed 
the predicted patterns. The theory of guided waves 
especially has developed as expected. Much atten- 
tion has been given to leaky and surface wave 
phenomena. Periodic structures and structures 
loaded with anisotropic media have been investigated 
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further. Military and space research requirements 
apparently have stimulated many investigations 
connected with the propagation and generation of 
electromagnetic waves in plasmas. Studies and ap- 
plications have been made of special features in- 
volved in the propagation of waves through a 
plasma, such as the possibility of obtaining a di- 
electric constant with a negative real part. Much 
work has been devoted to the propagation of electro- 
magnetic waves in lossy media because of an in- 
terest in communications between underground or 
underwater terminals. It is quite clear that these 
efforts will continue as long as people are discussing 
hardened sites, hardened communications systems, 
and as long as people are worried about nuclear 
attack. 

4. In the field of classical diffraction and scatter- 
ing, considerable efforts have been made involving 
passive and active reflectors and especially their 
applications to satellites. This is recognized in the 
present report by the inclusion of a paper of J. R. 
Burke of NASA who discusses passive communica- 
tion satellites and the Project West Ford launches 
of 1961 and 1902. The general field of passive and 
active reflectors has been covered in a paper originally 
submit ted in two separate versions by J. Kaiser of 
the Institute of Defense Analysis and I. Kay of 
Conductron Corporation. These have been com- 
bined into one paper of eliminate duplication. 

The many-body problem has continued to receive 
attention and this work is summarized within. 

In the field of diffraction and scattering we find 
that the usual problems are receiving ever increasing 
attention, and workers are obtaining new approxi- 



mate solutions by old methods and old solutions by 
new methods. The writer predicted the last time 
that this report was written that there would be a 
considerable effort made for the resonance region 
during the years 1960 to 1963. However, few new 
results yet exist for the resonance region, and it 
becomes apparent that the love for the asymptotic 
theories is an ever consuming one. Better and more 
accurate answers continue to accumulate for the 
region we already know a considerable amount about, 
but people are not generating methods for use in the 
region where little is really known. There is a ray 
of hope, however, in that more people are becoming 
interested in obtaining values for the currents on a 
body, and in fact this interest has led to the investi- 
gation of a new kind of problem during the period. 
This problem is based on introducing an artificial 
boundary condition to simulate such structures as 
fine wire gratings: a requirement that boundary 
currents be unidirectional. This procedure is similar 
to one previously used to analyze traveling wave 
tubes. 

The electromagnetic theory held is growing in the 
United States. The volume of information is in- 
creasing. The harder problems, however, are still 
not really being attacked; only in antenna theory 
can we see that really significant advances have been 
reported as compared with the last period. 
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Part 1. Progress in Antennas 1960-1962 

V. H. Rumsey, 1 R. C. Hansen, 2 and A. Ksienski 3 



Progress in the United States on antennas during 
the triennium 1960-62 is reported under three cate- 
gories: frequency independent and small antennas; 
arrays and electronic scanning, near fields, aperture 
antennas; data processing and synthetic aperture 
antennas. The survey is based on the definitive 
papers with a bibliography of supporting literature. 
All important U.S. journals have been covered. 

A number of new log periodic antennas have been 
invented consisting of TEM biaxial or coaxial lines 
which are loaded with monopole or dipole antennas. 
Conical and multielement spiral antennas have been 
developed giving circular polarization in all directions 
and an interesting variety of amplitude and phase 
patterns. Solutions of Maxwell's equations have 
been worked out for idealized spiral and periodic 
antenna types which show how the effective aperture 
is controlled by the geometry as well as predicting 
some unexpected and unverified effects. 

A more precise and thorough investigation of the 
characteristics of small antennas for use as elements 
of large arrays has been carried out. 

The ferrite phase shifter and frequency scanned 
arrays have declined in popularity, because multiple 
beam and digital arrays offer more versatility with 
little or no increase in complexity. Effort in the 
latter areas will probably expand, especially on dis- 
tributed transmitter arrays, and on receiving arrays 
with an amplifier/processor at each element. Digital 
array control may be the ultimate for multiple 
targets and target discrimination.. Nonuniform 
spacing of elements has come of age, but very few of 
the design techniques produce low sidelobes; a con- 
spicuous exception is the space tapering approach 
where 20 db has been achieved. A theoretical 
analysis and synthesis technique is badly needed to 
allow synthesis of low sidelobe and other practical 
designs. Automatic steering systems for arraying 
large dishes, for automatic angle return, for atmos- 
pheric scintillation compensation, and for beaming 
large amounts of RF power represent one of the most 
exciting developments. Much work remains to be 
done in this area. 

The analytical understanding of radiation and 
near fields, of interest particularly for power transfer 
and high power focusing, has improved markedly. 
Promise for millimeter, submillimeter, and laser 
applications is offered by the beam waveguide. 

Low noise antenna systems have made significant 
advances, especially in hardware techniques. Use 
of the Gassegrain geometry and of feed shaping has 
produced antenna temperatures well below 20 °K. 
These systems will continue to be important for 
space and satellite communications. Several differ- 
ent types of large radio astronomy instruments have 
come to fruition: the fixed cylindrical parabola with 
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scanning linear array feed; the standing parabola 
with tillable flat plane reflector; a Christianson 
array of dishes, in a cross form. The unfortunate 
cancellation of the Sugar Grove 600-ft dish project 
is a loss to the science of space communications. 

The last three years saw an increasing interest in 
unconventional techniques for satisfying the de- 
mands imposed on antenna systems. The synthetic 
array is an excellent example of what can be achieved 
if an antenna system is optimized for a particular 
application, i.e., ground mapping by means of a 
moving vehicle. This type of specialized solution 
will become more common as the conversion of 
electromagnetic (or acoustic) signals into required 
data is considered as a single continuous process 
which is to be optimized jointly for the desired final 
objective. The generalization of the analytical 
treatments to include at once all dimensions of 
interest and the inclusion of processing as part of 
the antenna proper are symptomatic of this develop- 
ment. The concepts of feedback and self-adaptive- 
ness will become familiar in antenna design, and as 
the characteristics of various data processing arrays 
are better assessed they will provide the versatility 
of solutions needed for the ever increasing variety of 
objects. 

U.S. National Committee reports on antennas for 
the previous two triennia have been published and 
may be of interest: Biekmore, R. W., and R. C. 
Hansen (Nov.-Dec. 1960), J. Res. NBS 64D (Radio 
Prop.), No. 6, 731-741; Cottony, H. V., et al. (Jan. 
1959), IRE Trans. Antennas Propagation, AP-7, 
No. 1, 87-98. 

1.1. Frequency Independent Antennas 

Since the ideal of frequency independent antennas 
was first mooted in 1954 at the University of Illinois 
and successfully demonstrated there in 1956, much 
progress has been made in the understanding and in 
the development of these antennas. They can be 
classed as spiral types or log periodic types. Roughly 
speaking the distinction is that the former give 
patterns which, while they remain constant, rotate 
around the spiral axis with frequency, while the latter 
give patterns which change but repeat when the 
logarithm of the frequency is increased by a fixed 
number. In practice the rotation in the former and 
the variation within a period in the latter are negligi- 
ble effects, so that the patterns and impedances are 
frequency independent over a band which is limited 
at the high end by the size of the feeding transmission 
line and at the low end by the size of the antenna. 
Bandwidths of 20:1 were obtained in the early models 
and these do not appear to have been exceeded since 
then presumably because there are no transmitters 
and receivers which could make use of such wide 
bandwidths. In this connection it is worth noting 
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that a beautiful solution to the problem of designing 
a very wide band balanced to unbalanced 1 ransformer 
has recently been published [Duncan and Minerva, 
I960]. 

Many valuable forms of log periodic antennas 
have been invented during- the last three years such 
as the log periodic dipole [1 shell, I960; Mayes and 
Carrel, 1961] and monopole [Berry and Ore, 1961; 
Greiser and Mayes, 19011 arrays, the tapered ladder 
which is fed by capacitive coupling to a transmission 
line [Wickersham, et al., 19611 and the zigzag which 
is fed directly from a parallel wire transmission line 
[Can, 1961]. Unlike the original log periodic an- 
tennas these all take the form of a conventional 
TEM biaxial or coaxial line which is log periodically 
loaded with radiating elements. The successful de- 
velopment of these inventions has been based on 
simple ideas combined with much practical experi- 
ment. Although no high gain forms have been 
found, gains of up to 17 db (± ?) have been reported 
for the dipole array type. 

There have also been some interesting measure- 
ments on the current distribution and field in the 
vicinity of typical log periodic antennas of the earlier 
inclined zigzag type [Bell et ah, I960]. They show 
the existence of a nonradiating wave traveling out- 
ward from the terminals which gets converted in the 
"active region" to a radiating wave traveling inward 
with the orthogonal polarization. A typical charac- 
teristic of the active region is that the currents in 
adjacent elements of the zigzag are 90° out of phase. 

There has also been a sizeable effort on the theory 
of log periodic antennas, although no solution of 
Maxwell's equations for even the most rudimentary 
log periodic structure has been discovered. The dipole 
array type has been analyzed by ignoring all hut a 
finite number of dipoles in the active region, and so a 
useful body of design information has been worked 
out [Carrel, 1961]. Another approach assumes thai 
the solution is approximately obtained by substitut- 
ing an expanding amplitude and period into the 
solution for a simply periodic structure, a reasonable 
approximation for slowly expanding log periodic 
antennas [Mayes et al., 1961 ; Oliner, 1963; Rumsey, 
1963]. This admits the application of the well 
developed periodic structure theory particularly the 
Brillouin diagram methods. A solution of Maxwell's 
equations has been obtained for the modes of a 
plane sheet of sinusoidal wires and, despite its far- 
fetched connection with log periodic antennas, shows 
remarkable agreement with amplitude and phase 
measurements of currents on log periodic zigzags. 
More general arguments based on the properties of 
periodic structures without regard to Maxwell's 
equations have contributed significantly to the 
understanding of log periodic structures [DuHamel, 
1963; Mittra, 1962].' In this way, it is seen that the 
attenuation of current with distance from the center 
of a log periodic antenna sets in when the propagation 
constant of the periodic structure becomes complex. 
At present there is a vigorous controversy on whether 
the attenuation effect is comparable to the slow wave 
stop band effect or the fast wave radiating effect 



associated with open periodic structures. It does 
appear that this approach gives a good idea of the 
current distribution in the active region, in much the 
same way that uniform transmission line theory gives 
a good idea of the current on straight wire antennas. 

It was recognized at the outset when frequency 
independent antennas were first bruited that not all 
log periodic structures could be truncated; i.e., the 
current was not certain to be attenuated. This still 
remains a fundamental question but the periodic 
structure method is getting very close to an answer. 
Another important question still more or less un- 
resolved is how to minimize the variation of pattern 
within a period. 

Turning now to the spiral types of frequency 
independent antenna, here also there have been 
several practical and theoretical discoveries in the 
last three years. Beginning with the theoretical side 
we note that Maxwell's equations have been solved 
for an anisotropic sheet of equiangular spiral wires 
[Cheo et al., 1961]. This is a far cry from the 
practical two element antennas, having instead an 
infinite number of filamentary elements, but it does 
represent many of the essential features of the an- 
tenna and the theoretical results agree amazingly 
well with pattern and polarization measurements. 
The solution shows that for large / the magnetic 
field at the antenna, i.e., surface current density, 
decreases as r~ 2 . The reduction in current with r 
depends markedly on the spiral curvature, but much 
less markedly on the distance along the spiral: typi- 
cally the current is reduced to 0.1 of the input current 
when r=0.3A and the spiral radius expands by a 
factor of two per turn. The polarization is circular 
in all directions; it is a fact equally as amazing as the 
frequency independent pattern, that this appeal's to 
hold also for the measured radiation from many plane 
and conical spiral antennas. 

One of the surprising results of the theoretical 
solution is this. When the excitation is such as to 
produce an outward current wave at the input, the 
current for large /' consists of an inward wave, but 
the pattern and current amplitude are the same as 
for the reverse excitation. Measurements of the 
phase in the distant field of two element antennas 
always correspond to the solution with an inward 
current wave near the input, which may thus be 
called the natural sense of the antenna. Attempts 
to drive four and six element antennas in the un- 
natural sense have produced some remarkable results 
[Dyson et al., 1961; Sussman, 1961]. Measurements 
of the continuous variation of phase around the axis 
showed that near the input only the unnatural mode 
was present but in the radiation field the direction 
of circumferential phase velocity was reversed and 
consisted again of only one mode, namely the lowest 
harmonic which would fit the given (discrete) set of 
phases at the input terminals [Deschamps, 1959]. 
An input variation of 2t thus gets converted to 
2ir(l — n) in the radiation field of an n element 
antenna. Why this is so is a most interesting and 
unresolved question. A variety of patterns has been 
obtained from four and six element antennas, some 
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of which give excellent omnidirectional coverage in 
the horizontal plane with maximum gain and circular 
polarization. 

Summary 

A number of new log periodic antennas have been 
invented consisting of TEM biaxial or coaxial lines 
which are loaded with monopole or dipole antennas. 
( )onical and multielement spiral antennas have been 
developed giving circular polarization in all direc- 
tions and an interesting variety of amplitude and 
phase patterns. Solutions of Maxwell's equations 
have been worked out for idealized spiral and periodic 
antenna types which show how the effective aperture 
is controlled by the geometry as well as predicting 
some unexpected and unverified effects. 

Small Antennas 

A considerable part of the effort on small antennas 
has gone into more precise evaluation of their effect 
on the performance of arrays. The cause of this 
appears to be the fact that sidelobes nominally 
designed to be 40, 50, °o db down rarely turn out to be 
more than 30 db down. Of the many obvious causes 
for this, mutual impedance is probably dominant. 
Theoretically, it can be taken into account, but the 
practical difficulties of doing so are formidable. For 
an infinite uniformly spaced array of identical 
elements, so energized as to launch a plane wave in a 
single direction, the analysis can evidently be reduced 
to consideration of a unit cell of the structure. In 
this way the variation of impedance with direction of 
the beam has been worked out [Edelberg and Oliner, 
1960a and 1960b]. Obviously this does not apply 
near the edge of a finite array. This edge effect has 
been thoroughly explored, and extensive quantitative 
results have been obtained, both theoretically 
[Carter, 1960] and practically [Kurtz, et al., 1961]. 
Application of these results to precise low sidelobe 
design has not yet been reported but should give 
much improved performance. 

Better or alternative methods of exciting slots in 
waveguides have been the subject of several contri- 
butions, such as the use of position and orientation 
of the slot to get independent control of amplitude 
and phase [Maxum, 1960], the use of asymmetrical 
irises to control the excitation of colinear slots 
[Dudley, 1961; Tang, 1960] and the variation of 
slot length [Rutz, I960]. The use of islands as 
natural slot antennas for VLF has also been studied, 
with rather negative results [Morgan, 1960; Staras, 
1963]. Chen and King [1961] have developed a dipole 
electromagnetically coupled to a two wire line; this 
is the two wire line analog of the waveguide slot. 

The question of what really is the current distri- 
bution on antennas of assorted shapes and sizes 
persists with pristine fascination, the cylindrical 
antenna retaining its preeminence. It still stands as 
an elusive mathematical problem to which numerous 
solutions have been worked out, but none of which 
unfortunately is exactly correct: hence the intract- 
able endurance of the problem. Some interesting 



new results for the current density near the input of a 
cylindrical antenna have been discovered [Chen et al., 
1962; Duncan, 1962; Harrison et al., 1961; Hasserjian 
et al., 1962]. 

Finally, let us note a practical application for super- 
gain arrays [Schildknecht, 1962]. The objection to 
these antennas is their extremely low efficiency, but 
for reception at VLF this may be tolerable because of 
the high background noise, so, if the frequency and 
excitation can be maintained with sufficient pre- 
cision, we may have an excuse to try out this extra- 
ordinary aspect of array theory. 

Summary 

A more precise and thorough investigation of the 
characteristics of small antennas for use as elements 
of large arrays has been carried out. 
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1.2. Arrays and Electronic Scanning 

The triennial period lias seen the appearance in 
the literature of papers describing several advanced 
electronic scanning systems which have been under 
investigation lor several years. In particular, the 
multiple beam array has received considerable 
attention. Shelton and Kelleher [1961] showed that 
independent beams can be formed by a lossless 
passive and matched network only for sin x/x 
patterns or combinations thereof. Beam orthogo- 
nality was shown by Allen [1961, and White 1962] 
to be independent of antenna structure, and to im- 
pose a restrictive relation between beam crossover 
level and sidelobe level. Those restrictions can be 
overcome by use of lossy or active combining net- 
works. Other papers include a matrix approach 
using cross coupling factors [Stein, 1962; Kahn and 
Kurss, 1962]. Hardware realizations include a 
trapezoidal assembly of waveguide couplers (due to 
Blass) with one set of guides connected to elements 
and the other set connected to receivers [Blass, I960], 
and a similar HF system using tapped delay lines 
[Brueckmann et al., 1962]. A more sophisticated 
multiple beam matrix system is the Butler array 
[Butler and Lowe, 1961; Delaney, 1962], which is a 
combining network consisting of rows of hybrids 
with appropriate fixed phase shifters between rows. 
The Butler array requires fewer hybrids or couplers 
to produce N beams from N elements than do 
matrix arrays; // levels of couplers are used where 
w=log 2 N. Effects of random errors in this multiple 
beam system have been investigated by Kiss [1962]. 
Another scheme uses a digital approach, limiting 
sampling and combining the element signals [Gore, 
1962]. The limiting and sampling process introduces 
new sidelobes, but their effect can be negated by 
proper element spacing. The high speed of the digi- 
tal operations allows essentially simultaneous mul- 
tiple beam performance. In addition, control over 
sidelobes and shaped beams is hearty unlimited. 
Other scanning techniques utilize IF heterodyning 
[Cottony and Wilson, 1961; Pels and Liang, 1962], 
and bridge circuits employing variable capacitance 
diodes [Vogt, 1962]. Very little has appeared on the 
formerly popular ferrite phase shifter array [Roush 
and Wiltse, 1961]. The frequency scanned array has 
been analyzed by Ishimaru and Tuan [1962] using 
oj-jS-diagrams; this approach gives insight into the 
interrelationships among the main lobes. 

General studies of electronic scanning effects and 
systems have been prepared by Von Aulock [19(H)], 
Ogg [1961], and Shnitkin [1961]. Some applicable 
material appears in a book by Jasik [1962]. A re- 
lated topic is the gain or effective aperture versus 
scan angle. An approximate value has been ob- 
tained [King and Thomas, 1960; Thomas, 1962], 
and has been extended to cases where polarization or 
impedance are not matched [Tai, 1961]. Allen 
[1962] has measured and calculated mutual imped- 
ance effects under the assumption that the element 
to element mutual impedance or element feed current 
do not change with scan but that the effects of scan 
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are artificially attributed to changes in element 
pattern and directivity. Regardless of the con- 
ceptual approach, the results are that array per- 
formance is insensitive to height of elements above 
a ground plane, and sensitive to element spacing. 
The physical parameters may be chosen to minimize 
the effects of scan changes. For scan angles other 
than broadside, isolation of the order of the sidelobe 
level between transmitted and received signals is 
provided by the fact that the transmitted and re- 
ceived phase fronts are at supplementary angles 
[Rubin and Rabinowitz, 1962]. 

Tracking accuracy has been analyzed from several 
viewpoints. Nester [1962; Leichter, I960], and 
others have related excitation errors to tracking 
accuracy. Brennan [1961; Kirkpatrick, 1962], 
treated a linear array with individual noisy ampli- 
fiers for both amplitude and phase comparison 
systems. His results extend those of Manasse 
[I960] who derived the tracking accuracy (in terms 
of beamwidth) as a function of S/iVand system time- 
bandwidth product. 

A new class of automatic steering arrays has ap- 
peared of which the simplest example is the Van 
Atta array [Sharp and Dial), 1960; Bauer, 1961] and 
subsequently others [Bahret, 1961; Wanselow, 1962] 
proposed passively modulating telemetry on such an 
array (in a satellite) thereby affording space com- 
munications with economy of vehicle power. The 
modulation removes a serious objection to the Van 
Atta array: the pattern nulls produced by inter- 
ference between scattered and reradiated energy. 
Active Van Atta arrays for communications satellites 
were studied by Hansen [1961a, 1961b] who showed 
that their advantage is for those situations where 
narrow beam widths can be utilized; i.e., for trans- 
mission back to the interrogation point, or for space- 
craft at distances greater than 20,000 nmi. 

Turning now attention to pattern synthesis, equal 
sidelobe line source synthesis has been extended to 
the circular aperture by Taylor [I960]. Here, as in 
the case of the line source, a function is constructed 
with the proper asymptotic behavior and with equal 
level sidelobes out to a specified point. Tables of 
this "approximate space factor" for circular apertures 
are available [Hansen, I960]. Cheng and Ma [1961] 
have treated the linear array as a sampled data sys- 
tem in which the pattern is the Z- transform of the 
array distribution envelope. Using this approach, 
equal sidelobe distributions were rederived [Cheng 
and Ma, 1961]. Ksienski derived maximally flat 
and ripple top synthesis techniques for sector-type 
beams, using filter synthesis ideas. However, the 
array polynomial has no poles so that the maximally 
flat synthesis is more difficult than for the Butter- 
worth filter where poles are used to produce a steep 
dropoff [Ksienski, I960]. Other synthesis work in- 
cludes a least squares fit with computer solution 
[Hoffman, 1961] and a shaped beam technique with 
coefficients adjusted for the number of terms used 
[Bricout, I960]. Harrington [1960] extended the 
work of Chu on limitations of antennas, specifically 
on bandwidth and Q as a function of excitation and 



size. Using spherical wave functions he showed that 
any antenna with gain much greater than that ob- 
tained with uniform excitation has a narrow band- 
width, and if lossy materials are utilized, also has 
very low efficiency. The curves presented show 
conclusively that for large antennas supergain is 
impractical. Hansen [1960] considered large antenna 
designs where sidelobes are equal level out to a 
specified point. For a sufficiently narrow main 
beam, the energy in the sidelobes may be comparable 
to that in the main beam, producing a limitation on 
gain. A sin x/x envelope will always obviate this 
problem; he shows where the sidelobe taper must 
start for a given aperture to avoid gain loss. Finally, 
the gain values for Tschebycheff arrays of Brown and 
Sharp's tables (to 40 elements) have been extended to 
120 elements [Stegen, I960]. 

Nonuniformly spaced arrays have aroused much 
interest of late, although the initial work was done 
by Unz in 1955. The objectives of nonuniform 
spacing are to reduce grating lobes, and to use 
fewer elements for a given beamwidth. Random 
removal, symmetrical about the center, has been 
studied by Maher and Cheng [1962]; however, most 
efforts have been devoted to finding an analytical 
or psuedorandom approach that will optimize 
parameters. Several techniques have been tried 
including perturbations from uniform spacing [Har- 
rington, 1961] numerical integration (quadratures) 
such as trapezoidal rule to produce nonuniform 
spacing [Bruce and Unz, 1962a and 1962b; Lo, 1962; 
Maffett, 1962], and computer search [Andreasen, 
1962]. A variety of analytic distributions such as 
log spacing have been tried with little success 
[King et al., I960]. The nonuniform spacing is 
equivalent to an amplitude taper [Willey, 1962], and 
an equivalent uniformly spaced array may be 
defined [Sandler, I960]. Ishimaru [1962] finds an 
equivalent line source using a sum of finite Fourier 
transforms, which for the uniformly spaced array 
reduces to a sum of sin x/x patterns. Each sin x/x 
represents a main beam. This latter approach 
provides some insight into the nonuniform phe- 
nomenon. For 5-element arrays, a solution has 
been obtained which for a given sidelobe level 
yields a shorter array than the Dolph-Tschebycheff 
array for the same sidelobe level [Brown, 1962]. 
However, the beamwidth is correspondingly greater 
in the 5-element nonuniform array. Lastly, Allen 
[1962] has observed that the nonuniformity of the 
mutual coupling in a nonuniform array will have a 
profound effect on the array illumination. 

Summary 

The ferrite phase shifter and frequency scanned 
arrays have declined in popularity, because multiple 
beam and digital arrays offer more versatility with 
little or no increase in complexity. Effort in the 
latter areas will probably expand, especially on 
distributed transmitter arrays, and on receiving 
arrays with an amplifier /processor at each element. 
Digital array control may be the ultimate for 



442 



multiple targets and target discrimination. Non- 
uniform spacing of elements has come of age, but 
very few of the design techniques produce low 
sidelobes; a conspicuous exception is the space 
tapering approach where 20 db has been achieved. 
A theoretical analysis and synthesis technique is 
badly needed to allow synthesis of low sidelobe 
and other practical designs. Automatic steering 
systems for arraying large dishes, for automatic 
angle return, for atmospheric scintillation compensa- 
tion, and for beaming large amounts of RF power 
represent one of the most exciting developments. 
Much work remains to be done in this area. 

Near Fields 

A significant advance in near field theory was 
made by Goubau and Schwering [1961; Christian 
and Goubau, 1961a] with the invention of the iris 
beam waveguide. At a point in the radiating near 
field region, each mode of the field can be recon- 
structed to the original aperture value by a phase 
transformation, thus the field reiterates. The phase 
transformer is a convex lens which retards the beam 
center with respect to the outside; i.e., it refocuses 
the beam. Thus the structure is a periodic beam 
waveguide, capable of carrying power with low loss, 
expecially for beams large in wavelength. The 
phase transformers may be approximated by irises: 
experimental measurements [Christian and Goubau, 
1961b] and higher mode analyses [Beyer and Scheibe, 
1962] have disclosed a loss per iris below 0.1 db for 
spacing of 1 /20 of D 2 /\. ( 'losely related to the beam 
waveguide is tin 4 work of Kay [1960] on the optimum 
phase and amplitude distribution for maximum power 
transfer between two large 4 apertures. In addition 
to the phase curvature needed, as in focusing for 
power transfer [Bickmore, I960] some amplitude 
changes are required for short separations. Jull 
[1962] has shown that measured sidelobes are 
generally higher than expected for antenna separa- 
tions less than 2D 2 /\, but that for certain combina- 
tions of distance and diameter measured, sidelobes 
may be lower. Near-field calculations by Hu 
previously available only in report form, have now 
been published [Ming-Kuei Hu, 1960 and 1961]. 
Hu uses the small angle Fresnel approximation for 
circular apertures of (1— p 2n ) distribution and obtains 
the field in terms of a pair of Lommel functions of 
two variables. Plots of field amplitude and phase 
are given for a number of fractions of D 2 /\. These 
may be compared with the general Fresnel results 
of Hansen and Bailin [1959]. Recently published 
Russian tables [Dekanosideze, 1960] of these Lom- 
mel functions will extend the usefulness of Hu's 
work. Other near field studies are on focused 
apertures [Sherman, 1962]. 

Summary 

The analytical understanding of radiation and 
near fields, of interest particularly for power transfer 
and high power focusing, has improved markedly. 
Promise for millimeter, submillimeter, and Laser 
applications is offered by the beam waveguide. 



Apsrture Antennas 

The Cassegrain reflector system as a large antenna 
has found considerable use as a low noise tempera- 
ture device and for applications requiring large or 
complex feeds. Hannan [1961] and others [Wilkin- 
son and Appelbaum, 1961; White and DeSize, 1961; 
Martin and Schwartzman, I960] have analyzed the 
Cassegrain antenna and have derived the minimum 
blockage criterion. Minimum blockage occurs when 
the feed and subreflector shadows are equal. Mono- 
pulse arrays have been designed in a Dolph-Tscheby- 
cheff sense with element values obtained by simple 
transformation from a Dolph array [Price and Ilyne- 
man, I960]. Wheeler [1962] has made a semiquan- 
titative analysis of monopulse and defocused feeds 
by considering an infinite Gaussian aperture distri- 
bution. Even with defocusing (for acquisition or 
for wide targets) the pattern is also Gaussian, and 
the difference pattern is the derivative or Rayleigh 
pattern. The maximum gain realizable in the differ- 
ence pattern has been shown to be —2.15 db below 
sum gain, for a rectangular aperture [Hannan, 1961a; 
Kinsev, 1962], Multihorn feeds with more than four 
horns offer advantages in optimizing difference pat- 
tern performance, and have been investigated by 
Hannan [1961b and 1961c; Hannan and Loth, 1961] 
and Ricardi [1961]. A two-dish phase-amplitude 
monopulse system has also been developed [Hausz 
and Zachary, 1962]. 

The geodesic or "tin hat" Luneberg lens has been 
studied analytically and experimentally by Johnson 
[1962] and Walter [1962]. The coma corrected zoned 
reflector proposed and analyzed by di Francia [1961] 
has been studied further both theoretically [Dasgupta 
and Lo, 1961] and experimentally. Such a reflector 
allowed scanning [Provencher, 1960] by feed move- 
ment over about ±10 beamwidths compared to ±4 
for a solid parabola. A stepped parabola has also 
been used in a satellite radiometer to reduce feed 
heating [Richter, 1962]. The steps were flat and 
cylindrical instead of conical (as in the zoned reflec- 
tor) and were closely spaced at the operating wave- 
length. Beam scanning of a parabola by feed move- 
ment has engendered additional work [Sandler, 1960; 
Lo, I960]. Work has continued [Altshuler, 1962] on 
line source feeds for scanning spherical reflectors, 
such as the 1000 ft bowl under construction at 
Arecibo, Puerto Rico. Several millimeter wave- 
length parabolas have been successfully produced 
by the liquid spin casting process [Dawson, 1962]. 
An X-Y or elevation-on-elevation mount has been 
developed for tracking near zenith [Rolinski et al., 
1962]. An important development in wide-band dish 
feeds for telemetry has been a conical lobing dual 
log spiral, with the two spirals symmetrically dis- 
posed about the lobing axis. The two spirals are 
decoupled through use of opposite sense circular 
polarization, and cover a 4:1 band [Jasik et al., 1961]. 
Low noise antennas have been of great practical in- 
terest for space communications [Pauling-Toth, 1962; 
Livingston, 1961]; the hog horn antenna developed 
by Crawford [1961] and others and used with 
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TELSTAR probably represents the highest achieve- 
ment in sidelobe and backlobe control. Schuster 
[1962] and Potter [1962] have achieved an antenna 
noise temperature of 15 °K at 2.4 Gc/s in an 85 ft 
dish using a shaped beam horn feed. The latter 
consists of a horn in ground plane, with a corrugated 
surface wave annulus on the outer portion of the 
ground plane. 

Radio astronomy instruments continue to occupy 
much engineering effort. Kraus and Ko [1961] have 
described a 360 ft standing parabola with ground 
plane and tilting reflector, which has just been com- 
pleted. Swenson and Lo [1961] describe a 600 ftX 
400 ft fixed parabolic cylinder with nonuniformly 
spaced scanning feed array. Bracewell [1961] has 
described a Christianson array of dishes built as a 
cross interferometer. Finally, Reber [1961] has 
written on the history of the cross antenna. 

Summary 

Low noise antenna systems have made significant 
advances, especially in hardware techniques. Use 
of the Cassegrain geometry and of feed shaping has 
produced antenna temperatures well below 20 °K. 
These systems will continue to be important for 
space and satellite communications. Several differ- 
ent types of large radio astronomy instruments have 
come to fruition: the fixed cylindrical parabola with 
scanning linear array feed; the standing parabola 
with tiltable flat plane reflector; a Christianson array 
of dishes, in a cross form. The unfortunate cancella- 
tion of the Sugar Grove 600 ft dish project is a loss 
to the science of space communications. 
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1.3. Data Processing and Synthetic Aperture 
Antennas 

A steadily increasing number of applications 
require extreme angular accuracies. These may be 
specified in terms of pointing accuracy for single 
target situations, such as in satellite tracking and 
communication. Other applications may require 
maximum resolution without high gain, such as in 
radio astronomy where integration time may offset 
possible loss in gain. In some applications, how- 
ever, such as ground mapping, one desires both high 
resolution and gain. In the conventional antenna, 
such as a lens, reflector, or linear array, gain, resolu- 
tion and pointing accuracy are rigidly connected and 
there is no way of optimizing one aspect of the per- 
formance at the expense of others. Thus, the con- 
ventional approach could solve these problems only 
in terms of very large apertures and/or very high 
frequencies, with the accompanying tolerance prob- 
lems. In the last decade several unconventional 
solutions were obtained by various researchers; these 
new techniques may be termed collectively data 
processing antennas. These techniques indeed ex- 
ceed performance of the conventional antenna, but 
most of them, e.g., nonlinear antennas, do so in re- 
gard to a specific characteristic while sacrificing 
others. Certain techniques introduce a new param- 
eter such as time or bandwidth and thus avoid major 
parameter tradeoffs. 

Contributions in the field which appeared in the 
literature up to 1960 were reported in the last URSI 
report. A relatively detailed description of these 
techniques is also available in the literature [Ksienski, 
1961a and 1961b]. During the last three years, a 
significant portion of the research activity, theo- 
retical and experimental, has been concentrated in 
the evaluation of these techniques under various 
operating conditions. Several new techniques also 
appeared which were particularly applicable to radar 
and communication, and most significantly, the 
synthetic array radar appeared for the first time in 
the open literature 10 years after its inception. The 
continuing effort to extract more information from 
the signal arriving at the antenna led to the analyti- 
cal investigations which attempt to determine the 
ultimate capability of the processing arrays as well 
as that of the conventional antenna. Information 
theory, detection theory, parameter estimation, and 
other techniques which were successfully applied in 
the past to communication and radar systems are 
being applied to the antenna in an attempt to include 
the angular variables with the range variable and 
formulate a more general set of performance criteria 
and to optimize the antenna performance as a part 
of the overall system. 



A rather interesting and useful concept introduced 
into the antenna system is that of adaptiveness. 
Thus, the antenna must be capable of changing some 
of its parameters in response to the signal and/or the 
environment. This idea was incorporated in a self- 
focusing antenna array [Breese et al., 1961 ; Schrader, 
1962], where each individual element of the array 
continually adjusts the phase shifter (or delay line) 
so as to combine its signal contribution coherently 
with those of the other elements. This scheme pro- 
vides optimum signal reception for a given number 
of antenna elements independent of their relative 
locations and spacing (except for mutual effects), 
and can be used to compensate for atmospheric ir- 
regularities. The signal reception is optimum, due 
not only to the fact that all signals add in phase, but 
that the detection process for the array as a whole 
is phase coherent. This antenna system is excellent 
for communication purposes where little interference 
is present; however, it is not capable of discrimi- 
nating against other sources (or returns) which oper- 
ate at the same frequency as the communicating 
source. Further, signal components within the 
phase lock loop pass bands may cause spurious per- 
formance. This technique was experimentally tested 
with a two-element array [Breese and Sferrozza, 
1962] with performance close to the predicted. The 
self-focusing technique is a good example of what can 
be achieved if one aspect of an antenna array is 
optimized at the expense of others. Here, resolution 
was ignored in favor of gain or SNR. 

The general problem of optimum excitation of a 
linear array which has often been investigated in 
the past has received further attention. Myers 
[1960a and 1960b] investigated experimentally the 
image quality produced by antennas of various 
excitations and found that the uniform illumination 
is close to the optimum in terms of the ease with 
which both a mechanical as well as a human observer 
can recognize the image. An information theoretic 
study [Young and Ksienski, 1961] of the same 
problems reached similar conclusions, although 
under the condition of zero image noise, which is 
equivalent to a noiseless receiver, the information 
rate was found to be unaffected by illumination. 
These results are similar to those obtained by 
Linfoot [1955] for optical images. One may ap- 
proach the problem of optimizing an antenna in 
terms of separability of returns from two closely 
spaced targets. Thus, Woodward's ambiguity func- 
tion [Woodward, 1953] may be generalized to include 
angular variables. This was carried out by Urko- 
witz et al. [1962] who also computed the generalized 
resolution constants in four dimensions: range, 
range rate, and azimuth and elevation angles. It 
should be pointed out, however, that the ambiguity 
function per se is only a measure of resolution and 
does not predict at what angle of separation (or 
range separation) two targets can be resolved. A 
concept which is most useful in evaluating the 
imaging properties of antennas is that of spatial 
frequency, which is equivalent to the temporal fre- 
quency [Bracewell and Roberts, 1954] in networks. 



446 



For the important question of what is the limit of 
the information obtainable from a finite aperture, 
the answer would be given (except for noise) in 
terms of the spatial frequency bandwidth of the 
antenna. In accordance with Bracewell [1961], 
this bandwidth corresponds to the size of the antenna 
in wavelengths, with zero response beyond it. It 
was shown, however, that no such bound or cutoff 
[Lo, 1961; Walter, 1961; Ksienski, 1962; Chisholm, 
L962] exists, as evidenced for example by supergain 
[Spitz, 1959]. For large antennas such as are used 
in radio astronomy, Bracewell's conclusion is prac- 
tically correct, but for small antennas it may not be. 
The response of nonlinear arrays can also be more 
conveniently evaluated in terms of the spatial fre- 
quence response [Young and Ksienski, 1961; Brace- 
well, 1961]. Multiplication of the outputs of two 
antennas convolves their transforms and widens the 
overall bandwidth, but this does not necessarily 
imply that more information can be passed through 
this filter than was possible before multiplication. 
This depends on the coherence properties of the 
object. Thus, for radio astronomy, where the signals 
arriving from different portions of the sky are 
almost totally uncorrelated, the spatial pass band 
does increase, as shown by Bracewell [1961]. For 
partially coherent or completely coherent signals 
the results are far more complex and resolution im- 
provement is limited. The improvement in reso- 
lution when obtained will still often involve loss in 
gain and SNR. Because of the great importance of 
the signal characteristics on product arrays, a large 
amount of work was done in analyzing the perform- 
ance of these arrays for various signal conditions. 
Drane and Parrent [1962] computed the response of 
a single product array in terms of the mutual co- 
herence function of the object distribution. The 
coherence function represents the correlation between i 
the signals emanating from two points on the object 
separated by arbitrary distances in time and space. 
The results of the paper show that on the incoherent 
limit, such as in radio astronomy, the single product 
array is appropriate lor mapping extended sources, 
since it is linear in power. Linder [1961] investi- 
gated the response of single product arrays to both 
perfectly coherent (monochromatic) and incoherent 
signals and computed the effect of integration time 
on the response in the case of single product array 
to uncorrelated targets. Mays and Cheng [1960] 
computed the effect of varying the RC constant of 
an averaging filter following a single product array. 
Linder, as well as others [Pedinoff and Ksienski, 
1962], conclude that multiple product arrays are of 
limited value for even completely incoherent ex- 
tended sources or multiple target situation. The 
response of various single product arrays to two 
perfectly coherent in-phase targets was considered 
in some detail by Pedinoff and Ksienski [1962], and 
it was shown that even for the restricted in-phase 
condition, the responses sometimes bear little re- 
semblance to the true target distribution, except for 
one particular array configuration. But even this 
configuration results in distortions when arbitary 



relative phases between targets are assumed [Mac- 
Phic, Pedinoff, and Ksienski, 1962]. These distor- 
tions, however, do occur in linear arrays as well 
when targets of arbitrary relative phase and ampli- 
tude are assumed. A detailed study of the effect of 
varying phase and amplitude in linear and several 
single-product arrays indicates that the average 
boresight error is about the same in most arrays 
considered [Ksienski et al., 1962a and 1962b]. A 
rather interesting experimental study by Welsby 
[1961] shows that at least for sonar application, a 
multiplicative array has resolution capability defi- 
nitely superior to that of a linear array. 

Experimental work on the response of a product 
array to a single target: was performed by Gabriel 
[1961] with results agreeing well with theory. A 
multiple-product array was experimentally investi- 
gated by Davenport and Drane [1961] and by 
Band [1961]. The response that was obtained to 
a single point target agreed with theoretical predic- 
tions. A numerical and experimental study of the 
response of several single-product arrays to two par- 
tially correlated objects was carried out [Ksienski 
et al., 1962a and 1962b] with results indicating that 
improved resolution is achieved by the product 
arrays, as compared to a uniformly illuminated 
square law detected linear array, for practically all 
correlation levels from zero to unity. 

A subject of considerable importance to the prod- 
uct array performance is the correlation of the 
background noise, e.g., clutter. Freeman [1960] 
investigated the systematic error in direction find- 
ing by means of correlation arrays when the back- 
ground noise is partially coherent. This noise co- 
herence also limits the useful integration time for 
(he array; thus, the relation of the maximum useful 
averaging time to the noise coherence is obtained. 
The effect of mechanical versus electronic scanning 
of the antenna elements and/or the array on the out- 
put SNE of a, correlation array was investigated by 
Jacobson and Talham [I960]. In another paper 
[Jacobson and Talham, 1961] these authors consider 
the effect of scanning errors on correlation array 
output SNR. Jacobson [1962] also investigated the 
effect of beam steering by both time and space 
translation on the output of product arrays due to 
background noise correlation. 

In the product arrays discussed above, the illumi- 
nation of the various elements was kept uniform. 
One may try, however, to modify that illumination 
in order to further improve the resolution of the proc- 
essing array. Jacobson [1961] suggested the use of 
two adjacent linear arrays whose elements are so 
weighted that when the outputs are multiplied and 
averaged, one obtains a Tchebycheff pattern. Price 
[1960] proposed a scheme which is essentially equiv- 
alent to a single product and produces a Tchebycheff 
pattern corresponding to a linear array twice the 
physical size of the processed array. Mattingly 
[1960] suggested a nonreciprocal scheme applicable 
to two-way patterns only. The transmit and receive 
patterns differ from each other such as to form a 
product which yields a Tchebycheff pattern of the 
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order of 2n instead of obtaining a two-way pattern 
corresponding to an nth order squared. This 
process is completely linear, and its main disadvan- 
tage compared to reciprocal patterns is some loss in 
gain. MacPhie [1962] combined the two-way mode 
of operation with a multiplicative process on receive 
to obtain a Tchebycheff of order An, compared to a 
square law detected reciprocal array whose response 
corresponds to TV- 

In the above discussion, an implied constraint on 
the antenna was its physical size; thus, one by as- 
sumption would sample only the energy falling in 
the area occupied by the antenna. If that constraint 
is removed, a completely new field is provided — that 
of synthetic arrays. 

A synthetic array, in the present context, simu- 
lates the performance of a conventional antenna 
array by sequentially moving a single element of the 
array to the various positions which would be simul- 
taneously occupied by the elements of the con- 
ventional array. 

There are 2 main areas to which the synthetic array 
principle has been applied. One field is radio 
astronomy, where Ryle and Hewish [1960] success- 
fully synthesized both one- and two-dimensional 
apertures. In this case, the sources are essentially 
incoherent and the synthetic aperture is passive, 
i.e., operates in the receive mode alone. The second 
is in reconnaissance or ground mapping, in which 
the successive returns of an active radar on a moving 
vehicle were coherently combined. The resolutions 
obtained substantially exceeded that obtained from 
the same radar operating in the conventional mode. 
There is little doubt that this technique, wherever 
applicable, such as in ground mapping application, 
yields results far more spectacular than those 
obtainable from other data processing antennas and 
does not involve serious tradeoffs in important per- 
formance characteristics such as SNR. The great 
merits of the technique are probably the main 
reason for the 10-year delay in their appearance in 
the published literature. In fact, some of the 
material on signal processing is still classified. 

The possibility of obtaining synthetic apertures 
was recognized at least as early as 1951 [Sherwin et 
al., 1962], and in 1952 the principle was experi- 
mentally demonstrated, with actual high resolution 
maps being obtained in 1953. The resolution was 
improved by a factor of 10 beyond the one obtain- 
able in the conventional mode of operation. Most 
of the work since that time involved the develop- 
ment of data processing equipment for the imple- 
mentation of synthetic arrays. One of the most 
important requirements is the preservation of 
phase over the whole integration period or over the 
total synthetic array length. Various sources con- 
tribute to deviations from the exact phase required 
for optimum processing. These include equipment- 
instability, erratic vehicle movements, propagation 
distortions due to the medium, and Fresnel zone 
effects because of the extreme length of the syn- 
thetic aperture. Some of these effects may be com- 
pensated by appropriate processing, while others, 



particularly the random phase errors (except for 
those caused by vehicle motion which can be com- 
pensated), will degrade the synthetic array per- 
formance in various ways [Greene and Moller, 1962]. 
The most significant phase deviation is caused by 
the Fresnel zone effect and requires a focusing 
process to prevent destructive interference for 
arrays beyond a certain length. In fact, the vari- 
ous synthetic arrays are divided into two broad 
classes, the focused and unfocused arrays. The 
focused arrays require much more complicated 
processing systems than the unfocused ones, par- 
ticularly since the focusing function is range de- 
pendent. The resolution obtained, however, when 
focusing is performed is theoretically unlimited (up 
to the uncertainty limit) and even practically 
far in excess of that of the unfocused arrays. An 
excellent paper presenting the patterns produced 
by both focused and unfocused synthetic arrays 
is by R. C. Heimiller [1962]. General expressions 
are derived, and numerous performance curves are 
given for various ranges, frequencies of operation, 
and array lengths. A somewhat simplified com- 
parison between the resolution of conventional and 
synthetic arrays is presented by Cutrona and Hall 
[1962]. The discussion and results given in this 
article as well as the one by Cutrona et al. [1961], 
are mostly qualitative, but helpful in understanding 
the principles involved in a synthetic array. The 
synthetic array performance may be described, as 
was done above, in terms of an interference pattern 
produced by the coherent addition of various re- 
turns. An alternate explanation of the performance 
of a synthetic array is in terms of the Doppler 
characteristics of the return signal. It can be 
seen quite easily that the return from a target which 
at a given instant is located exactly on a normal to 
the flight path of the radar will have at that instant 
no velocity component towards the radar and thus 
will experience no Doppler shift. A target which is 
in the same range at the same interval but at a 
slightly different angle to the flight path will return 
a Doppler shifted signal. The ability to resolve the 
two returns depends on the smallest separation 
bandwidth detectable by the system, which in turn 
is inversely proportioual to the integration time 
allowed. Thus, the Doppler resolution is directly 
proportional to the total synthetic array length 
which coincides with the considerations based on 
interference patterns. A paper discussing various 
representations of the synthetic array process is by 
H. L. McCord [1962]. The early development of 
the concepts and the experimental demonstrations 
including the resulting high resolution maps are pre- 
sented by Sherwin et al., [1962]. Also, the signal 
processing used and some of the actual equipment 
are described both by Sherwin et al., [1962] and 
Cutrona et al., [1961]. The article by Sherwin is 
very helpful in understanding the implementation 
of some of the original processing technique. As 
mentioned earlier, uncompensated phase errors may 
be expected to present the major limiting factor in 
the resolution capability of the synthetic array. 



448 



This problem is investigated by Green and M oiler 
[1962], and various effects of random phase errors 
of prescribed statistical characteristics are presented 
graphically. The effects computed include degra- 
dation of gain, beam width and sidelobe structure, 
and beam-pointing scanning effects. The effect 
of random tropospheric phase errors on synthetic 
array is considered by Rondinelli and Zeoli [1962]. 
Both theoretical and experimental results are 
presented. 

An interesting application of the synthetic array 
principle is in the achievement of simultaneous 
scanning. Jf the elements of a fixed linear array are 
successively excited, simulating a traveling source, 
the Doppler effects discussed above are produced, 
i.e., one can distinguish the returns from different 
angular locations on the basis of their varying 
frequency characteristics. This technique was suc- 
cessfully implemented both in the electromagnetic 
field (1961) and in acoustics (1962). An important 
shortcoming of the system is its low gain, but it also 
provides advantages over other scanning methods 
[Hewish and Shanks, 1962; Kuinmer et al., 1962]. 

Summary 

The last three years saw an increasing interest in 
unconventional techniques for sat isfying t he demands 

imposed on antenna systems. The synthetic array 
is an excellent example of what can be achieved if an 
antenna system is optimized for a particular applica- 
tion, i.e., ground mapping by means of a moving 
vehicle. This type of specialized solution will 
become more common as the conversion of electro- 
magnetic (or acoustic) signals into required data is 
considered as a single continuous process which is to 
be optimized jointly for the desired final objective. 
The generalization of the analytical treatments to 
include at once all dimensions of interest and the 
inclusion of processing as part of the 1 antenna proper 
are symptomatic of this development. The con- 
cepts of feedback and self-adaptiveness will become 
familiar in antenna design, and as the characteristics 
of various data processing arrays are better assessed 
they will provide the versatility of solutions needed 
for the ever increasing variety of objectives. 
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Part 2. Statistical Problems in Electromagnetics 

2.1. Rough Surface 
W. S. Ament ! 



The archtypal, rough surface is that of a tideless 
sea. The sea receives plane-wave illumination as 
from a distant elevated radar at grazing angle 0, 
at a surface point observed from above by a distant 
receiving antenna at grazing angle 6 ' , the relative 
azimuth and the polarizations being left implicit. 
A description of the surface is given in terms of 
appropriate statistics, the simplest being the rms 
surface height H above the mean (sea level) plane, a 
length L characteristic of the dominant ocean wave- 
lengths, time-lagged correlations of surface heights, 
slope skewness due to surface winds, and other 
second-order statistics. Basic problems are to 
calculate, as averages over the surface statistics, (a) 
the interferometrically measurable specular reflec- 
tion coefficient R=R (0) with phase referred to the 
mean plane, (b) the closely related surface boundary 
conditions effective for normal modes propagating 
over the surface, (c) the (dimensionless) scattering 
parameter a (6, 6 '), or bistatic radar area per unit 
area of mean plane, which goes into the radar 
equation for predicting mean received power, (d) 
radar and geometric parameters and data-processing 
methods permitting best inferences of oceanograph- 
ically interesting data. This is also the problem in 
radar terrain surveillance and planetary radar design. 

Ideally, one would solve Maxwell's equations for 
each conformation of the surface, from the appro- 
priate field or power flux at the receiver, and average 
over the statistics controlling the conformations. 
The intractable initial problem leads one to deal 
with simplified model surfaces and to use various 
methods of premature averaging. The chief models 
are: 

(1) The actual surface is taken as perfectly con- 
ducting and the statistics are simplified. 

(2) The mean plane is assumed perfectly con- 
ducting, with particular objects or " bosses' ' ran- 
domly distributed over the plane to resemble waves. 

(3) The mean plane is the mathematical locus of 
randomly distributed objects, this two-dimensional 
random array being isolated in space. A major 
subclass of this is (3a), in which the objects are semi- 
infinite parallel structures, such as wires or parallel 
plates with ends or edges in the mean plane. 

(4) The mean plane bounds a half -space filled with 
a statistically described multiply scattering medium. 
Major advances with this model are described else- 
where in this report. 

In 2, 3, 4 the scattering parameters of the indi- 
vidual objects are assumed known and the problem 
is one of properly accounting for interactions before 
or during the eventual averaging. Some phen- 
omenological models highly useful for organizing and 
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explaining radar data in terms of oceanographic 
surface parameters are discussed in the Commission 
II report. 

Dimensionless numbers important in the problem 
are kH, kL, H/L, Ra=kHsmd, Z=sin 2 d/kL, Z' = 
sinds'md' /kL, where &=27r/(radar wavelength in air). 
When both kH and kL are small, quasi-static methods 
apply. This is the case for estimating losses in 
rough waveguides. When Ra is definable and small, 
the surfaces satisfies the "Rayleigb criterion" for 
smoothness, and specular reflection often dominates 
the scattered energy. Radar resolution is high 
and the scattering problem difficult for values of 
Z, 7J near unity; here diffraction over one wave crest 
becomes important in the field illuminating the 
next crest. The useful range of phenomenological 
models is sometimes extended to regions of small Ra 
by considering a wave reflected from the mean plane 
or from the locally smoothed surface as part of a 
coherent illuminating field, but this practice breaks 
down for Z~l. This review deals principally with 
the theoretically interesting and intractable prob- 
lems in which the excitations and the scattering are 
controlled by interactions among surface elements. 

Methods 

One can (in principle) write exact integral equa- 
tions with Green's function kernels for particular 
conformations of model surfaces, and variational 
expressions [Meecham, 1956] of the Schwinger type 
for scattering into particular an ten in as, but the taking 
of averages poses unsolved difficulties. A varia- 
tional scheme retaining its efficacy for optimum field 
or power estimates throughout the necessary averag- 
ing would form a major advance. One is driven to a 
self-consistent calculation: a parameterized position- 
dependent excitation is assumed for all parts of the 
surface. This establishes the scattered field, in the 
neighborhood of a surface point p, as a superposition 
of scatterings by individual objects or parts of the 
surface. This superposition is then averaged and 
added to the originally incident field to give the 
average excitation about p. The parameters of the 
beginning excitation's representation are then ad- 
justed so that the represented and calculated exci- 
tations are identical. (Recent improvements in this 
self-consistent method are outlined in the section on 
multiple scattering. Wiener-Hopf methods are use- 
ful here with models 3a and 4.) 

Having estimated a plausible excitation of, and 
therefore the scattering from, an object at p, the 
problem becomes one of averaging the scattering to 
get a and 11. With no nonreciprocal elements in 
the surface, reciprocity theorems apply: imagine a 
dyadic Green's function G appropriate to a particular 
fixi (I surface with the object or portion at p removed. 
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With this G one can show that the partial contribu- 
tion ApB of an object at p to a field at B due to a 
radar at A is the same as BpA, the partial field meas- 
ured on interchanging transmitter and receiver. It 
follows that this "local reciprocity" is true for 
average fields and powers, and that the scattering 
from a surface point partakes of interactions of the 
same type and importance as the excitation. 

General criteria jor evaluating proposed solutions to 
rough surface problems in "statistical reciprocity" : 
Reversing incident and scattering directions should 
produce no change in a or in R, for surfaces built of 
reciprocal elements. A more general criterion is 
energy conservation, which demands that the inci- 
dent minus the specularly reflected power, per unit 
area of the mean plane, be accounted for in the total 
randomly scattered power implied in a, plus that 
calculated as absorbed in the surface under the unit 
area. A third criterion, applying immediately to 
model 1 and to 3a and 4 when the objects have finite 
losses, is that coherent and random power fluxes die 
out with distance beneath the mean plane. This is 
so obvious with model 1 that it should be examined 
for suggestions as to calculation methods. 

The second law of thermodynamics is violated 
unless a (6, 6') falls off as 6' or faster for 6'— >0; solu- 
tions satisfying reciprocity and within some measure 
of meeting the energy criterion would appear neces- 
sarily to have the correct behavior here [Anient, 
1956]. Next, fields and currents estimated or as- 
sumed on parts of the surface should satisfy locally 
applicable criteria of diffraction theory such as 
Meixner's [1949] edge condition and Raisbeck's 
passivity criterion [Richards, 1959, p. 97]. (For 
model 3a and Z^l, the edge condition should have 
a central importance.) Finally, extrapolations to 
nature must agree with observations; some of the 
empirical data are found in the Commission II report. 

The foregoing discussion is concerned with aver- 
ages taken over regions large compared with L 2 ; the 
more detailed averages necessary for the radar 
resolution, and for time-dependence backscatter, 
seem to have been treated only phenomenologicaly. 

We turn to some conjectures and unsolved prob- 
lems (beyond the methodological ones implied above) 
that should be within grasp. Calculations for ver- 
sions of model 1, 3a, and 4 yield j?(0)— > —1 and the 
factorization a ^j(6)f(d') in the small 0,0' limit. 
Plausibility arguments suggest that this "grazing 
factorization" is universally true when R(9)^> — 1, 
but can one prove it? Chandrasekhar's exact albedo 
calculations [1960] for dilute versions of model 4 
"factor" differently, but, for these, /i = 0. Next, 
assume an ocean with dielectric constant fractionally 
less than that of the air above, permeability being 
unity everywhere. For dead calms, reflection at 
the mean plane is then total for 6 less than some small 
critical angle O . Assuming rounded waves of large 
kL, kH, and kx (rms radius of curvature) , what a do 
we get for O <0<20 O , knowing that double reflections 
can take place, both with total reflection by ray 
theory? This problem is basic for wave-layer 
theories of tropospheric scatter and there are iono- 
spheric equivalents. 



We turn to specific papers. Ament [1960] intro- 
duces the foregoing reciprocity concepts and raises 
the problem of "grazing factorization." For a tilted 
parallel-plate version of model 3a, he uses a self- 
consistent method to find approximate currents in 
the plates, and local reciprocity to find a reciprocity- 
satisfying a obeying "grazing factorization." But 
the edge conditions and energy conservation are not 
satisfied, and R changes on reversal of incident and 
specular direction. For a model 4 surface with 
isotropic scatterers and scalar waves, the method 
yields better results which compare informatively 
with Chandrasekhar's exact albedo calculations 
[I960]. His final paragraph, arguing qualitatively 
that grazing radar backscatter should always be 
independent of polarization, uses an assumption 
suggesting that radar resolution is also independent 
of polarization; this is contrary to observations with 
oceanographic radars. 

Marsh [1961a, b] writes the scattered field as a 
superposition of upgoing plane waves of various real 
horizontal components of propagation vector, and 
asks that the superposition cancel the downgoing 
plane incident wave at each point of the model 1 
surface, which is specified by a spectrum and rms 
surface height H. He expands in powers of H, finds 
a process, based on Fourier transform theory and 
expressed in the form of iterated operators, for 
obtaining coefficients of H n in the scattered-wave 
spectrum. Putting the results back together in a 
formal expression involving operators, he has a result 
amenable to machine methods for finding specular 
and random field components to arbitrary preassigned 
accuracy. The results satisfy reciprocity, energy 
conservation, and grazing factorization [private 
communication]. While the validity of neglecting 
downgoing components of the scattered field between 
waves is suspect [Lippmann, 1953], the neglect may 
be rectifiable; meanwhile the paper appears the first 
to develop the full power of transform methods in 
rough surface problems. 

The only known rough surface model for which all 
scattered waves are everywhere upgoing is one 
having an impedance boundary condition which 
varies randomly over the mean plane. Hessel 
[1960a, b] and Oliner have treated a sinusoidally 
varying reactive impedance boundary condition and 
obtain exact results in full detail. Their study is 
directed at traveling-wave antenna applications so 
they do not consider the random case. Here one 
would have the difficulty of assigning statistics so 
that the surface is everywhere electromagnetically 
passive, i.e., satisfies Raisbeck's criterion everywhere. 

Senior [1960a] discusses an impedance boundary 
condition effective at imperfect and certain rough 
surfaces in terms of the impedance ratio 77 of local 
impedance of a medium to the free-space impedance. 
He shows that for plane and slowly curving bound- 
aries between free space and a medium of high 
refractive index, having small relative variations of 
dielectric constant only, the average impedance 
boundary is best described in terms of an average 77. 
Senior's second paper [1960b] treats a model 1 surface 
of small-scale isotropic roughness. His main result 
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is a local 2X2 impedance matrix connecting the 
components of E and H tangential to the mean 
plane. The matrix, expressed in terms of an average 
77, diagonalizes when preferred axes are chosen; the 
results seem not unlike those of Biot [1957, 1958a, 
1958b, 19601 f ()r a small-scale model 2 surface. But 
the merit of Senior's form is that it applies to curved 
surfaces. To estimate the effect of small-scale 
roughness on backscatter from large spheres, Hiatt, 
Senior, and Weston [1960] average Senior's matrix 
boundary condition over all incidence angles and 
obtain a scalar impedance boundary condition 
expressed in terms of an effective 77: 77^- ikH 2 /L. From 
this condition they estimate the coherent back- 
scatter field through terms linear in 77. They com- 
pare the obsolute-squared result with measured 
backscattered powers, but give no theoretical esti- 
mate of the nonspecular scattered components in 
their measured backscatter or elsewhere. 

Although introducing no new diffraction theoretic 
ideas, Chen and Peake [1961] develop the basic con- 
cepts in the closely related problem of the angular 
spectrum of thermal radiation from a rough ground 
surface of temperature T which is illuminated by 
prescribed thermal sources distributed over the sky. 
Insofar as the surface elements are lossy, absorbed 
radiation must reappear in the form of heat radia- 
tion, and the problem is correspondingly more com- 
plex than the diffraction-theoretic one of accounting 
for the remnants of the incident radiation only. 

Twersky [1962a, b, c, d] gives for models 2, 3, and 
4 several formalisms for self-consistent calculations 
of the specular reflection coefficient R and the mean 
coherent excitation of general objects of the surface. 
Treating the resulting scattered waves as traveling 
in object-free regions, he finds energy conserved, but 
his o-'s appear valid only at, high angles. His R's 
appear as the natural consequence of his self-con- 
sistent approach and are conjectured to be auto- 
matically correct when the scattering by the individ- 
ual objects is correctly calculated. 

Twersky's papers and earlier summary [I960] list 
most of the known earlier papers in the rough surface 
field. Other bibliographies have been prepared by 
Bachynski [1959], Hagn et al. [1961], Lysanov [1958], 
Marsh [19621, and Wolff [I960]. 
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2.2. Electromagnetic Wave Propagation in a Random Medium* 



W. C. Hoffman* 



1. Introduction 

The phenomenon of wave propagation through a 
random medium has, historically speaking, two main 
origins. One stems from statistical mechanics, and 
accordingly, as classical or quantum models are 
assumed, is capable in principle at least, of explaining 
all electromagnetic phenomena in a material medium. 
In practice this approach tends to be limited to one 
single aspect of "randomness," viz, configurations of 
randomly placed scatterers, and to the determination 
of one single parameter, usually the average value of 
some field quantity associated with the incident plus 
scattered fields. Only this approach seems to be 
adequate for analyzing resonance phenomena, such 
as fluorescence, Doppler broadening, resonance ab- 
sorption, anomalous dispersion, and the like. 

The other approach, essentially macroscopic in 
character, emerged from the problem of scattering 
of acoustic and electromagnetic radiation by turbu- 
lent media. The constitutive parameters of the 
medium are not derived from first principles but are 
taken as given, i.e., to constitute a stochastic process 1 
of some sort. In this context, therefore, the electro- 
magnetic field is a functional of the stochastic process 
representing the random refractive index of the 
medium. The theoretical problem thus has both a 
probabilistic as well as an electromagnetic character. 
In much of the past work the probabilistic side has 
often languished in comparison to the electromag- 
netic. A smooth though random variation of the 
medium is an essential characteristic of this model. 

The coverage of this portion of the report is 
restricted to the macroscopic model, the so-called 
"perturbed continuum" problem [Bremmer, 19631 2 , 
while the electromagnetic phenomena associated with 
random configurations of scatterers are described in 
the section on multiple scattering by J. E. Burke and 
V. Twersky. By the same token, only problems 
involving the field within an unbounded medium are 
taken up here, the companion subject of scattering 
from random surfaces being discussed by W. S. 
Anient in another section of this report. Essentially 
then this part of the report deals with the theory of 
wave propagation in still, cold, random, continuous 
media at radio frequencies not greater than 3 Gc/s, 
and the contributions discussed below are by and 
large restricted to those made during the last 
triennium. The greater part of the work prior to 
the current triennium has been referenced in the 
summary "On Multiple Scattering of Waves," pre- 
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1 A stochastic process is a one parameter family of random variables, the 
parameter being time (as in random noise phenomena) or points in space or 
space-time, etc. 

2 The debt that the present writer owes to H. Bremmer for his definitive report 
on "Scattering by a Perturbed Continuum" will be clear to anyone who has 
read the latter. 



pared by V. Twersky [1960], for the USA National 
Committee report to the XT II General Assembly, to 
which the interested reader is referred. 

Section two of this report deals with the mathe- 
matical theory of electromagnetic waves in a random 
continuum. This is followed by sections on ray- 
tracing in a random medium, transport equations of 
various types, and what has come to be known as tur- 
bulent scatter. Reflection and transmission prob- 
lems associated with random stacks of slabs have re- 
ceived considerable attention of late and an assess- 
ment is made of these developments under the same 
head as transport processes. Finally, we take up 
certain mathematical and statistical results which 
may have significance for the analysis of /'-/radiation 
in random media. 

2. Theory of Wave Propagation in a Random 
Continuum 

Nearly all treatments of wave propagation in 

a random continuum have assumed harmonic time 
dependence, so that the reduced wave equation 
V 2 [7+F(x)t7(x)=0 can be employed rather than 
the full wave equation in x and t. This condition 
will be taken to apply in the sequel except where 
specific mention is made to the contrary. As pointed 
out by Silver [1963], this assumption requires that the 
phenomenon be essentially monochromatic, and such 
things as Doppler diffusion of frequency in hot, ener- 
getic media are thereby excluded. In general, use 
of the reduced wave equation appears valid whenever 
the time variation in the properties of the medium is 
much slower than the propagation time of the wave. 
A medium which is swirling just as rapidly as the 
wave is traveling will clearly have frequency shifts 
which invalidate the nionochromaticity assumption. 
The probabilistic nature of a slowly swirling medium 
can be treated either through the time dependence of 
the stochastic process or by limiting consideration to 
realizations which are sufficiently far apart in time 
to be independent. However, as pointed out by 
Silver [1963], an ergodic theorem for relating time 
averages to expectations (ensemble averages) is re- 
quired in either case. 

Perhaps the most important single work on wave 
propagation in a random medium that has appeared 
during the last triennium is that due to J. B. Keller 
[1962]. The first portion of Keller's paper is note- 
worthy for a very lucid discussion of the problem and 
its context. He subsequently covers ray-tracing in a 
randomly inhomogeneous medium, and finally the 
exterior problem for the wave equation in a random 
medium. The latter portion of Keller's analysis is 
based on perturbation theory and contains many 
new and important results for the field and its first 
and second moments. The principal result is that 
to terms of the 3d order in the perturbation parameter 
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€ the mean value of the field propagates as if in a 
medium of effective complex refractive index defined 
by 

n *=\ 1+€ 2 <ju 2 >(1-2;& f (e 2ikr -l)N(r)dr 

V. Jo 

where €m(x) is the random component of the refrac- 
tive index n(x) and iV(|x— x'|) is the correlation 
function </x(x)/z(x')>/<m 2 >. 

Another treatment of the problem of radiation in a 
randomly inhomogeneous medium has been given 
by W. C. Meecham [1961]. Through appropriate 
use of Green's functions Meecham was able to deter- 
mine the field generated by a point source in a 
medium having small point to point fluctuations of 
dielectric constant. Energy is found to be trans- 
ferred from the average field to the fluctuation field 
at a rate dependent upon the low frequency com- 
ponents of the random fluctuation of the medium, 
and it is shown that the phase velocity of the average 
solution is decreased. The method is applied to the 
case of a medium fluctuation governed by a Gaussian 
correlation function. The average solution thus 
formed is in accord with the result of Foldy and Lax 
for the random configuration problem. For the one- 
dimensional case the formal procedures of the Foldy- 
Lax method have been put on a rigorous basis by 
J. Bazer [1959], who derived sufficient conditions 
for the approximation of the solution of the random 
configuration problem by the corresponding solu- 
tion for a continuous distribution of scattering 
material. These approaches are also closely related 
to the work of Twersky [1962], discussed elsewhere 
in the Commission 6 report under the head of mul- 
tiple scattering. 

A rigorous probabilistic treatment of the electro- 
magnetic field in a randomly inhomogeneous medium 
has been given by Hoffman [1959], based upon 
proper spectral representations for the stochastic 
process representing the refractive index. A gen- 
eralized 3-dimensional Riccati equation is obtained 
and used to relate spectral representations of field 
and gradient of refractive index. On the other hand 
the reduced inhomogeneous wave equation is shown 
to be equivalent to an inhomogeneous Fredholm 
integral equation. The latter can be solved by a 
mean-square convergent Neumann series for the 
case of a Gaussian n 2 (x) provided the Neumann 
series for the nonrandom medium converges in the 
ordinary way. Thus multiple scattering of all orders 
can be taken into account in this formulation. 

A theoretical treatment of the effect of medium 
irregularities of the order of a wavelength in size 
upon the field of an oscillating dipole has been given 
by E. C. Barrows [1963]. Barrows' technique con- 
sists in first expressing the spatial fluctuations in 
permittivity as a Fourier integral, introducing the 
result into the Green's dyadic solution of the vector 
wave equation, and evaluating the result to within 
the Born approximation. The probabilistic nature 
of the Fourier integral representation is not men- 
tioned, although an attempt is made to introduce 
"randomness" by means of a time average correla- 
tion matrix for the electric field itself. 



3. Ray Paths in Random Media 

A substantial portion of Keller's paper mentioned 
above [Keller, 1962] is devoted to determination of 
ray paths in a medium containing small random 
inliomogeneities. After determining the perturba- 
tion equations which govern such a ray path, Keller 
obtains formulas (valid for not too large path lengths) 
for the average ray position, equivalent path length, 
mean-square distance between the endpoints of a ray, 
the mean-square transverse displacement, and the 
ray diffusion coefficient. He then treats the propa- 
gation of rays as a Markov process starting from the 
appropriate Fokker-Planck equation and obtains 
results which agree with those of the perturbation 
equation approach within their mutual ranges of 
validity. Phase and amplitude fluctuations are then 
studied by means of normal congruences of rays 
and the method of characteristics, and the respective 
means and variances are calculated. 

W. S. Anient [1960] has analyzed the effect of 
multiple scattering upon tropospheric forward scatter 
signals in terms of ray tracing through successive 
" blobs." An analog computer specifically designed 
for ray-tracing was employed and perturbations 
consisting of delta-function voltages of random am- 
plitudes and zero means were introduced with random 
spacing. Anient proposed several problems of fun- 
damental importance that bear directly, in the 
further treatment of the multiple scattering ap- 
proach, upon the neglected phase coherence aspect 
of u blob" scattering. 

4. Transport Equations for Field Quantities in a 
Random Medium 

Although electromagnetic waves must obey a 
vector wave equation, certain quantities associated 
with a randomly scattered field, notably the prob- 
ability density function for the ray position after n 
scatterings and the spectral density function for the 
scattered energy, must obey a diffusion equation. 
This approach assumes that the scattered electro- 
magnetic field constitutes a stationary (strict sense) 
Markov process. The probability density function 
must therefore satisfy a Fokker-Planck differential 
equation (or alternatively, in integral form, the 
Chapman-Kolmogorov equation). In the United 
States this approach to the problem of scattering by 
a random continuum has been carried forward by 
D. S. Bugnolo [1960a, 1960b, 1961a]. The basic 
theory of Bugnolo 's approach is set forth mainly in 
the first of the above references, in which he develops 
a transport equation for the spectral density function 
of a multiply scattered electromagnetic field and gives 
a detailed solution for the case of dielectric noise. 
The result is applied by Bugnolo [1960a and 1960b] 
to forward scatter in the troposphere and to laying 
down a statistical criterion for the validity of the 
usual single-scattering hypothesis. The spectral 
width increases with distance in the transport for- 
mulation, whereas this realistic feature is lacking in 
the usual single scattering theories of forward scatter. 
Bugnolo [1961a] has also employed a transport 
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equation for the photon density function to analyze 
the effects of ionospheric multiple scattering upon 
radio star scintillation. 

The method of "invariant embedding" [Bellman 
and Kalaba, 1960; 1961 ; Ueno, 1 96 1 1 constitutes another 
aspect of the approach via transport equations. In 
this method functional equations are obtained for 
the reflection and transmission coefficients of a plane 
wave incident normally upon one or more slabs 
whose refractive indexes are random functions of 
distance from the interface. The functional equa- 
tions in the single slab case turn out to be the well- 
known Riccati equations for the reflection and trans- 
mission coefficients of an inhomogeneous transmission 
line. Reflection from a stratified slab, each stratum 
having a different random wave number, has also 
been considered, and recurrence relations for the 
transmission and reflection coefficients obtained. 
Some attempt has been made to give a probabilistic 
cast to the results, but no real advance has yet been 
accomplished for wave propagation problems. The 
method has enjoyed notable success with respect to 
such transport phenomena as neutron scattering and 
radiative transfer in stellar atmospheres, but seems 
inherently incapable of taking into account diffrac- 
tion effects or non-plane-parallel geometries. 

Significant results with respect to the one-dimen- 
sional scattering problem for a random stack of 
dielectric slabs have been achieved by Kay and 
Silverman [1958]. These authors obtain a Neumann 
series for the mean-square transmission and reflection 
coefficients resulting from multiple scattering. An 
upper bound for the strength of scattering of any 
order from the random medium is used to show that 
these Neumann series converge much more rapidly, 
for a large number of slabs, than would be expected 
on the basis of non-random scattering. On the other 
hand the randomness of the scattering medium does 
not effectively reduce multiple scattering when the 
order of the scattering exceeds the number of slabs. 

5. Turbulent Scatter 

Interest in the problem of electromagnetic radia- 
tion in a randomly inhomogeneous continuum has 
been motivated principally by the phenomenon of 
turbulent forward scatter, it is not our intention 
to review here the theory of turbulent forward scat- 
ter. This subject properly belongs to the domain 
of Commissions 2 and 3, and several excellent review 
articles are in existence [Wheelon, 1959; Bremmer, 
1963]. However, those aspects of turbulent forward 
scatter theory which are relevant to areas of interest 
to Commission 6 will be briefly surveyed here. 

Three theoretical approaches to the problem of 
the turbulent forward scatter have been pursued: 
(i) extensions of the Booker-Gordon single scatter- 
ing theory, (ii) partial reflection from atmospheric 
"platelets", and (hi) diffusion theories for the mul- 
tiple scattering of signals. The latter have been dis- 
cussed under sec. 4 above, and the only U.S. con- 
tributions to the second category have been made 
by Friis, Crawford, and Hogg [Bremmer, 1963], prior 
to the triennium covered in this report. There re- 



mains category (i): generalizations of the classical 
single scattering theory. Contributions in this area 
by Barrows [1963] and Wheelon [1959] have already 
been noted, in addition, Bowhill [1961a] has studied 
the case of a medium with inhomogeneities having 
different scales along three space axes. In this event 
the effect of the medium is not equivalent to a succes- 
sion of thin phase-changing screens, which are inde- 
pendent, probabilistically speaking, of one another. 
The diffractive changes as the wave passes from one 
inhomogeneity to one arbitrarily nearby correlated 
with it renders invalid such an extension of the theory 
of ernerging angular power spectrum. In a compan- 
ion paper Bowhill [1961b] discusses the effect of a 
"random shallow phase screen" upon the transmis- 
sion of radio waves. Bowh ill's analysis applies only 
to a stationary Gaussian process. The results are 
expressed in terms of spatial spectra and correlation 
functions for the emergent wave. The effects of 
anisotropy of the diffracting screen and oblique inci- 
dence of the waves are also included in the analysis. 

Yeh [1962] has investigated the propagation of 
spherical waves through a medium containing anis- 
tropic (elongated) inhomogeneities. The correlation 
function of the refractive index is taken to be a 
Gaussian function of position with ellipsoidal sym- 
metry. A perturbation analysis applied to the 
scalar reduced wave equat ion then leads to first order 
expressions for the mean-square deviations and 
correlation functions of the phase and logarithmic 
amplitude. 

A sort of turbulent scatter which depends upon the 
velocity distribution functions of the electrons and 
ions in a plasma has received much consideration 
recently in connection with "incoherent scatter" of 
high-power radio signals in the ionosphere. A dis- 
cussion of American contributions to this field may 
be found in the section of the Commission 6 report 
devoted to plasma propagation. 

6. Miscellaneous Theoretical and Statistical 
Complements 

In a comprehensive discussion of the subject of 
scattering by a randomly inhomogeneous medium 
Silver [1963] brought out many of the difficulties 
which beset an adequate general theory. To gen- 
eralize single particle scattering to a random config- 
uration, a solution of the many body problem must 
be available and certain questions having to do with 
the ergodic hypothesis must be settled. If the 
medium is turbulent, it may well not be in the 
thermodynamic equilibrium, and this feature poses 
an additional complication for theories based on 
statistical mechanics. If on the other hand a ran- 
domly fluctuating continuum is assumed, then an 
approach to a generally valid theory of constitutive 
parameters is required, e.g., in such media as stellar 
atmospheres. Silver advocates an approach to prob- 
lems of this second type through the polarization 
vector P and the probabilistic properties of the 
Fourier amplitudes representing the distribution of 
scattering elements. It is also noted that the 
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assumption of monochromatic time dependence 
usually made requires that no Doppler diffusion of 
frequency take place. An ergodic theorem, relating 
temporal behavior to the ensemble of stationary 
configurations, thus seems required. Such a theorem 
would probably involve a hypothesis that time 
variations of the constitutive parameters are of a 
much lower order than the propagation time of the 
wave. Alternatively, from a theoretical standpoint, 
one could restrict attention to independent realiza- 
tions, that is, thinking of the medium as a continu- 
ously swirling fluid, one would "freeze" the medium 
at times far enough apart so that the resulting con- 
figurations (realizations) of refractive index are 
independent in the probability sense. 

The application of the ergodic theorem to problems 
of wave propagation through a random time-varying 
medium has been considered by T. J. Skinner [1961]. 
He showed that the time-varying medium can be 
replaced in the solution of the scalar problem by a 
time independent ensemble provided at the boundary 
of the time varying medium the scalar field V(x, t) = 
Afat) exp {i<t)(x,t)- J r2Tiat} is such that 

(i) A(x,t) and <l>(x,t) are ergodic; 

(ii) The absolute value of any frequency compo- 
nent of either A(x, t) sin </>(x, t) or A(x, t) cos 
0(x, t) is less than a, the mean frequency of 
the source; 

(hi) V(x, t) is quasi-monochromatic. 

D. S. Bugnolo [1961b] has given a heuristic dis- 
cussion of the question of stationariness (in the prob- 
ability sense) when the dielectric constant is a 
stochastic process in both space and time. He con- 
cludes in the tropospheric case that the permittivity 
must of necessity be a nonstationary process and 
proposes a redefinition of the expected value of the 
permittivity which would make the corresponding 
homogeneous vector wave equation consistent with 
the purely deterministic problem for an inhomogene- 
ous atmosphere. 

C. I. Beard [1961] in a clever experimental simula- 
tion of the idealized conditions of statistical me- 
chanics has shown that in the "mid-field" region the 
phase quadrature components of the incoherent 
scattered field are not Rayleigh distributed. 

Samuels and Eringen [1959/60] have studied the 
behavior of the nt\\ order linear differential equation 
with random coefficients and find a phenomenon of 
mean-square instabilit}^ under certain conditions. 
The problem of the dynamics of first order linear and 
nonlinear oscillators with random coefficient has 
been exhaustively studied by Kraichnan [1961] by 
means of an auxiliary set of random coupling 
coefficients. Bourret [1961] has investigated the 
angular response of a linear array to a scalar signal 
from a point source when the medium is slightly 
randomly inhomogeneous. 

Significant statistical contributions (in the sense of 
statistical inference) to the propagation of electro- 
magnetic waves in random media have been made by 
S. S. Siddiqui [1960, 1961, 1962]. M. J. Beran [I960] 



has derived the equations governing the propagation 
of the mth order correlation function for a scalar wave 
field, and studied the propagation of the second order 
correlation function between parallel planes in free 
space. 
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2.3. Partially Coherent Electromagnetic Fields 

Francis J. Zucker * 



This summary report begins where the previous 
URSI report on coherence theory [Parrent, 1960] left 
off. Activity during the past three years has been 
brisk. A conference held in 1960 [University of 
Rochester, 1961; O'Neill and Bradley, 1961] paid 
much attention to partial coherence, and the Pro- 
ceedings of the 1962 URSI Copenhagen Symposium 
on Electromagnetic Theory and Antennas (Pergamon 
Press, Oxford, 1963) will contain six papers in the 
field. The proceedings of two 1963 meetings, the 
Third International Conference on Quantum Elec- 
tronics (held in Paris) and the Symposium on Optical 
Masers (held at the Polytechnic Institute of Brook- 
lyn) should be watched for articles on laser coherence 
and on the connection between coherence and 
quantum mechanics. 

The basic reference is still Wolf [Born and Wolf, 
1959]. An introduction to the present state of the 
art was written by Zucker [1963], and two compre- 
hensive surveys [Mandel, 1963a; Gamo, 1963b] cover 
recent developments. A book on statistical optics 
[O'Neill, 1963] devotes two chapters to coherence 
topics, and another book, now in press [Beran and 
Parrent, 1963b] promises a detailed treatment of the 
entire subject. 

We first consider applications of partial coherence 
theory to problems in microwaves and optics. 

Aperture diffraction. The diffraction patterns of 
infinite slots [Parrent and Skinner, 1961] and circular 
apertures [Shore, 1963] have been calculated for a 
continuous range of aperture size to coherence inter- 
val ratios; as might be expected, the pattern ap- 
proaches that of a dipole as the coherence interval 
becomes small compared to the aperture diameter. 
Schell [1961] pointed out that the partially coherent 
power pattern is given by the convolution of the 
coherent pattern with the Fourier transform of the 
degree of coherence, which suggests rapid graphical 
methods for antenna and aperture pattern determina- 
tions. The polarization aspects of partially coherent 
diffraction were examined by Karczewski and Wolf 
[1963], with the result that approximate evaluations 
known to give almost identical results for mono- 
chromatic radiation turned out to differ considerably 
in the polychromatic case. 

Propagation in random media. By the ergodic 
hypothesis, it should be possible to establish an 
equivalence between the aperture pattern in a clear 
atmosphere for a partially coherent signal, and the 
aperture pattern in a turbid atmosphere for a coher- 
ent signal, provided the statistics of the partially 
coherent signal and of the turbid atmosphere are 
identical; the conditions under which this equivalence 
holds, however, have not been rigorously established. 
A formal attack on the propagation of partially 
coherent radiation through a turbid atmosphere (i.e., 
randomness in signal as well as in medium) was made 
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by Parrent, Shore, and Skinner [1962], and specific 
applications by Skinner are in the offing. 

Partial polarization. The original work in optics 
on the relation between degree of polarization, Stokes 
vector, and coherency matiix [Born and Wolf, 1959; 
Parrent and Roman, 1960] has been capably re- 
viewed by Marathay [O'Neill, 1963], and adapted 
for antenna usage by Ko [1961a, b, c; 1962]. The 
response of an elliptically polarized antenna to quasi- 
monochromatic, partially polarized radio waves from 
an extended source has been worked out in detail. 
Ko clearly shows the advantages of this language in 
the design and evaluation of radioastronomical 
polarization measurements. The effect of partial 
polarization on the effective antenna aperture has 
also been considered [Tai, 1961; Ko, 1962]. 

Imaging. The optical theory of image formation 
has been generalized to include the effects of partial 
coherence [Parrent, 1961], which in turn has sug- 
gested a re-examination of the concept of resolution 
[Parrent and Rojak, 1963]. As a most useful appli- 
cation in the microwave field, it has been shown 
[Drane and Parrent, 1962] that data-processing 
antennas of the Drane-Covington type map extended 
incoherent sources (such as the radio sky) linearly 
in power, even though they employ nonlinear ele- 
ments. Because of the finite aperture size, a 
sampling theorem originally due to Gabor [Miya- 
moto 1960, 1961] can be applied to the object and 
image, leading to a formulation of imaging in terms 
of finite matrices or equivalently, because of the 
finite number of degrees of freedom, to a formulation 
in terms of entropy transformations [Gamo, 1960; 
O'Neill and Asakura, 1961; O'Neill, 1963; Gamo 
1963b]. 

The coherence aspects of two important non- 
Gaussian radiators have been investigated, the 
blackbody and the maser (and laser). An initial 
discussion of the blackbody by Bourret [1960] was 
given a broader coherence-language base by Kano 
and Wolf [1962] and Mehta [1963], and then related 
to its quantum mechanical background by Sarfatt 
[1963]. That the spatial coherence of maser light 
can be explained in terms of the general property of 
incoherent radiation becoming progressively more 
coherent as it passes through a periodic medium was 
discovered by Wolf [1963]. That maser radiation 
differs widely from thermal light is well known 
[Mandel, 1961a; Smith and Williams, 1962]; it is 
therefore a significant advance that at the Third 
International Conference on Quantum Electronics 
(Paris, 1963), Mandel was able to present a prob- 
ability distribution for the radiation of a well- 
stabilized single-mode maser, from which the 
comolete coherence properties should now be deducible. 

The concept of coherence itself has been clarified 
and extended in several directions. 

Clarifications of coherence. Skinner [1961] showed 
that arbitrarily narrow-band (but not monochro- 
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matic) radiation can be very incoherent, and Mandel 
and Wolf [1961b] demonstrated the converse, that 
to be coherent, radiation does not necessarily have 
to be very narrow-band. If the coherence function 
can be factorized into a purely space-dependent and 
a purely time-dependent term, Mande] [1961b] refers 
to the radiation as "spectrally pure"; he shows that 
light from an extended incoherent source has this 
property. The limiting cases of complete coherence 
and incoherence have been carefully analyzed, with 
the result that neither condition is physically 
realizable [Parrent, 1963; Beran and Parrent, 1963a, 
1963b]. In preparation for the treatment of non- 
v Gaussian sources, the concepts of bandwidth and 
coherence time have been carefully re-examined 
[Mandel and Wolf, 1962]. 

Generalizations oj coherence. To describe the 
coherence properties of electromagnetic fields, the 
coherence function of scalar fields (or the coherency 
matrix of partially polarized plane waves) must 
now be replaced by three-by-three electric and 
magnetic coherency tensors [Roman and Wolf, 
1960; Roman, 1961b; Ivarczewski, 1963]. The phys- 
ical meaning of some of the differential equations 
connecting these tensors is still obscure. Source 
terms can be included in the field equations [Beran 
and Parrent, 1962; Roman, L961a], but attempts at 
extending the equations to nonstationary ensembles 
[Kano, 1962] have not yet succeeded. 

Transient interference. The limitation of the co- 
herence concept to infinite time averages has inhib- 
ited its application to all cases in which time-variable 
interference effects might occur. The discussion on 
this point has been brisk [Golay, 1961; Bracewell, 
1962; Neugebauer, 1962; Ready, 1962]. Mandel 
[1962b] contributed the fundamental observation 
that in the case of Gaussian sources, transient inter- 
ference can be handled if the concept of "photon 
degeneracy" (the number of photons falling on a 
coherence area in the coherence time) is added to 
that of coherence, while for non-Gaussian sources, 
higher order correlation functions will have to be 
taken into account as well. 

So far, only second-order correlation (i.e., coher- 
ence) effects have been discussed; we now turn to 
fourth-order effects, i.e., intensity correlations. It 
has been known for some time (H anbury Brown 
and Twiss experiments) that intensity correlations 
can yield information about the spatial and spectral 
intensity distributions of sources; Wolf showed that 
they also yield polarization information [Wolf, 1960; 
Mandel and Wolf, 1961a; Mandel, 1963c]. Com- 
pared with the traditional (Michelson) second-order 
method, intensity correlation has the advantage of 
relative independence from atmospheric scintilla- 
tions, but the disadvantage of yielding only the 
magnitude, not the phase of the coherence function; 
unless the intensity, spectral, or polarization dis- 
tributions are known to be symmetric, they cannot 
therefore be fully determined. Gamo has devised 
an ingenious optical method for the direct measure- 
ment of the phase of the coherence function by 
injecting a coherent background; the method works 
best with intense and narrowband sources [Gamo, 



1961a, 1961b, 1963a]. More attention has to date 
been paid to this method in optics [Mandel, 1963a] 
than in microwaves; the physics of optical square-law 
detectors and photon-coincidence counters has there- 
fore been carefully investigated [Mandel, 1960; 
Mandel and Wolf, 1961a; Mandel 1963a, 1963c]. 

In intensity-correlation spectroscopy [Mandel, 
1963b], Wolf showed that analytic properties of the 
coherence function often lead to a full determination 
of the spectral profile, even when Gamo's phase 
determination is not available [Wolf, 1962]. For- 
rester [1961] proposed an optical superheterodyne 
technique for improving the signal-to-noise ratio, 
with which Javan et al. [1961] then measured laser 
spectra. Spectral modulation in superposed, co- 
herently modulated light beams, first discussed by 
Alford and Gold, has been extended to coherent but 
unmodulated beams by Givens [1961, 1962 1, and 
was shown by Mandel [1962a] to be a fourth-order 
correlation effect, like that of Hanbury Brown and 
Twiss. 

For a complete stochastic description of electro- 
magnetic radiation, the total probability distribu- 
tion must of course be specified. In the Gaussian 
case, the probability distribution depends only on 
the phase coherence and on the degree of polarization 
| Mandel, 1963c], but in general all higher-order 
probability densities, and thus correlations, must be 
known [Beran and Parrent, 1963b]. Wolf showed 
at the Symposium on Optical Masers (Polytechnic 
Institute of Brooklyn, 1963) that a quantum-mechan- 
ical transcription by means of the correspondence 
principle of the higher-order correlation leads to the 
quantum-theoretical coherence concept already for- 
mulated by (dauber [1963; Mandel and Wolf, 1963]. 
The question as to the relation between coherence 
theory and quantum mechanics, raised some time 
ago by Purcell, Dicke, and Fano (see for example 
University of Rochester [1961]), was carefully 
examined by Sudarsham at the Brooklyn Sympo- 
sium; it appears that, by working with the analytic 
signal and admitting negative as well as positive 
probability distributions, one obtains a stochastic 
theory of coherence which, at least for all linear 
effects, is isomorphic with the quantized description 
of radiation in terms of the density matrix (see also 
Sudarshan [1963]). 

References 
Beran, M. and G. B. Parrent (1962), On the equations 

governing the coherence of the radiation emitted by 

random charge and current distributions, J. Opt. Soc. Am. 

52, 98. 
Beran, M. and G. B. Parrent (1963a), The mutual coherence 

of incoherent radiation, Nuovo Cimento 27, 1049. 
Beran, M. and G. B. Parrent (1936b), Theory of partial 

coherence (Prentice-Hall, New Jersey). 
Born, M. and E. Wolf (1959), Principles of Optics, Ch. 10 

(Pergamon Press, London and New York). 
Bourret, R. C. (1960), Coherence properties of blackbody 

radiation, Nuovo Cimento 18, 347. 
Bracewell, R. N. (1962), Defining the coherence of a signal, 

Proc. IRK 50, 214. 
Drane, C. J. and G. B. Parrent (1962), On the mapping of 

extended sources with nonlinear correlation antennas, IRE 

Trans., Ant. Prop. AP-10, 126. 
Fano, U. (1961), Quantum theory of interference effects in 

the mixing of light from phase-independent sources, An. J. 

Phys. 29, 539. 



461 



Forrester, A. T. (1961), Photoelectric mixing as a spectro- 
scopic tool, J. Opt. Soc. Am. 51, 253. 

Gamo, H. (1960), An aspect of information theory in optics, 
[RE L960 International Convention Record, pt, 4, 189. 

Gamo, H. (1961a), On the intensity interferometer with co- 
herence background, Advances in Quantum Electronics, 
252 (Columbia University Press, New York). 

Gamo, H. (1961b), An intensity interferometer for partially 
polarized and partially coherent beams, J. Opt. Soc. Am. 
51, 1474. 

Gamo, H. (1963a), Phase determination of coherence func- 
tions by the intensity interferometer, Copenhagen Sympo- 
sium on Electromagnetic Theory and Antennas, June 1962, 
ed. E. C. Jordan, pt. 2, 801-810 (Pergamon Press, Oxford). 

Gamo, H. (1963b), A matrix theory of partial coherence, 
Progress in Optics 3 (North-Holland Publishing Co., 
Amsterdam, in press). 

Givens, M. P. (1962), Photoelectric detection of interference 
between two light beams having a large path difference, 
J. Opt. Soc. Am. 51, 1030 (1961) ; also J. Opt. Soc. Am. 52, 225. 

Glauber, R. (1963), Photon correlations, Phvs. Rev. Letters 
10, 84. 

Golay, M. J. E. (1961), Note on coherence vs narrow banded - 
ness in regenerative oscillators, masers, lasers, etc., Proc. 
IRE 49, 958. 

Javan, A., W. R. Bennett, Jr., and D. R. Herriott (1961), 
Population inversion and continuous optical maser oscilla- 
tion in a gas discharge containing a He-Ne mixture, Phys. 
Rev. Letters 6, 106-110. 

Kano, Y. (1962), Equations governing the propagation of 
second-order correlations in non-stationary electromagnetic 
fields, Nuovo Cimento 23, 328. 

Kano, Y., and E. Wolf (1962), Temporal coherence of black 
body radiation, Proc. Phys. Soc. (London) 80, 1273. 

Karczewski, B. (1963), Degree of coherence of the electro- 
magnetic field, Physics Letters 5, 191. 

Karczewski, B. and E. Wolf (1963), Polarization properties 
of the electromagnetic field diffracted from an aperture, 
Copenhagen Svmposium on Electromagnetic Theory and 
Antennas, June 1962, ed. E. C. Jordan, Pt. 2, 797-799 
(Pergamon Press, Oxford). 

Ko, H. C. (1961a), Theoretical techniques for handling 
partiallv polarized radio waves, with special reference to 
antennas, Proc. IRE 49, 1446. 

Ko, H. C. (1961b), The use of the statistical matrix and the 
Stokes vector in formulating the effective aperture of an- 
tennas, IRE Trans. Ant. Prop. AP-9, 581. 

Ko, H.C. (1961c), On the response of a radio antenna to 
complex radio waves, Proc. Nat. Electron. Conf. 17, 500. 

Ko, H. C. (1962), On the reception of quasi-monochromatic, 
partiallv polarized radio waves, Proc. IRE 50, 1950. 

Mandel, L. (1960), Detection of interference effects with 
incoherent light beams by photon counting, J. Opt. Soc. 
Am. 50, 1131. 

Mandel, L. (1961a), Photon degeneracy in light from optical 
maser and other sources, J. Opt. Soc. Am. 51, 797. 

Mandel, L. (1961b), Concept of cross-spectral purity in 
coherence theory, J. Opt. Soc. Am. 51, 1342. 

Mandel, L. (1962a), Interference and the Alford and Gold 
effect, J. Opt. Soc. Am. 52, 1335. 

Mandel, L. (1962b), Transient coherence in optics, J. Opt. 
Soc. Am. 52, 1407. 

Mandel, L. (1963), Fluctuations of light beams, Progress in 
Optics 2 (North-Holland Publishing Co., Amsterdam). 

Mandel, L. (1963b), The determination of spectral density 
distributions from measurements of radiation fluctuations, 
Copenhagen Svmposium on Electromagnetic Theory and 
Antennas, June 1962, ed. E. C. Jordan, Pt. 2, 811-817 
(Pergamon Press, Oxford). 

Mandel, L. (1963c), Intensity fluctuations of partially 
polarized light, Proc. Phys. Soc. (London), in press. 

Mandel, L. and E. Wolf (1961a), Correlation in the fluctu- 
ating outputs from two square-law detectors illuminated by 
light of any state of coherence and polarization, Phys. Rev. 
124, 1696. 

Mandel, L. and E. Wolf (1961b), Some properties of coherent 
light, J. Opt. Soc. Am. 51, 815. 

Mandel, L. and E. Wolf (1962), The measures of bandwidth 
and coherence time in optics, Proc. Phys. Soc. (London) 
80, 894. 



Mandel, L. and E. Wolf (1963), Photon correlations, Phvs. 
Rev. Letters 10, 276. 

Mehta, C. L. (1963), Coherence-time and effective band- 
width of blackbody radiation, Nuovo Cimento 28, 401. 

Mivamoto, K. (1961), On Gabor's expansion theorem, J. Opt. 
Soc. Am. 50, 856 (1960); also J. Opt. Soc. Am. 51, 910. 

Neugebauer, H. E. J. (1962), Coherence time of a maser, 
J. Opt. Soc. Am. 52, 470. 

O'Neill, E. L. (1963), Introduction to Statistical Optics 
(Addison- Wesley, Boston). 

O'Neill, E. L. and T. Asakura (1961), Optical image forma- 
tion in terms of entropy transformations, J. Phys. Soc. 
Japan 16, 301. 

O'Neill, E. L. and L. C. Bradley (1961), Coherence properties 
of electromagnetic radiation, Physics Today 14, No. 6, 28. 

Parrent, G. B. (1960), A bibliographv on coherence theorv, 
J. Research NBS 64D (Radio Prop.) 742. 

Parrent, G. B. (1961), Imaging of extended polychromatic 
sources and generalized transfer functions, J. Op. Soc. Am. 
51, 143. 

Parrent, G. B. (1963), Some fundamental questions of 
coherence theory, Copenhagen Symposium on Electro- 
magnetic Theory and Antennas, June 1962, ed. E. C. 
Jordan Pt. 2, 785-786 (Pergamon Press, Oxford). 

Parrent, G. B., and F. Rojak (1963), Effects of coherence on 
resolution, J. Opt. Soc. Am. 53 (in press). 

Parrent, G. B. and P. Roman (1960), On the matrix formula- 
tion of the theory of partial polarization in terms of ob- 
servables, Nuovo Cimento 15, 370. 

Parrent, G. B., R. A. Shore, and T. J. Skinner (1962), On 
the mutual coherence function in an inhomogeneous 
medium, J. Math. Phys. 3, 678. 

Parrent, G. B., and T. J. Skinner (1961), Diffraction of 
partially coherent light by a plane aperture, Optica Acta 8, 
93. 

Ready, J. F. (1962), Time coherence in rubv lasers, Proc. 
IRE 50, 1695. 

Roman, P., (1961a), Correlation theory of stationary electro- 
magnetic fields. Part III — The presence of random sources, 
Nuovo Cimento 20, 759. 

Roman, P. (1961b), Correlation theory of stationary electro- 
magnetic fields. Part IV — Second order conservation laws, 
Nuovo Cimento 22, 1005. 

Roman, P. and E. Wolf (1960), Correlation theory of sta- 
tionary electromagnetic fields, part I: the basic field 
equations; part II: conservation laws, Nuovo Cimento 17, 
462 and 477. 

Sarfatt, J. (1963), Quantum mechanical correlation theory of 
electromagnetic fields, Nuovo Cimento 21, 1119. 

Schell, A. C. (1961), Doctoral dissertation, Massachusetts 
Institute of Technology. 

Shore, R. A. (1963), Partially coherent diffraction by a circular 
aperture, Copenhagen Symposium on Electromagnetic 
Theory and Antennas, June 1962, ed. E. C. Jordan, pt. 2, 
787-795 (Pergamon Press, Oxford). 

Skinner, T. J. (1961), Incoherent source with an arbitrarily 
narrow power spectrum, J. Opt. Soc. Am. 51, 909. 

Smith, A. W., and G. W. Williams (1962), On the detection of 
maser signlas bv photoelectric mixing, J. Opt. Soc. Am. 52, 
337. 

Sudarshan, E. C. G. (1963), Equivalence of semiclassical and 
quantum mechanical descriptions of statistical light beams, 
Phys. Rev. Letters 10, 277. 

Tai, C. T. (1961), On the definition of the effective aperture 
of antennas, IRE Trans. Ant. Prop. AP-9, 224. 

University of Rochester (1961), Report on Conference on 
Coherence Properties of Electromagnetic Radiation, 
AFOSR-583, Report No. 5 under Contract AF49(638)-602, 
University of Rochester, N.Y. 

Wolf, E. (1960), Correlation between photons in partially 
polarized light beams, Proc. Phys. Soc. (London) 76, 424. 

Wolf, E. (1962), Is a complete determination of the energy 
spectrum of light possible from measurements of the degree 
of coherence?, Proc. Phys. Soc. (London) 80, 1269. 

Wolf, E. (1963), Spatial coherence of resonant modes in a 
maser interferometer, Physics Letters 3, 166. 

Zucker, F. J. (1963), Introduction to partially coherent 
electromagnetic waves, Proc. 1962 Copenhagen Symposium 
on Electromagnetic Theory and Antennas (Pergamon 
Press, Oxford). 



462 



Part 3. Radiation 



3.1. Electromagnetic Fields in Lossy Media 



James R. Wait 1 



In recent years the propagation of electromagnetic 
waves in dissipative or lossy media lias received con- 
siderable attention. This work has been motivated 
by a need to communicate between underground or 
underwater terminals. The general problem has 
been discussed by Ghose [1960a] and Hansen [1963] 
in two valuable papers. 

One proposed mode of propagation is based on the 
existence of a zone with very low conductivity within 
the crust of the earth [Wait, 1954; Wheeler, 1961; 
de Bettencourt and Carolan, 1963]. According to 
G. V. Keller [1963], it is reasonable to expect that 
s, such a zone exists, since the pressures at depth tend 
to reduce the pore space in which water may exist. 

Unfortunately, no measurements have been re- 
ported in the literature which confirm the existence 
of low-conductivity layers in the earth's crust. Some 
of the first experiments indicating low attenuation 
[Barrett, 1949, 1952 ] have been open to serious 
criticism [Pritchett, 1952; Wait, 1954]. It was ac- 
tually suggested by Barrett that conventional elec- 
tromagnetic theory was not valid for geological 
conductors. On the other hand, Pritchett's experi- 
ments were carried out under idea] conditions, and 
it was assured that the signal traversed the medium 
and was not conducted along pipes, wires, and other 
buried conductors; his measured attenuations were 
of the order of 1 db per foot at 1 Mc/s and are in good 
accord with plane wave theory. The more recent 
experiments of Saran and Held [19(H)] also confirm 
the classical exponential depth attenuation predicted 
by simple skin-depth theory; they lowered a rope- 
supported receiving antenna into a fresh-water lake 
about 100 miles from a VLF transmitter. Other 
related experiments have been carried out by Harmon 
[1961], de Bettencourt and Frazier [1963], and Ghose 
[1960b]. 

Scale-model measurements of the fields of a 
horizontal magnetic dipole immersed in a conducting 
fluid were carried out by Kraichman [1960], who 
found good agreement between theory and experi- 
ment. Similar measurements have been carried out 
recently by Blair [1963], employing horizontal 
electric dipoles immersed in the solution. 

Theoretical work on fields within lossy media is 
very extensive. In this category we may include 
ground wave propagation with the antenna above 
or below the air /ground interface. A number of 
papers dealing with this topic have appeared [Gerks, 
1962; Wait/ 1961a; Walters and Johler, 1962; 
Wait and Walters, 1963]. In some of these, atten- 
tion has been given to mixed path conditions such 
as occur in propagation from land to sea. 



i National Bureau of Standards, Boulder, Colo. 



Theoretical treatments dealing with dipoles sub- 
merged in a conducting half-space have been par- 
ticularly numerous. Wait [1961b] has given a 
unified derivation for all distance regions, the termi- 
nal antennas being above or below the interface. 
The derivation made direct use of his earlier pub- 
lished work in the period L951-1954. Moore and 
Blair [1961] have treated the problem as an extension 
of Sommerfeld's classical solutions. Their results, 
which are based partly on Moore's [1951] thesis, 
cover near, intermediate, and far zones but omit the 
quasi-stal ic region. 

Biggs [1962] has considered the radiation field in 
air from a submerged horizontal antenna. The 
final results appear to be identical to those of Norton 
[1937], apart from the expected exponential depth 
factor. Ghose [1960a] has also considered the buried 
horizontal dipole antenna and, at the same time, he 
has given some useful design information. Two 
papers by Durrani [1962 a, and b] treat essentially 
the same problem and contain some new informa- 
tion. The fields in a three-part lossy media have 
been considered by Anderson [1962]; this is the 
seawater-air-ionosphere problem with the source 
dipole located in the seawater medium. 

It is expected that a very comprehensive book on 
the subject of dipoles in the presence of a, half-space 
will be published soon by Professor A. Bafios of the 
University of California in Los Angeles. This will 
be based in part on his earlier (1953) unpublished 
reports coauthored with el. P. Wesley. In fact, 
most of the published work mentioned above is 
contained implicitly in these two monumental 
reports. 

The radiation from pulse-excited dipoles in 
homogeneous conducting media has received some 
attention. Solutions of this type for unbounded 
media were given by Anderson and Moore [1960], 
Galejs [I960], Zisk [1960], and Burrows [1962]. 
The corresponding transient solutions for dipoles 
immersed in a half-space were carried out in an 
extensive paper by Wait [1960], who showed how 
pulse waveforms are distorted when they propagate 
into the half-space. Some extensions were pub- 
lished by Mijnarends [1962], who treated dipoles 
excited by step-function currents. In the results 
given by Anderson and Moore, Galejs, Zisk, and 
Burrows mentioned above, displacement currents 
were neglected. Solutions taking these into account 
have been given by Wait [1960] and Burrows [1963]. 
Their effect is important when dealing with poorly 
conducting media and at small times in the transient 
response. 

Propagation of groundwave pulses over the surface 
of the earth was reviewed by Johler [1962]. Certain 
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ejects of early published results on this subject 
were discussed by Wait [1962al in a short note. 

Two recent survey papers which contain material 
on fields in lossy media were published by Logan 
[1962] and Wait J1962bl. Both of these papers deal 
with methods for computing diffraction fields for 
large convex lossy surfaces. 

In many problems, involving the influence of the 
earth on electromagnetic fields, it is permissible to 
utilize the concept of surface impedance. This 
approach has been utilized extensively in the study 
of magnetotelluric fields [Wait, 1962c; Price, 1962; 
Jackson, Wait, and Walters, 1962]. 

Galejs [1962a] has published an interesting paper 
dealing with scattering from a conducting sphere 
which is embedded in a conducting half-space. The 
induced electric and magnetic dipole moments depend 
on the incident surface wave and the interactions 
with the interface. The fields which are scattered 
from the sphere into the lossless medium are deter- 
mined from known solutions of this dipole problem 
[Wait, 1961b]. Galejs [1962b] has also considered 
the excitation of slots above a lossy dielectric half- 
space. This problem has application to the design 
of peninsula and island antennas. The latter 
problem has also been considered by Staras [1962]. 

Scattering from lossy coated objects in free space 
has been investigated in a thorough manner by Hiatt, 
Siegel, and Weil [1960a and b]. For electrically 
small spheres, it is shown that the coating can have 
little effect on the Rayleigh cross sections. They 
also found that scattering from large objects had a 
strong forward lobe which was actually enhanced by 
the coating. This was true e\ r en for coatings which 
are "radar absorbers." 

The input resistances of dipoles in the presence of 
a conducting half-space have been calculated by 
Vogler [1963]. Similar results have been given by 
Wait [1962d], who gave his results in graphical form. 
The latter author also considered the extension to 
multiple boundaries and to anisotropic lossy media. 

A survey of electromagnetic waves in stratified 
lossy media has been published by Wait [1962e] in a 
recent text. Other related topics contained in this 
book are phase integral methods for dissipative 
media, modes in waveguides with lossy walls, propa- 
gation along a curved (lossy) surface, connections 
between modes and rays in lossy media, influence of 
curvature in terrestrial waveguides, etc. 
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3.2. Antennas in Lossy Media 

C. T. Tai i 



Since this is the first time that U.S. Commission 
VI has selected the above title for a specific report, 
some early studies will be mentioned, mainly to 
compile a useful bibliography. This report will cover 
principally problems concerning the characteristics 
of antennas in a lossy medium of either infinite or 
finite extent. In regard to the nomenclature, one 
normally would interpret a "lossy medium" either as 
one with ohmic loss due to conduction, such as salt 
water, or as one ionized with collision loss, such as 
the ionosphere. However, we shall not adhere to 
such a strict interpretation; problems dealing with 
antennas in an ideal plasma without collision loss, or 
in a lossless anisotropic medium, will also be 
mentioned. 

During the preparation of this report many authors 
have been kind enough to supply material describing 
their work. Unfortunately, a large amount of the 
data quoted by some authors is contained either in 
unofficial company memoranda or technical reports 
to which this reviewer does not have access. It is 
therefore unavoidable that some valuable informa- 
tion may have been overlooked. 

Early Studies 

The early investigations by Campbell [1923], Car- 
son [1926], and Foster [1931] concerning the effects 
of earth conduction on transmission systems repre- 
sent the typical problems that attracted attention 
during the thirties. Much of this work is well 
documented in the excellent book by E. D. Sunde 
[1948], who also made an outstanding contribution 
in this field. The theoretical technique used by these 
authors follows closely the method originated by 
A. Sommerfeld [1909] in formulating the dipole and 
earth problem. No particular emphasis, however, 
was placed on radiating power and antenna charac- 
teristics in these early studies. The radiating power 
of a dipole in a lossy medium was touched upon by 
Sommerfeld [Frank and von Mises, 1935] in formu- 
lating a reciprocity theorem for two dipoles in a lossy 
medium. The master, however, did not define 
clearly the radiating power in this situation. The 
fact that an infinitesimal dipole in direct contact 
with a lossy medium cannot maintain a finite amount 
of power was later pointed out by Sommerfeld and 
Renner [1942]. In the latter part of the forties, 
several groups in the United States were also engaged 
in research on antennas in lossy media. The main 
interest at that time was stimulated by problems 
connected with submarine antennas. 

The author of this report, under the direction of 
Professor R. W. P. King at Harvard, wrote two 
technical reports on this subject [Tai, 1947, 1949]. 
By investigating the power relationship pertaining to 
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an infinitesimal dipole immersed in a lossy medium r 
he confirmed the conclusions of Sommerfeld and 
Renner, namely that the input power is infinite An 
insulated dipole was then introduced as a more 
realistic model and Sommerfeld' s reciprocity theorem 
was modified accordingly, based upon this model. 
The radiation field of a spherically stratified dipole 
was studied by J. B. Keller and H. B. Keller [1951]. 
J. R. Wait [1957] derived an expression for the radi- 
ation resistance of an insulated circular loop. In 
the case of antennas of finite dimensions, the theory 
of cylindrical antennas based upon the integral 
equation method and the theory of biconical anten- 
nas were formally extended to lossy media [Tai, 1949]. 
Except for the case of an insulated biconical antenna, 
no numerical calculations were made according to 
these formulations owing to a lack, at the time, of 
adequate tables of functions, particularly the sine- 
and cosine-integrals and the spherical Bessel func- 
tions with complex arguments. The behavior of 
submarine antennas was also carefully studied by 
Moore at that time. The results of his investigation 
were contained in his Cornell dissertation [Moore, 
1951]. A part of his work deals with antennas in a 
lossy medium of infinite extent. The major portion 
of his dissertation is concerned with the effectiveness 
of wave propagation across an air-water boundary 
with an antenna placed in the water. His work 
describes the fundamental difference between the 
characteristics of an antenna in air and in a lossy 
medium [Moore, 1962]. In the mid fifties, not many 
activities were reported in this field, except for one 
basic work by H. A. Wheeler [1958] on small antennas 
where the importance of the dissipation in the neigh- 
borhood of an antenna placed in a lossy medium is 
emphasized. In recent years, the advancement of 
satellites brought with it a number of problems re- 
lating to the satellite induced plasma and hence a 
renewed activity in this field. 

Recent Studies 

A paper surveying the recent studies on the topic 
which Moore investigated previously was given by 
Hansen [1962] at the Copenhagen Symposium. 
Hansen introduced several useful merit factors in 
analyzing the efficiency of transmission. The char- 
acteristics of low-frequency subsurface radiating 
structures have been investigated in great detail, 
both theoretically and experimentally, by Guy and 
Hasserjian [1962]. Their theoretical results are pre- 
sented in the form of generalized curves which may 
be used to design an optimum antenna for any given 
set of uncontrollable variables such as ground con- 
ductivity and frequency. 

Extensive research has been conducted by King 
and his associates on cylindrical antennas in a con- 
ducting medium [King, 1962; King and Harrison, 
1960; King, Harrison, and Denton, 1961; Iizuka 
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and King, 1962a, b; Gooch, Harrison, King, and Wu, 
1962]. The experimental work by Iizuka and Kino- 
is particularly interesting as it gives a very com- 
prehensive picture of the current distribution on the 
dipole. The measured data are in fair agree hum it 
with the theoretical results obtained by King and 
Harrison [I960]. Larson [1962] has obtained some 
experimental data on the input impedance of a 
plasma imbedded dipole over a wide range of colli- 
sion and operating frequencies. It is found that the 
theory of King and Harrison for short and half- 
wavelength dipoles is generally followed, and fair 
agreement is obtained with a Langmuir probe on the 
values of electron density. 

The input impedance of a circular loop in a dis- 
sipative medium has been computed by Kraichman 
[1962] and by Chen [1962]. The latter is an exten- 
sion of Storer's work [1955] on the impedance of a 
thin-wire loop antenna in air. Numerical values 
have been compiled by Chen and Tai [1962a] for the 
input impedance of thin biconical antennas in a dis- 
sipative medium. The same two authors [1962b] 
have also computed the radiation patterns of some 
linear arrays in a lossy medium. The relationship 
between the impedance function of a class of an- 
tennas defined in a dissipative medium and in a 
lossless medium is pointed out by Deschamps [1962]. 
For media with moderate losses, lengthy computa- 
tions can be avoided by using a graphical extrap- 
olation based upon the existing data obtained for a 
lossless medium. Polk [1959] lias studied the 
resonance and supergain effects in small ferro- 
magnetically or dielectrically loaded biconical an- 
tennas. The effect of dielectric loading of electric 
dipole antennas has also been investigated by 
Galejs [1962]. The results show that Wheeler's 
general conclusion is not applicable to some cases. 

In a number of technical reports, Katzin and his 
associates [1957, 1958, 1959, 1960] have investi- 
gated the radiation pattern and impedances of several 
radiating systems under the influence of a plasma 
sheath. 

Since the original formulation by Bunkin [1957] 
on the problem of radiation in an anisotropic medium, 
some advancements have been made [Arbel, 1960; 
Kogelnik, 1960; Kuehl, 1961; Mittra and Deschamps, 
1962; Wu, 1962]. The work of Arbel was particularly 
elegant. It supplied many details on the power 
relationship and the radiation patterns of electric 
dipoles in such a medium. Kogelnik and Kuehl 
have evaluated the radiation resistance of electric 
dipoles in a magnetoionic medium for the case of a 
very high gyromagnetic frequency. Mittra and 
Deschamps have studied the near-zone field. Wu 
has supplied some details about the far-zone fields 
due to magnetic dipoles. Very little is yet known 
about the input impedance of an antenna placed in 
such a medium. Cohen [1962] has obtained the 
radiation resistance of a current element in a plasma 
medium which takes into consideration the acoustic 
source as well as the electromagnetic source. 

Upon the completion of this report, we were in- 
formed of an extensive research program conducted 



at Raytheon Company on antennas and fields in 
rock strata. The details are contained in a report 
entitled "Studies in Deep Strata Radio Communica- 
tions" by J. T. deBettencourt, D. A. Hedlund, 
R. A. Sutcliffe, L. Ames, J. F. Frazier, and A. Orange, 
issued by Raytheon Company, Norwood, Mass. 
[Oct, 1962]. 
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3.3. Inhomogeneous Media and Guided Waves 
R. E. Collin " and A. A. Oliner 2 



This section is concerned with progress during the 
past three years in those aspects of guided waves 
which are relevant to radiation phenomena. The 
topics chosen for review are waves along interfaces, 
including surface (trapped) waves and leaky waves, 
the excitation of such waves due to sources and to 
geometric discontinuities, propagation phenomena in 
unbounded periodic media, and guiding and radiation 
effects on periodically-modulated surfaces. 

Waves Along Interfaces 

Electromagnetic waves guided by interfaces may 
possess purely real wavenumbers, in which case they 
are completely bound to the interface and are called 
surface or trapped waves, or their wavenumbers may 
be complex. In the latter case, such waves may be 
associated with anisotropic media within closed wave- 
guides, or may be found under various circumstances 
on open-boundary structures; in open regions, such 
complex waves are also called leaky waves, and may 
be either spectral or nonspectral. A systematic dis- 
cussion of these various wave types in open regions, 
including the different kinds of leaky waves, their 
spectral or nonspectral character, and the treatment 
of surface waves as a special case of complex waves, 
has been presented by Tamir and Oliner [1962a]. A 
corollary paper [Tamir and Oliner, L962b] treats the 
manner in which radiation patterns are influenced by 
the presence of these waves. The discussion below 
considers the surface waves and the complex waves 
separately. 
± 

Surface (Trapped) Waves 

Despite the fact that the basic features of surface 
waves are well known, and that a considerable litera- 
ture already exists, a surprising amount of activity 
still continues on this topic. Most of this activity is 
associated with the investigation of the electro- 
magnetic properties of new structures, including those 
consisting of anisotropic ferrites, dielectrics with trans- 
verse variation, or plasmas with negative dielectric 
constant. 

Two excellent reviews of surface wave structures 
and their properties have appeared during this period. 
We should mention first the notable chapter, written 
by Zucker [1961], which appears in the Antenna 
Engineering Handbook, edited by Jasik. Zucker 
[1961] presents a thorough survey of the properties 
of surface waves, methods of measuring these prop- 
erties, characteristics of the 4 most important struc- 
tures which support these waves, and the use of these 
structures as antennas including design principles for 
single antennas and arrays. The second review, by 
Harvey [1960], includes descriptions of various types 
of surface wave structures and their properties, and 

1 Case Institute of Technology, Cleveland, Ohio. 

2 Polytechnic Institute of Brooklyn, Brooklyn, N.Y. 



also considers multiple media, effects of curvature, 
and questions of launching. Both reviews also 
present extensive bibliographies. A paper with a 
general tone, by Lengyel and Mitzner [1961], exam- 
ines the conditions under which surface waves of the 
type which decrease exponentially on both sides of an 
interface can be guided by a single plane interface, 
but does not treat any specific structures. 

A new phenomenon has also been noted in recent 
years that a backward surface wave can exist on an 
interface in a uniform structure. Furthermore, this 
backward wave is a true mode and not a space 
harmonic as occurs in periodic structures. The 
existence of such a backward wave has been observed 
previously in anisotropic media, such as ferrites or 
plasmas in the presence of magnetic fields, and on 
isotropic plasma columns possessing certain inhomo- 
geneities in their cross section. More recently, 
Trivelpiece, Ignatius, and Holscher [1961] have pro- 
posed the consideration of backward waves on a, 
magnetized ferrite rod for use in a backward wave 
oscillator at millimeter wavelengths, thereby elimi- 
nating the need for a delicate periodic structure at 
these wavelengths. They have computed the prop- 
erties of the wave by using a quasi-static approxima- 
tion, and have experimentally verified the computed 
wave characteristics. Thompson [1961] points out 
the distinction between body and surface backward 
waves, and indicates that while the calculations of 
Trivelpiece et al. [1961] apply to the body modes, 
their experiments relate to the surface waves. A 
reply [1962] agrees with these criticisms but states 
that further measurements will be taken under 
different conditions to finally resolve the question. 

Backward surface waves were also found to exist 
on isotropic uniform plasma slabs [Tamir and Oliner, 
1961 ; Oliner and Tamir, 1962] or plasma-coated 
cylinders [Paik, 1962] in the range of negative 
plasma dielectric constant, e p . It was also noted 
[Tamir and Oliner, 1961; Oliner and Tamir 1962] 
that while forward surface waves on plasmas may be 
present only when e p is negative and |e P Kl, these 
backward surface waves occur essentially in the 
range 0<C|e p |<O- The properties of forward surface 
waves on plasma cylinders were also investigated 
[Rusch, 1962]. 

Hessel, Marcuvitz, and Shmoys [1962] have 
examined the properties of forward surface waves 
on an interface between air and a compressible 
plasma. The compressibility of the plasma permits 
the presence of acoustic modes which are coupled to 
the electromagnetic modes at the interface. The 
resulting surface wave includes both the acoustic and 
electromagnetic contributions, and differs from the 
purely electromagnetic solution only at the higher 
frequencies, where the former wave does not exhibit 
the customary high-frequency cutoff. 

A number of studies have been conducted on 
structures which may be regarded as variants of 
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structures with known properties. The V-line 
[Diament et al., 1961] has been devised, and the 
properties of the line have been explored theoretically 
and experimentally; the V-line is a variant of the 
dielectric image line, but the paper also considers 
modifications of the V-line which could be of prac- 
tical interest. A thorough theoretical and experi- 
mental study [Colin, 19601 has also been made of the 
dielectric-loaded through line, which may be viewed 
as a bisected H -guide. The properties of surface 
waves on symmetrical three-layer dielectric sand- 
wiches have been examined theoretically [Richmond, 
19601 and studies were conducted [Schlesinger et al., 
I960] of higher-order hybrid modes on dielectric 
cylinders. The propagation characteristics of ellip- 
tical dielectric waveguides have also been investi- 
gated [King et al., 1960; Yeh, 1962] theoretically and 
experimentally, with the motivation being the pos- 
sible use of this guiding structure at millimeter wave- 
lengths since it was noted experimentally that the 
attenuation of a circular dielectric rod is reduced 
considerably when the rod is flattened. 

Several other investigations of the properties of 
structures which can support surface waves have 
been conducted with specific applications in mind. 
With the object of obtaining a transmission means 
with low attenuation at millimeter and subinillimeter 
wavelengths, detailed studies [Sobel et al., 1961; 
King et al., 1962] were conducted by Wiltse and 
associates of the properties of various types of 
single-wire lines at these frequencies. Among the 
lines examined were uncoated and dielectric-coated 
metal wires of cylindrical cross section, and un- 
coated wires of elliptical cross section. The prop- 
erties of dielectric cylinders encased in a dielectric 
cladding of different dielectric constant were studied 
by Snitzer [19631 to understand phenomena asso- 
ciated with fiber optics. These properties have been 
examined in considerable detail, including those of 
higher modes, small differences between the dielectric 
constants of the cladding and the core, and mode 
coupling between neighboring fibers. Snitzer [19611 
has also proposed the use of such optical fibers as a 
mode selector in optical masers. 

The application of surface wave structures to 
Luneberg lenses has been proposed by Walter [1960] ; 
he considers as possible structures a dielectric sheet, 
a bed of metal posts, and a plate filled with holes, 
all over a ground plane, derives their indices of re- 
fraction, and applies the results to the design of 
several types of Luneberg lenses. Lichtenberg and 
Woodyard [1962] have shown that a waveguide com- 
posed of an air gap between two semi-infinite plasma 
regions can support a mode with an attenuation con- 
stant significantly lower than that of parallel plate 
waveguide with the usual metallic walls, even when 
plasma collision losses are taken into account. The 
mode is a fast wave w T ith a sinusoidal field distribu- 
tion in the air region, and an exponential decay into 
the plasma regions. Yen [1961] has employed the 
coupling between surface waves as a means of deter- 
mining the radiation pattern from a broadside array 
of endfire antennas. 



Studies have been made of surface waves guided 
by structures with a transverse variation in the 
properties of the structure. Vigants and Schlesinger 
[1962] present an approximate method, employing 
the programming of a computer, for the determina- 
tion of the propagation characteristics of surface 
waves on dielectric cylinders with arbitrary radial 
variation in the dielectric constant. Hirsch and 
Shmoys [1962] examine the modifications in the 
properties of a surface wave guided by an air-plasma 
interface when the single-step interface is replaced 
by a double-step transition. 

The properties of surface waves supported by 
anisotropic sheets consisting of ideally closely spaced 
wires have also been examined. Rumsey [1961] has 
developed a new way of solving Maxwell's equations 
which is ideally suited to such unidirectional aniso- 
tropic sheets; his solution leads to circularly polar- 
ized surface waves along the sheet. His method has 
also been applied by Cheo, Rumsey, and Welch [1961] 
to a sheet with spiral anisotropy as an explanation 
for the behavior of frequency-independent equiangu- 
lar spiral antennas. The solution is highly success- 
ful and yields radiation patterns which agree closely 
with measured results on two- and four-element an- 
tennas. Karal and Karp [1963], and Seshadri [1962], 
consider a field solution which includes a surface 
wave propagating along a plane unidirectionally con- 
ducting screen; Seshadri's phrasing of the problem is 
more general in that for his case the waves travel 
obliquely across the wires, but the approach em- 
ployed in both cases by these authors differs from 
Rumsey's [1961]. 

The surface waves supported by a ferrite slab in a 
rectangular waveguide with a transverse d-c mag- 
netic field were examined by Bresler [19601 in an 
attempt to resolve the so-called " thermodynamic 
paradox." He found that one-way transmission 
was indeed obtained when the ferrite slab is located 
against one wall, but that bidirectional propagation 
always occurs when the slab is removed from the 
wall, no matter how small this distance. The para- 
dox is therefore to be resolved by choosing as the 
correct solution that which is obtained by removing 
the slab a finite distance from the wall and then re- 
ducing this distance to zero as a limit, rather than 
choosing the distance to be zero initially. Ishimaru 
[1963] also treated this question, but disagreed with 
Bresler's [1960] approach; he preferred instead, 
SeidePs concept of "intrinsic loss," and presented 
arguments in its defense. Ishimaru [1962] has also 
demonstrated that a unidirectional surface wave 
can exist along a perfectly conducting plane which 
bounds a semi-infinite plasma with a d-c magnetic 
field. 

The last type of guided wave to be mentioned 
here is not strictly a surface wave, but is the so- 
called wave beam, recently devised for application 
at millimeter and subinillimeter wavelengths. These 
wave beams have been treated theoretically b}^ 
Goubau and Schwering [1961], Beyer and Scheibe 
[1962], and by Goubau [1963]; experimental data 
were presented by Christian and Goubau [1961]. 
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These wave beams are reiterated and guided by 
reconstituting the cross-sectional phase distribution 
at certain regular intervals by means of lenses; the 
beams are shown to satisfy orthogonality relations 
like mores in a waveguide. 

Complex Waves, Including Leaky Waves 

The chapter by Zucker [1961], referred to above 
in connection with surface waves, contains an equally 
fine review of the common structures which support 
leaky waves and the characteristics of these struc- 
tures, the design of leaky wave antennas, and an 
extensive bibliography. A paper by Oliner [1963] 
discusses the distinction between spectral and non- 
spectral leaky waves and reviews the role played 
by leaky waves in a variety of electromagnetic 
phenomena in open regions. Among these phe- 
nomena are radiation from plasma layers, radiation 
from periodically-modulated surfaces, log periodic 
antennas, Cerenkov radiation, and Wood's anomalies. 
The properties of leaky waves guided by slits in 
rectangular waveguide are well known. Goldstone 
and Oliner [19(51] have derived the characteristics 
of such waves on slitted cylinders, for both circular 
electric and circular magnetic excitation, by the use 
of radial transmission line relations and a perturba- 
tion technique. r \ ne effects of coupling on the leaky 
wave properties due to the presence of a neighboring 
slit on a plane and on a cylinder have been investi- 
gated theoretically by Nishida [1960a, 1960b]. 

A surface wave on a plane interface is purely 
bound, and propagates without radiation; on a 
curved surface, however, it is known that leakage is 
associated with the wave. The quantitative influ- 
ence of curvature on the attenuation constant of a 
surface wave was considered by Wait [1960], who 
treated the waves excited by a vertical electric 
dipole on an inductively layered sphere. He dis- 
agrees with a previous paper by Barlow, and corrob- 
orates the original conclusions of Elliott. In reply, 
Barlow [1960] indicates that the low attenuation 
values which follow from Elliott's original analysis 
are not in agreement with some preliminary experi- 
ments, and that the question will be explored further. 
Logan and Yee [1962] have derived a simplified 
mathematical expression which permits one to obtain 
more easily the propagation properties of surfaces 
of the type treated by Wait [I960]. 

Some studies have been conducted which demon- 
strate the presence of leaky waves on plasma slabs. 
Tamir and Oliner [1962] showed that E-type leaky 
waves are very strongly excited on isotropic plasma 
slabs in the range for which the plasma dielectric 
constant is positive. A study by Ishimaru [1962], 
referred to above, showed that a unidirectional sur- 
face wave can be supported in a semi-infinite plasma 
by a metal plate in the presence of a d-c magnetic 
field. When the plasma height becomes finite, this 
surface wave is transformed into a unidirectional 
leaky wave. In both of these studies, the leaky 
wave is shown to influence strongly the radiation 
pattern. 



Complex waves of the spectral type also exist; 
they are characterized by the facts that they occur 
in degenerate pairs and do not carry any real power. 
The presence of such modes in anisotropic ferrite 
regions has been recognized previously; a demonstra- 
tion of their existence based on the asymptotic form 
of the characteristic equations has been given by 
Tai [I960]. Tamir and Oliner [1961] have found that 
such complex waves are also present on isotropic 
plasma slabs in the range for which the plasma 
dielectric constant is negative. 

Excitation of Waves on Open Boundary Structures 

In the previous section various structures capable 
of supporting surface and leaky waves were discussed. 
In practice the ultimate usefulness of such structures 
for antenna applications depends on the relative ease 
and efficiency of launching or exciting the desired 
wave types. It is this particular application that 
often motivates the analysis of wave excitation on 
surface wave structures. The problem is essentially 
solved once the fields excited by a unit source of 
electric (or magnetic) current has been found since 
other sources may be considered to be a suitable 
weighted distribution of unit strength sources. 
Interest in problems of this type has continued 
during the past three years despite the fact that the 
basic techniques of analysis were developed earlier, 
together with results for a number of specific struc- 
tures. The following section presents a survey of the 
problems in this category that have been investigated 
in the period since January 1960. 

A second class of problems concerned with surface- 
wave excitation are those 4 in which surface waves and 
a radiation field are excited by a surface wave 
impinging on a discontinuity along the surface wave 
structure, or by a combination of a source and a 
discontinuity. For the most part only those prob- 
lems that can be classified as two-part boundary 
value problems have received extensive treatment, 
the obvious reason being that these can be solved 
by the Wiener-Hopf technique. Studies of specific 
structures lend insight into the mechanism of radia- 
tion from surface-wave end-fire antennas, and hence 
antenna applications are again a strong motivating 
factor for work in this area. Needless to say, these 
two-part boundary value problems with excitation 
by an incident surface wave or line source are much 
more intricate to solve, from a mathematical point of 
view, than those pertaining to excitation of waves 
on a uniform unbounded structure. The section on 
Discontinuity Effects discusses the progress achieved 
in this area over the past three-year period. 

Not all structures that have been analyzed are 
necessarily appropriate structures for antenna appli- 
cations, for example, a plasma slab. Nevertheless, 
it is of considerable interest to examine the types of 
waves that may be excited on these general open 
boundary structures because of the overall insight 
into radiation and diffraction phenomena that is 
gained. Discussion of these is included along with 
those that have potential application in the antenna 
field since there is no reason to differentiate open 
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boundary structures on the basis of possible applica- 
tions. Mathematically they fall into the same 
category. 

Source Excitation 

a. Uniform Planar Structures With Line Sources 

The canonical problem consists of a line source 
located above a reactance plane or a slab of material 
(dielectric, ferrite, plasma, etc.). A solution in the 
form of an integral is found for the field component 
parallel to the line source in terms of transmission in 
the direction perpendicular to the slab. The poles 
of the integrand determine the surface and leaky 
waves that may be excited; the integrand also has 
one or more branch cuts which yield the radiation 
field. The integral is usually evaluated by deforming 
the integration contour into a steepest descent con- 
tour. During this deformation some of the poles may 
be crossed, in which case part of the solution consists 
of discrete modes. The remainder of the solution is 
determined by an asymptotic evaluation of the inte- 
gral along the steepest descent contour to yield 
results valid for the far zone radiation field. 

The above basic procedure was applied by Colin, 
Cassedy, and Kott [1960] to determine the TE waves 
excited on a dielectric slab above a ground plane. 
These authors present a number of computed curves 
for the power radiation patterns and also for the bi- 
directional launching efficiency for the dominant 
TE 2 o surface mode. A range of values for slab thick- 
ness and line source position is considered. Experi- 
mental data were obtained by placing two conducting 
planes perpendicular to the slab and a small distance 
apart; the fields were then not disturbed, since the 
source is an electric current. The agreement with 
theory was excellent and verified that as much as 
97 percent of the total power could be excited as a 
surface wave. 

In a second paper by two of the same authors, 
Cassedy and Colm [1961], an experimental check of 
the importance of the leaky wave as a partial, but 
dominant, representation of the near zone field at 
the guiding surface was made. Using the same 
dielectric-loaded trough waveguide as in the earlier 
paper, it was possible to demonstrate, by suitable 
amplitude and phase measurements of the total field 
near the guiding surface, that the leaky wave can be 
used to respresent a major portion of the total 
radiation field (not including the surface wave) near 
the surface. 

The case of surface wave excitation by a magnetic 
line source in the presence of two parallel, but 
separated, dielectric slabs has been studied in great 
detail by Angulo and Chang [1961]. 

A problem, somewhat related to the above, where 
the dielectric slab is replaced by a plasma slab and 
the source is a magnetic or electric current line 
source located on the ground plane on which the 
plasma slab rests was treated by Tamir [1961] and 
Tamir and Oliner [1961, 1962b]. The plasma is 
characterized as a medium with a dielectric constant 
which may lie between zero and unity or be negative. 



Tamir [1961] has explored this geometry in great 
detail, and has found that for the electric line source 
no surface waves or spectral complex waves are 
present, but only a large number of nonspectral 
leaky waves. For the magnetic line source case, 
many interesting results were obtained, and some 
of these were reported in the previous section. In 
the negative dielectric constant range, surface 
waves [Tamir and Oliner, 1961; Oliner and Tamir, 
1962] of the backward and forward type, and spectral 
complex waves [Tamir and Oliner, 1961], are present. 
When the dielectric constant is positive, only leaky 
waves are present, but it is significant that the near 
field is then dominated by the presence of a single 
leaky wave [Tamir and Oliner, 1962b]. 

The importance of this last result is that the radia- 
tion field may then be determined by a Kirchhoff- 
Huygens integration over the slab surface using the 
leaky wave field as the source distribution. This is 
in actual fact the procedure used in calculating the 
radiation field from a leaky wave antenna. Kadia- 
tion patterns computed in this manner exhibit one or 
more sharp peaks at angles closely related to the 
leaky wave angles. These patterns agree very well 
with those computed directly from saddle point 
evaluations. The analysis by Tamir and Oliner 
[1962b] provides strong analytical justification for 
this approach to evaluating leaky wave antenna 
radiation patterns. 

Along the same lines, an analysis of a specific 
idealized model of a leaky wave antenna (an induc- 
tive grid structure) was given by Collin [1962]. 
It was demonstrated that for this particular struc- 
ture the use of the leaky wave field as the aperture 
distribution yields a radiation pattern in excellent 
agreement with a more rigorously computed pattern 
(the saddle point evaluation). Further considera- 
tions regarding the relation between leaky wave 
poles and radiation patterns have been given by 
Hessel [1962]. 

In the structures discussed above, the impedance 
presented to the incident waves at the plane of the 
slab was isotropic. The surface impedance can 
instead be made anisotropic in various ways. For 
example, this plane may be replaced by the unidirec- 
tionally conducting plane [Karal and Karp, 1963; 
Seshadri, 1962] discussed in the section on Surface 
Waves. We saw there that along such a plane a 
surface wave mode of propagation is possible. 
The launching of this surface wave mode by an 
electric line source has been examined, and it was 
found [Seshadri, 1962] that excitation efficiencies 
exceeding 90 percent were possible when the line 
source is located close to the surface of the screen. 

A situation that could be expected to have features 
in common with the general anisotropic impedance 
surface mentioned above is that of a plasma slab 
with an applied magnetostatic field. In the general 
case, TE and TM modes are coupled together by the 
anisotropy, but if the line source is directed parallel 
to the applied magnetostatic field, the TE and TM 
modes are not coupled together. This problem is 
therefore considerably easier to analyze. Two 
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independent investigations Ishimaru [1963],Seshadri, 
[1962b] of the excitation of waves in a grounded 
plasma halt-space, or on a grounded plasma slab, 
by a magnetic line source oriented parallel with the 
magnetostatic field have been carried out. Seshadri 
[1962b] considered a half-space filled with plasma 
and bounded by a conducting plane. Ishimaru 
[1962] considered the case of several plasma layers 
bounded at the bottom by a conducting plane. An 
interesting feature of the solution (discussed briefly 
under the Surface Waves section) is the existence, 
for certain plasma parameters, of a surface wave 
guided along the perfectly conducting plane. In 
an isotropic medium such a surface wave does not 
exist. The presence of an anisotropic medium makes 
it possible for the magnetic field to satisfy a boundary 
condition of the type 



dx 



v -+\H y =0 



at the conducting plane. This boundary condition is 
of the same type as that encountered for a reactive sur- 
face and it thus predicts the existence of a surface wave. 
A further property of the surface wave is that it is 
unidirectional, i.e., it propagates in one direction 
only at right angles to the source. Seshadri [1962b] 
also examined the nature of (he radiation field and 
the efficiency with which the surface wave could be 
excited. In the slab problem analyzed by Ishimaru 
[1962], a much wider range of plasma parameters 
was considered. In the opaque region under ap- 
propriate conditions it was found that the surface 
wave pole was, in actual fact, a leaky wave pole 
and furthermore that it was located very near to the 
real axis. This resulted in a sharply peaked radiation 
pattern in accord with the usual behavior of leaky 
wave antennas. In the unbounded plasma the 
surface wave is a fast wave, as in a waveguide. The 
presence of an upper boundary perturbs this surface 
wave into a leaky wave that radiates strongly at a 
particular angle to the surface. 

Another modification of the basic problem dis- 
cussed at the beginning of this section occurs in the 
case of a compressible plasma slab. Using a fluid 
model of a plasma, Hessel, Marcuvitz, and Shmoys 
[1962] have demonstrated that the field in a com- 
pressible plasma can be decomposed into two un- 
coupled sets which are called optical and plasma 
modes. In the presence of a finite-thickness slab or 
an air-plasma interface these modes become coupled 
together through the boundary conditions at the 
surface. The formal solution for the fields excited 
by a magnetic line source located above the interface 
of such a compressible plasma medium and air was 
given. The surface wave portion of the solution was 
discussed in the section on Surface Waves. A still 
different modification of the canonical geometry 
occurs in the problem solved by Kornhauser and 
Keller [1963] that of the field due to a point source 
in a stratified medium with a refractive index that 
varies as n=no (l-a 2 2 2 )K, which is a case of duct 
propagation. The held solution is obtained rigorously 



in terms of tabulated functions, and the ray solu- 
tion has also a particularly simple form. 

b. Cylindrical Structures With Line Sources 

Although the majority of work pertaining to sur- 
face wave excitation has been for planar structures, 
some attention has also been given to other geom- 
etries, in particular, cylindrical structures. 

The efficiency of launching the dominant HEn 
(dipole) mode on a circular dielectric rod by a mag- 
netic ring source and that of launching the TM i 
(Sommerfeld-Goubau mode) wave on a finitely con- 
ducting circular rod by a similar ring source was 
investigated by Cohn and King [1962]. For the HE U 
mode the source had a sinusoidal angular amplitude 
variation, while for the TM i mode there was no 
azimuthal variation. Launching efficiencies of 80 
percent or more were demonstrated. 

As was mentioned above in the section on Complex 
Waves, if the surface which supports a surface wave 
is curved in the direction of propagation, the surface 
wave mode is perturbed into a complex mode that 
loses some energy by radiation. For surfaces with 
large radii of curvature the concept of a surface 
wave mode is still a useful one and enables certain 
properties of the curved surface to be described in a 
satisfactory physical way. A corrugated cylinder, 
with waves excited such that they propagate along 
the angular coordinate, is an example of a curved 
surface along which a perturbed surface wave mode 
may be excited. The resonance excitation of a 
corrugated cylinder excited by an axial electric line 
source was studied by Cullen [I960]. The cases of 
an axial magnetic line source and axial slots and 
dipoles were considered by Wait and Conda [1960 
and 1961]. The radiation held was found to have 
an angular dependence factor which contains a 
denominator that can be made almost zero by an 
appropriate choice of surface impedance. For such 
a choice, the power in a given mode may become 
several hundred times greater than that in any of the 
other modes. Consequently, the radiation pattern 
will vary essentially as cos n<t>, if the nth mode is the 
"resonant" one. 

Although simple line sources and dipole sources 
are generally considered when studying surface wave 
excitation, these are not always the most useful 
sources because of limitations such as narrow band- 
width, high Q, etc. Chu and Kilcoyne have, instead, 
used a helix to excite a dielectric rod antenna in an 
efficient manner [Chu and Kilcoyne, 1961]. In a 
short note by Ravid-Weissberg, the use of distributed 
slots along a waveguide was shown to also provide 
efficient surface wave excitation [Ravid-Weissberg, 
1961]. 

c. Surface Waves on Wedges and Cones 

Although planar and cylindrical structures have 
been most frequently the type of structure considered 
in surface wave studies, other geometries such as 
wedges and cones have not by any means been 
bypassed. Simple solutions for these structures, 
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under the assumption of an impedance type boundary 
condition, are possible only if the surface impedance 
varies linearly. If the surface impedance does not 
have this type of variation then, although the wave 
equation is separable, the boundary conditions are 
not possible and mode coupling takes place. The 
excitation of surface waves on wedges and cones 
having such a linearly varying surface impedance has 
been systematically examined in detail by Felsen 
[I960]. Other work pertaining to surface waves on 
wedges in the presence of discontinuities is discussed 
in the next section. 

Discontinuity Effects 
a. Planar Structures 

The most obvious planar discontinuity problem 
that can be posed and which is amenable to solution 
is that classified as a two-part boundary value 
problem. Such problems have a rather long history 
as far as fundamental analytical techniques go. In 
recent years, the well-known Wiener-Hopf method for 
the solution of such problems has been applied to 
obtain solutions to a number of specific surface wave 
structures where the guiding surface is either semi- 
infinite or has different electrical properties in the two 
regions z <^0 and z >0. Excitation is either by an 
incident surface wave, a line source, or waveguide 
aperture. These solutions are of interest in con- 
nection with surface wave antennas which may, at 
least for a first approximation, be treated as a semi- 
infinite surface waveguide with uniform properties. 
The basic study in this area was published a few 
years ago by A. F. Kay; an outstanding feature of his 
solution is that, in spite of the nonelementary func- 
tions required to describe the field structure, the 
expressions for radiated power involve only elemen- 
tary functions and are readily evaluated. 

The case of a reactive half-plane in the region 
z <0 and a perfectly conducting half-plane in the 
z >0 region excited by a magnetic line source at the 
function was solved by Kane [I960]. His analysis is 
very similar to that in Kay's work. The power 
radiated as a surface wave was shown to reach values 
of 80 percent or more as the surface reactance was 
increased. Above 80 percent the rate of increase of 
surface wave power with reactance was small, 
however. 

The radiation field at the reactive surface is zero 
to at least the first order in (^v 9 ) -1 , where r is the 
distance from the source. This principle, referred 
to as the Karp-Karal lemma by Kane, permits the 
interaction of two or more sufficiently spaced dis- 
continuities on the reactive surface to be evaluated 
as in a conventional waveguide with single mode 
propagation. The reason is that along the surface 
the radiation field arising from each discontinuity 
does not interact, to order (kor)~ l , with the ad- 
jacent discontinuities. Utilizing this lemma, Kane 
applied the results of the analysis of the half -plane 
problem to obtain some approximate results for a 
surface wave antenna of finite length [Kane, 1962a, 
1962b, 1963]. The effect of antenna length, re- 
actance, etc., on gain was determined. 



A two-part boundary value problem involving 
unidirectional conducting screens was treated by 
Sheshadri [1962b]. Three cases were considered, 
involving respectively a semi-infinite unidirectional 
conducting screen, a semi-infinite unidirectional 
conducting screen and a perfectly conducting half- 
plane, and finally two different unidirectional con- 
ducting screens. The analysis was again very similar 
to that in the problem considered earlier by Kay. . 

Excitation of surface waves on a reactive plane by 
a wave incident in a parallel-plate transmission 
line was treated by Wenger [1962]. By choosing 
the extension of the reactive surface into the par- 
allel-plate line properly, the coupling to the surface 
wave can be maximized. The author has delayed 
the publication of the formal analysis in a journal 
until numerical results can be computed. The 
solution is similar in several respects to that of Kay 
and others, although in the present problem no 
simple expression for radiated power seems possible. 

b. Wedges 

It was pointed out earlier that the general wedge 
problem with mixed boundary conditions is readily 
solved only if the surface reactance is linearly 
varying. However, for the special case of a right- 
angle wedge it is possible to find a solution when the 
surface impedance is a constant. The trick used 
here is the introduction of an auxiliary function that 
is a linear combination of the wave function being 
sought and its cartesian derivatives. The auxiliary 
function satisfies simple Neuman or Dirichlet 
boundary conditions and, hence, can be determined 
in a straightforward manner. 

By using the above transformation a number of 
authors examined the problem of excitation of surface 
waves on right angled wedges and the radiation 
produced when a surface wave is incident along one 
face of the wedge [Chu et ah, 1962; Karal et al., 
1961]. This work is an extension of earlier work by 
Karal and Karp dealing with surface wave propaga- 
tion along a wedge. 

c. Approximate Methods 

Many of the open boundary structures of interest J 
for surface wave antenna applications are not 
amenable to rigorous solution. For such structures 
approximate methods of analysis must be devised. 
Many of the rigorous solutions for idealized structures 
are of great value for testing the accuracy of approxi- 
mate solutions. This is one reason why it is worth- 
while to obtain solutions for idealized structures 
even though at times it might appear that such 
problems are only mathematical exercises. 

For surface wave problems involving dielectric 
slabs, Jones [1961] has demonstrated that, by re- 
placing the dielectric slab by a generalized mixed 
linear boundary condition involving a linear combina- 
tion of the wave function and its normal and tangen- 
tial derivatives (a technique first introduced by 
Karp and Kane), approximate solutions may be 
obtained that agree quite well with those special 
cases for which rigorous solutions are available. 
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The technique may thus he extended with reasonable 
confidence to slab problems for which exact solutions 
cannot be found, e.g., a semi-infinite slab. 

Another approximate theory has been developed 
by Keller and Kara! [I960]. This is a geometrical 
optics or ray theory and yields excellent results for 
the high-frequency solutions. The theory enables 
the amplitude of the excited surface wave to be 
found. When applied to the problems of a line 
source and wedge with a variable impedance surface, 
and a line source in the presence of a circular cylinder, 
the results were found to agree exactly with the high 
frequency asymptotic values of the exact solutions. 

Oliner and Tamil* [1962b] have used the fact that 
the near field on a plasma slab is dominated by a 
single leaky wave to determine the influence of 
different terminations on the plasma slab. The 
radiation patterns corresponding to two different 
types of discontinuities were evaluated numerically 
for various lengths between the source and the 
discontinuity. It was found that for short lengths, 
about one or two wavelengths, the effect of the 
discontinuity is significant, while for distances as 
large as twenty wavelengths the influence is negli- 
gible. This dependence on length is readily ex- 
plainable from the properties of the leaky wave, 
since for large distances most of the power has 
already been radiated away before the termination 
is reached. It should be noted that this dependence 
on length would have been very different if a surface 
wave, rather than a leaky wave, were involved. 

Felsen has applied the results of his study of 
surface waves on wedges with linearly varying 
surface reactance [Felsen, 1900a] to the study of a 
model for a surface wave antenna [Felsen, 1960b]. 
The structure considered consisted of a constant 
reactance plane extending oxer — °o<^<a, a section 
with reactance decreasing as x~ l over a<ix<ib, and a 
plane with zero reactance from x=b to infinity. 
Estimates of the effects of taper impedance and 
length on gain, beamwidth and sidelobe level, were 
obtained. 

Periodic Structures 

This section is concerned with propagation in 
unbounded periodic media, and in guiding and 
scattering effects associated with periodically modu- 
lated surfaces. The discussion is restricted to those 
phenomena and those specific structures related to 
radiation or to propagation in open regions. It does 
not include, for example, slow wave structures in 
closed waveguides, parametric effects due to time 
modulation, periodicity introduced by nonlinearities 
induced by high-power incident waves, or mutual 
coupling effects in two-dimensional antenna arrays. 
It is assumed that such treatments will appear else- 
where. 

a. Unbounded Structures 

Activity in tin 4 area of periodic structures of the 
unbounded type has not been particularly great in 
recent years. General properties of periodic struc- 
tures within closed boundary waveguides and also 



unbounded structures in free space are well known 
for the most part. On the other hand, only a 
relatively few reasonably accurate solutions for 
specific structures are a Available. Also, a general 
investigation into the anisotropic properties of peri- 
odic structures, and in particular, periodic media 
with other than cubic or tetragonal lattice structures 
has not been adequately carried out to date. Recent 
work on periodic media, directed in part towards 
the latter-mentioned aspects of the problem, is 
briefly reviewed below. 

The well-known periodic medium consisting of 
thin circular disks arranged in a tetragonal lattice 
has been analyzed many years ago by using the dual 
of Bethe's small aperture theory to compute the 
shunt susceptance of a single plane of disks. This 
shunt susceptance may then be used in a conven- 
tional analysis for loaded transmission lines to obtain 
the propagation constant for the unbounded medium. 
Unfortunately, as shown by recent measurements 
[Gardner, 1960], this approach yields accurate 
results only for disks with k a less than about 0.4, 
where k =2ir/\ and a is the disk radius. Eggimann 
[1961] has improved the original calculations by 
extending the Bouwkamp power series expansion to 
obtain higher order expressions for the electric and 
magnetic dipolo moments, as well as higher order 
moments, induced in a disk by an arbitrary incident 
field. Results accurate to order (k a) 3 were obtained 
for the dipole moments. 

In an earlier attempt to improve on the results 
based on Bethe's theory and static interaction fields 
between disks, Collin and Eggimann [1962] evaluated 
the dynamic interaction constants for a planar array 
of disks. However, the use of these dynamic field 
interaction constants, together with zero-order ex- 
pressions for the dipole moments, yields results that 
are even further away from the experimental values. 
This is perhaps not too surprising since the dynamic 
interaction constants are accurate to much higher 
order in k () a than the Be the expressions for the dipole 
moments, and hence the use of other than static 
interaction fields is not justified without first obtain- 
ing more accurate expressions for the dipole moments. 

By using third-order expressions for the dipole 
moments, dynamic interaction constants computed 
for interacting dipole fields (not plane waves), and 
quadripole radiation terms, it is found that computed 
results for the shunt susceptance agrees well with 
the measured results for k^a up to one [Eggimann 
and Collin, 1962]. Beyond this point the analysis 
is too involved to be practical. 

A problem which is the dual of that of a planar 
array of disks is a conducting plane with a doubly 
periodic array of similar apertures. Such a screen 
comprising rectangular apertures has been considered 
by Kieburtz and Ishimaru [1962]. An approximate 
solution based on an integral equation formulation 
was obtained. Numerical results for the aperture 
field and scattering cross section for k d up to 2 were 
computed, wdiere d is the width of the square aper- 
ture. The scattering cross section shows a resonance 
effect for k d~1.5, at which point the transmission 
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through the screen is a maximum. These authors 
have also obtained the equivalent circuit parameters 
and scattering coefficients [Kieburtz and Ishimaru, 
1961] for such a periodically apertured screen by 
means of a variational expression which they re- 
formulated to include higher mode propagation. 

Another classical periodic medium that has been 
reexamined is the two-dimensional strip medium. 
This medium exhibits pronounced anisotropic effects. 
The early analysis of this medium treated each plane 
of strips as a capacitive grating which could be 
considered to load a transmission line periodically 
with shunt capacitive susceptances. As long as there 
is only dominant mode interaction between successive 
planes of strips good results are obtained. The ex- 
tension to higher mode interaction is cumbersome. 
By considering a different modal expansion, Kolettis 
and Collin [1961] were able to apply the results of 
the Carlson-Heins theory for the parallel-plate inter- 
face problem to overcome the problem of close 
spacing betweeu adjacent gratings. A combination 
of the two approaches gave theoretical results in 
excellent agreement with measured values for a wide 
range of spacings. The analysis has also been 
extended to the case of strips arranged in a non- 
orthogonal lattice pattern [Case Institute of Tech- 
nology, 1962]. No extensive numerical results have 
been computed to date, however. 

Kurss [1961] has examined the formal solution for 
eigen waves in a periodic structure composed of 
scatterers in a lattice defined by translational vectors. 
The same author also formulated the transmission 
line approach so as to be applicable to semiconductor 
band calculations [Kurss, 1961]. There are other 
lattice structures besides those defined by a transla- 
tion. In particular, a periodic medium can be built 
rip from planar arrays of scatterers with adjacent 
planes being obtained by a combination of rotation 
and translation. The properties of such structures 
are presently being investigated by H. Haskal at 
Case Institute of Technology. To date it has been 
established that the modes for such media are 
elliptically polarized waves. 

The scattering of waves at the interface between 
free space and an infinite array of thin dielectric 
sheets has been treated by Collin [I960]. The 
dielectric sheets were placed by infinitely thin 
polarization current sheets, and a rigorous solution 
was obtained for the scattering by the use of bilateral 
Laplace transforms. Numerical results for some 
special cases compared favorably with an approxi- 
mate solution obtained previously via the Rayleigh- 
Ritz method. 

b. Periodically Modulated Surfaces 

The studies during this period on the guiding and 
radiation associated with periodically modulated 
surfaces continue to be motivated primarily by 
antenna applications. The problems posed fall into 
two groups. The first, motivated by an earlier 
paper by Thomas and Zucker, considers the relation 
between the radiation pattern and the modulated 
amplitude and phase distribution over an equivalent 



aperture. The second problem, following an earlier 
paper by Oliner and Hessel, involves the relation 
between either the radiation field or the aperture 
amplitude and phase distribution and physical 
structures which are modulated periodically. 

The former of these problems was considered by 
Ishimaru and Bernard [1962], who presented a closed 
form of analysis rather than the series type of solution 
discussed by Thomas and Zucker. The radiation 
pattern is expressed as the product of an array factor 
and the radiation from one complete cycle; for 
sinusoidal modulation, Anger functions are involved. 
The case of constant amplitude, but with sinusoidal 
phase modulation, was discussed in detail and was 
applied to a modulated slow wave structure. 

Bolljahn [1961] presented an exact synthesis pro- 
cedure for the design of a surface which would sup- 
port a prescribed group of surface waves simul- 
taneously. The surface not only requires a periodic 
variation of surface impedance along its length, but 
it would itself have to be modulated. Some general 
properties of the surfaces and the composite waves 
they support were discussed, and it was indicated 
that the method is applicable only to corrugated 
surfaces and not to dielectric slabs because the 
analysis postulates that no energy flow occurs across 
this surface. 

Ishimaru and Tuan [1962] extended the earlier 
analysis of Oliner and Hessel to discuss the frequency 
scanning of antennas in terms of k versus ft diagrams. 
Modifications were introduced to make the results 
applicable to antennas with a finite number of ele- 
ments, and the case of Tchebycheff arrays was 
treated in detail. Further investigations at the 
Polytechnic Institute of Brooklyn on radiation from 
periodically-modulated surfaces were described by 
Oliner [1963]. The "mode-coupling" regions of the 
k versus diagrams were examined in detail, par- 
ticularly in the leaky wave range. It was shown 
that at broadside radiation, the periodic structure 
approach predicts a sharp drop in radiated power and 
severe reflections, in agreement with experimental 
evidence; other interesting effects were found to 
occur at the backward and forward endfire positions. 
The propagation characteristics of waves on sinusoid- 
ally modulated corrugated rods were considered by 
Wang and Kornhauser [1962]. They measured the 
guide wavelengths on two types of modulated rods, 
and performed an analysis based on the Mathieu 
equation. The measured and theoretical values 
agreed rather well with each other, but both measure- 
ment and theory were confined to the range in which 
the guided waves were purely bound. 

The scattering of electromagnetic waves by a 
sinusoidally varying air-dielectric interface in the 
limit of small wavelength has been investigated by 
Jacobson [1962]. He applies the Luneberg- Kline 
analysis to a general air-dielectric interface, and then 
treats in detail the case of a sinusoidal interface. 
The solution is valid only in the small wavelength 
range, and has been used to evaluate quantitatively 
the accuracy of the geometric optics approximation 
for certain parameter values. In the range con- 
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sidered, geometric optics serves as an excellent 
approximation. 

A new theory of Wood's anomalies on periodic 
structures has been developed by Ilessel and Oliner 
[1903 J. Wood's anomalies are resonance effects 
which occur in the scattering of plane waves by 
periodic structures; the effects have been known for 
many years, and the new contribution lies in the 
recognition of their relation to leaky waves. This 
theory based on leaky waves is more complete than 
the previous, multiple-scattering theory, and also 
permits the prediction of resonance effects in other 
scattering phenomena. Wood's anomaly resonances 
were also found by Li and Oliner [1963] on a structure 
composed of a strip grating backed by a metal plate; 
the behavior of these resonances was explained 
readily by the above-mentioned leaky wave theory 
[1963]. The multiple-scattering theory of Wood's 
anomalies has been developed by Twer sky and others 
in a series of earlier papers; this approach was applied 
in detail by Twersky [1962] to the case of an infinite 
grating of circular cylinders. A discussion of other 
features of his paper and other papers concerned 
with multiple-scattering is presented elsewhere in 
this report . 

The properties of periodic structures are being 
found useful in an explanation of the operation of 
log-periodic antennas. Such antennas are members 
of the class of frequency-independent antennas and 
are characterized by a linear variation in the element 
spacing and element length along the structure. If 
this variation is slow, any given small section of the 
structure may be regarded as locally constant. 
Mayes, Deschamps, and Pat ton [1961] were among 
the first to appreciate this viewpoint, and to propose 
the application of periodic structure theory to this 
class of antennas. A recent report [Mayes et al., 
1962] by these authors applies these viewpoints in 
detail to a variety of such antennas, including a 
helix-fed monopole array, the bifilar zigzag antenna, 
and the log-periodic dipole array. Radiation pat- 
terns are presented and some of the features of 
operation are explained in terms of k versus /3 dia- 
grams. Such diagrams have also been employed by 
Oliner [1962, 1963] in an explanation of some of the 
symmetry requirements associated with proper 
operation of these antennas. Mayes and Ingerson 
[1962] have taken extensive near field measurements 
on a periodic analog of the log-periodic dipole array 
in order to achieve a better understanding ot these 
structures. A study of logarithmically periodic 
circuits, including the effect of loss which is intro- 
duced to represent radiation resistance, has been 
conducted by Mittra [1962] and Jones and Mittra 
[1962]. Log-periodic circuits have also been investi- 
gated by Du Hamel [1963], who has considered the 
general properties of such transmission lines. Rumsey 
[1963] has analyzed an idealization of a zigzag 
antenna consisting of an anisotropic sheet which is 
perfectly conducting only along a sinusoidal path. 
He has found the dispersion characteristics of the 
structure and lias demonstrated various similarities 
between the operation of this structure and that of 



the corresponding antenna. Gans [1962] has con- 
tinued this study by considering the behavior of a 
pair of such anisotropic sheets. 
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3.4. Radiation From Plasmas 
L. B. Felsen 1 



This part of the report deals with the radiation 
from applied electromagnetic sources in the presence 
of bounded or unbounded plasmas; the sources may 
be either stationary in space and oscillate in time 
(as in the antenna problem), or they may comprise 
moving charges (as in Cerenkov radiation). A 
distributed source may be regarded as a superposi- 
tion of vector point sources, whence the problem of 
radiation from the latter is fundamental and has 
received some detailed attention during the past 3 
years. Since explicit closed-form solutions for point- 
source radiation problems can be found only for a 
few selected cases, it is generally necessary to syn- 
thesize the point-source field as a superposition of 
wave fields which individually satisfy the source-free 
field equations. A knowledge of the propagation 
characteristics of plane waves (or their equivalent) 
in the medium is therefore essential to the study of 
radiation from localized source distributions. Be- 
cause of the functional complexity of the formal 
solutions resulting after this superposition, detailed 
investigations of radiation problems have been 
carried out so far only for the simplest characteriza- 
tions of the plasma medium. In most of the results 
to be reported below, temperature and ionic effects 
have been neglected and the plasma is assumed to be 
describable solely in terms of an equivalent tensor 
or scalar dielectric constant, depending on whether 
or not an external magnetic field is present. This 
assumption has reasonable validity for certain types 
of plasmas and externally located sources, but is 
bound to be inaccurate for sources in direct contact 
with the plasma since pressure effects, and the 
associated acoustic waves, may then not be ignored. 
These latter phenomena have been included in 
several studies. Because of the complexity of 
wave couplings in inhomogeneous anisotropic regions, 
the analysis of radiation problems has been limited 
primarily to homogeneous media or to inhomo- 
geneous media simulated by piecewise constant 
layers. 

While the representability of a plasma via the 
above-sketched simplified models is open to question, 
the study of electromagnetic radiation phenomena 
for such "artificial" media can shed considerable 
insight on some of the effects to be anticipated under 
more general and more realistic conditions. 

A. Stationary Sources in Isotropic Plasmas 

1. Incompressible Medium 

A number of contributions have appeared on 
the radiation from time-harmonic electromagnetic 
sources in a bounded isotropic plasma characterized 
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where u v is the (angular) plasma frequency, co the 
applied frequency, and v the electron-ion collision 
frequency. Most of the formal results are identical 
with those obtained in previous studies of ordinary 
dielectrics but the structure of the radiated fields 
may differ appreciably from the ordinary case since 
the real part of e is less than unity and may be 
negative. An important difference arises, for ex- 
ample, from the possibility of supporting a true 
undamped surface wave on a plane interface between 
two media having real dielectric constants e and e , 
respectively, provided that e< — e and that the 
wave has an electric field component normal to the 
interface (E mode). As a result, the spectrum of 
waves supported by a plasma layer may differ 
significantly from that of an ordinary dielectric 
layer [Tamir and Oliner, 1961; Oliner and Tamir, 
1962al and the excitation by localized sources of 
interface surface waves or of other wave types 
peculiar to a medium with negative dielectric 
constant is of interest. These aspects have been 
studied by Tamir and Oliner [1962] in connection 
with the problem of radiation from a line source of 
magnetic current covered by a homogeneous plasma 
slab, and they emphasize the role played by the 
leaky waves in determining some of the prominent 
features of the radiation field exterior to the slab 
region. Some related investigations have been made 
by Larson [1962] and Yee [1962]. Yee [1961] has 
also dealt with the coupling between two narrow 
slots covered by a thick plasma layer with 0<e<O, 
and he presents curves for the input and mutual 
admittance of a slot as a function of slot separation. 
Calculations of the fields radiated through a spheri- 
cal plasma shell by a small electric or magnetic 
antenna have been carried out by Raemer [1962], 
who has likewise considered the perturbation in 
antenna impedance caused by the plasma shell, 
and by Marini [1961], and by Katzin and Koo 
[1962]. The excitation of interface surface waves 
by a ring source of magnetic currents concentric 
with a homogeneous cylindrical plasma column 
was treated by Samaddar [1962a], who has also 
considered the surface wave problem for a plasma 
column with continuous radial variation of dielectric 
constant [1962b]. Hasserjian and Higgins [1962] 
have made calculations of the fields excited by a 
slot on a large circular cylinder surrounded by a 
concentric plasma layer. Methods of simulating 
plasmas by artificial media have been discussed by 
Rotman [1962]. A formal initial value solution of 
the time dependent Maxwell field equations has 
been given by Gerwin [1962]. 
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2. Compressible Medium 

In a more sophisticated description of the plasma, 
the assumption of incompressibility is dropped and 
an isotropic pressure gradient of the electron gas is 
taken into account. In this instance, an electro- 
magnetic source will excite both electromagnetic and 
acoustic waves; the relative strength of excitation of 
these waves is a matter of some concern since calcu- 
lations based on the incompressible model are likely 
to be inaccurate under conditions which favor the 
strong excitation of acoustic waves. Comparison 
with experimental measurements has shown that the 
observed impedance of an antenna surrounded by a 
plasma cannot sometimes be satisfactorily explained 
without taking into account the compressibility of 
the medium [Whale, 1962]. The radiation from an 
oscillating electric current element in an infinite, 
homogeneous, collisionless, compressible electron 
plasma has been calculated by Hessel and Shmoys 
[1961]. They obtain a simple closed form solution 
for the electromagnetic and acoustic fields and find 
that the ratio of power in the acoustic wave to that 
in the electromagnetic wave is proportional to 
(c/a) 3 (w /; /a>) 2 , where c and a are the electromagnetic 
and acoustic wave speeds, respectively, in the same 
gas without electronic charge. A similar calculation 
for a quarter wavelength dipole was made by Cohen 
[1962b], who also found that the acoustic contribu- 
tion to the radiation resistance predominates, espe- 
cially when oo~a) p . Since (c/a) is usually a large 
quantity, these results show that compression al 
effects do play a role, although the assumed model 
(which neglects nonlinear effects, collisions, motions 
of heavy particles, Debye shielding of the antenna) 
is likely to be unrealistic very near the source where 
strong field gradients exist. Cohen [1962b] has also 
considered the boundary condition at a perfect con- 
ductor in a compressible plasma and has made cal- 
culations of the positive sheath required to assure 
that no net current flows into the boundary. 

Calculations have also been made of the radiation 
from currents on a sphere embedded in the plasma, 
or from a dipole located at the center of a spherical 
void in the plasma [Hessel and Shmoys, 1961]. If 
the sphere radius is small compared to both the 
electromagnetic and acoustic wavelengths, the power 
ratio is essentially that given above while a large- 
radius sphere reduces the acoustic wave excitation 
substantially; hence, as noted previously, acoustic 
wave effects are expected to be pronounced when the 
plasma extends into the induction field of the source. 
Cohen [1962a] has dealt with the general source 
problem, discussed the boundary conditions at dis- 
continuities in the medium, and formulated equiv- 
alence theorems for the acoustic and electromagnetic 
fields. He has applied these results to the calcula- 
tion of scattering by a tenuous spherical in homo- 
geneity in the plasma. 

The excitation of a compressible plasma half-space 
by a magnetic line current placed in the exterior 
medium has been considered by Hessel, Marcuvitz, 
and Shmoys [1962]. They find that the interface 



can support a surface wave, and that the energy in 
an electromagnetic plane wave incident almost 
normally on the plasma half space can be converted 
entirely into the acoustic wave if co is slightly higher 
than ov The efficiency of conversion is extremely 
sensitive to the angle of incidence and to the wave 
frequency, and the effect is not expected to be signifi- 
cant if the line source is far enough above the 
interface. 

B. Stationary Sources in Anisotropic 
Plasmas 

1. Incompressible Medium 
a. Unbounded Region 

If an incompressible plasma is subjected to an 
external steady magnetic field H oriented along the 
2-axis, the electromagnetic properties of the medium 
are describable by a normalized dielectric tensor 

e, ie, (Pi 
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where e z has the form previously given and ei, e 2 are 
functions of w, co p , v, and w with co r representing the 
electron-cyclotron frequency. Even in the collision- 
less case v=0, with e u e 2 , e z real, the (ordinary and 
extraordinary) plane waves capable of propagating 
in this anisotropic medium may be separated into 
some eight categories distinguished by the values of 
the ratios co/co /; and co/co c ; see Allis [1961], Haskell and 
Holt [1961], Deschamps and Weeks [1962], for recent 
contributions in this area, whence the detailed study 
of the radiation from localized sources embedded in 
such a medium poses considerable difficulties. Two 
principal approaches have been employed to analyze 
the radiation from current elements in an infinite, 
homogeneous medium described by the above di- 
electric tensor: the first is based on the formulation 
of the dipole field in terms of a spectrum of three- 
dimensional plane waves ("cavity" mode approach) 
while the second utilizes a two-dimensional repre- 
sentation in a plane perpendicular to a preferred 
space direction ("guided" mode approach). The 
first procedure, employed first by Bunkin [1957] in 
the U.S.S.R. and subsequently in the United States 
by Kogelnik [1960], Meecham [1961], Ford [1961], 
Kuehl [1961, 1962], Mitra and Deschamps [1962], 
leads naturally to a triple Fourier integral representa- 
tion while the second, used by Abraham [1953], 
Arbel [1960], Arbel and Felsen [1962a], and by Wu 
[1962], involves double Fourier integrals. The cavity 
mode procedure is well suited to the calculation of 
radiated power in an infinite medium [Kogelnik, 
1960] and also permits the representation of the 
dyadic Green's function in an appealing symmetrical 
form; it has recently been employed by Chow [1962] 
to study radiation in media having both a tensor 
permeability and a tensor dielectric constant. The 
guided wave formalism is easily adapted to the 
analysis of stratified anisotropic media comprising 
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piecewise constant layers whose boundary is per- 
pendicular to the symmetry axis [Arbel, 1960; Arbel 
and Felsen, 1962a], and is also more convenient for 
an asymptotic evaluation of the radiation fields since 
it involves double, rather than triple, integrals. 

An important consideration for radiation problems 
in unbounded anisotropic regions is the phrasing of 
a radiation condition. While the requirement of 
bounded fields at infinity assures the proper formu- 
lation when losses are present, it is necessary in a 
lossless medium to replace the usual "outgoing wave" 
condition by an " outgoing energy" condition [Bresler, 
1959]. This is not a trivial distinction, since the 
directions of phase and energy propagation in an 
anisotropic medium generally do not coincide. Arbel 
and Felsen [1962a] have shown how the plane wave 
dispersion (refractive index) curves for the medium 
can be used in imposing the radiation condition, and 
that each of the constituent propagating plane waves 
must carry energy away from a localized source. 

Asymptotic calculations of the radiation fields have 
been made for various possible values of the plasma 
parameters co 6 ./o> and co^/co and for different orienta- 
tions of the dipole source [Arbel, 1960; Arbel and 
Felsen, 1962b; Meecham, 1961; Ford, 1961; Kuehl, 
1961; Wu, 1962]. Because of the anisotropy, the 
detailed field structure depends markedly on the 
plasma parameters and on the angular location of the 
observation point with respect to the magnetic field 
direction. For a restricted range of observation 
angles in certain plasmas, the field may comprise as 
many as four spherical waves which propagate with 
different phase velocities and individually carry 
energy in the direction from the source to the obser- 
vation point; for other conditions, no propagation is 
possible. The number and intensity of these wave 
contributions can be predicted from the refractive 
index curves or their equivalent. In plasmas having 
a refractive index which becomes singular along- 
certain directions, the fields exhibit strong singulari- 
ties and the total radiated power is infinite (yielding 
an infinite radiation resistance for the dipole antenna) 
[Kogelnik, I960]. These nonphysical singularities 
disappear if the discontinuous source function is 
replaced by a distributed source whose strength goes 
to zero gradually, or if loss is included. Enhanced, 
but bounded, field strengths are also encountered when 
the refractive index curves possess inflection points 
[Arbel and Felsen, 1962b]. 

If the applied magnetic field is infinite (co c =oo), 
e 2 = in (2) and the dielectric tensor is diagonal 
(uniaxial medium). For this special case, solutions 
for the point source or line source excitation problem 
can be obtained in closed form [Arbel, 1960; Kuehl, 
1961; Mittra, 1962; Arbel and Felsen, 1962c; Tuan 
and Seshadri, 1962]. Arbel and Felsen [1962b] have 
shown how these closed form solutions can be used 
to study the previously mentioned field singularities 
when o) c is finite. Mittra and Deschainps [1962] have 
expressed the dyadic Green's function for finite w c in 
terms of a modification of the closed form results 
for co c =oo plus a correction term; since the closed 
form expression is dominant in the near field of the 



source, this formulation appears attractive for an- 
tenna impedance calculations. A variational ap- 
proach to the antenna impedance problem has been 
presented by Anient, Katzin, and Katzin [1962]. 

Some general observations about radiation prob- 
lems in transposed anisotropic media have been made 
by Tai [1961] and by Villeneuve [1961]. 

b. Bounded Regions 

Problems of radiation in the presence of confined 
anisotropic plasmas have also received attention. 
Ishimaru [1962] and Seshadri [1962a] have con- 
sidered a magnetic line current in a plasma bounded 
by a perfectly conducting plane, with the external 
magnetic field H directed parallel to the line source. 
They note the asymmetry of the radiation pattern, 
and find that there may exist a unidirectional surface 
wave which is launched only to one side of the source 
and decays away from the bounding plane; the 
asymmetry is reversed upon a reversal of the direc- 
tion of H . Seshadri has calculated the efficiency of 
excitation of the surface wave as a function of the 
source location. Ishimaru has also dealt with a 
magnetic-line-source-excited plasma slab of finite 
width, as have Hodara and Cohn [1962a] and Shore 
and Meltz [1962]; in the latter analysis, the direction 
of H is perpendicular to that of the line current. 
Via the two-dimensional Fourier transform pro- 
cedure, Hodara [1962] has also discussed, and pre- 
sented numerical data for, the radiation from a finite 
slot in a ground plane covered by a lossy anisotropic 
plasma slab but has omitted any consideration of 
surface or interface waves which may be excited. 
Wu [1962] has performed similar calculations for the 
radiation from an electric dipole source located inside 
or outside the slab. Some of the two-dimensional 
surface waves capable of propagating in an aniso- 
tropic plasma slab have been investigated by Hodara 
and Cohn [1962b]. For a slab of negligible thickness, 
Wait [1960a] has studied the possible surface wave 
solutions and has also calculated (1960b) their 
influence on the fields radiated by vertical magnetic 
or electric dipoles. Imposition of the previously 
mentioned radiation condition in the plasma can be 
circumvented for the slab problem, in view of the 
finite width of the anisotropic region; the radiation 
condition in the isotropic half -space is the con- 
ventional one. 

There have been several studies of radiation from 
a plasma half-space. Arbel [1960] has dealt with the 
case where H is perpendicular to the interface, 
and where an electric dipole parallel to H is located 
either inside or outside the anisotropic region. He 
has made an asymptotic calculation of the far field 
and obtained an interpretation in terms of geo- 
metrical rays which are reflected, refracted and 
coupled at the interface, but he did not consider 
possible additional contributions arising from sur- 
face waves or from lateral (or head) waves. The 
latter, intimately connected with the phenomenon 
of critical refraction, are excited by an incident wave 
which is refracted into the exterior region in a direc- 
tion parallel to the interface, and sheds energy back 
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into the first medium by the reverse process. 
Analytically, the lateral waves usually arise from 
branch cut integral contributions which enter into the 
asymptotic field expressions. Their effects were 
studied by Tyras, Ishimaru, and Swarm [1962] for 
the case of a weak external magnetic field H Q directed 
parallel to the interface, with a magnetic dipole 
source parallel to H . They find a type of lateral 
wave which arises from the interface coupling of the 
ordinary and extraordinary rays and seems to have 
no counterpart in isotropic media. Felsen [1962a, b] 
has analyzed radiation from a uniaxially anisotropic 
half-space (oo c =co) y with the magnetic field oriented 
arbitrarily with respect to the interface. For co<o> p , 
he finds that the lateral waves constitute the 
dominant field contribution in those geometrical 
shadow regions of the plasma which cannot be 
reached by the waves radiated directly from the 
source or by waves reflected from the interface. 
The lateral waves in strongly anisotropic plasmas 
may therefore be of greater importance than in iso- 
tropic regions where they generally constitute a 
second-order effect . 

Samaddar [1961] lias obtained modal expansions 
for the fields excited by a ring source concentric 
with a shielded cylindrical column with gvrotropic 
dielectric constant and (or) permeability, and he has 
calculated the dispersion characteristics for the 
plasma case, with emphasis on possible slow surface 
waves. Wait [1961] has given the harmonic series 
solution for the field excited by a, line source in 
the presence of a radially stratified plasma cylinder 
subjected to an axial magnetic field. 

2. Compressible Medium 

The radiation from a magnetic line source in an 
infinite, homogeneous, compressible plasma sub- 
jected to a steady magnetic field //„ parallel to the 
source direction lias been investigated by Seshadri 
[1962b]. He obtains a closed form solution for the 
radiated field in terms of electromagnetic and acoustic 
wave contributions, presents their dispersion char- 
acteristics, calculates the power coupled into each, 
and finds that the power in the electromagnetic wave 
is always smaller than when the source is located in 
free space. 

C. Radiation From Moving Sources 

An electric charge moving along a straight line 
with constant speed v produces electromagnetic 
radiation in a medium which can support waves 
having phase speeds smaller than v (Cerenkov effect). 
In a plasma without an externally impressed mag- 
netic field, the dielectric constant e is less than that 
of vacuum whence Cerenkov radiation cannot occur. 
However, a moving charge may excite acoustic oscil- 
lations provided that v^>a, where a is the acoustic 
wave speed in the electron gas. Calculations to this 
effect were carried out by Abele (1961) and Cohen 
[1961]. Abele calculated the electromagnetic and 
acoustic fields, as well as the spectral distribution of 
radiated energy, for a line charge moving either in 
an infinite compressible plasma, or exterior and par- 



allel to the surface of a plasma half-space. Cohen 
carried out a similar treatment for a point charge in 
an infinite compressible plasma. Tuan and Seshadri 
[1962] have calculated the radiation from a moving 
line charge in a uniaxially anisotropic medium 
(a> c =oo) and have obtained expressions for the spec- 
tral distribution of the radiated energy; see also 
Johnson [1962]. In contrast to the case of an iso- 
tropic medium, a radially expanding spherically 
symmetric shell of charge may radiate in an aniso- 
tropic plasma. For a high-velocity burst, Ford [1961] 
has calculated the radiation fields and spectral dis- 
tribution of energy in the low-frequency range 
co<C<Cw c , co p , and he finds that the radiation is con- 
fined within a narrow cone about the direction of 
the external magnetic field. 

D. Diffraction by Objects in Anisotropic 
Plasmas 

The study of diffraction by objects embedded in 
an incompressible isotropic plasma characterized by 
the aforementioned dielectric constant proceeds as 
to I* ordinary dielectrics and has therefore received 
little attention per se. However, the effects of 
anisotropy cause considerable modification and have 
stimulated some interest. While diffraction prob- 
lems in a tensor dielectric medium are generally 
exceedingly complicated, there exist special cases 
which result in substantial simplification. For the 
uniaxial plasma (H =oo ? e 2 = above), it can be 
shown that the analysis of diffraction by perfectly 
conducting cylindrical scatterers of arbitrary cross 
section can be simply related to the corresponding 
scattering problem in an isotropic medium, provided 
that the cylinder axis is parallel to H [Arbel and 
Felsen, 1962c; Friedman, 1962]. Thus, formal solu- 
tions for such two-dimensional and three-dimensional 
diffraction problems involving a perfectly conducting 
half-plane, wedge, circular cylinder, etc., can be 
found with comparative ease; the real task is the 
meaningful interpretation of these solutions (or of 
their asymptotic approximations) so as to highlight 
in the broadest possible manner the effect of medium 
anisotropy upon the diffraction field. If H in a 
uniaxially anisotropic medium is perpendicular to 
the cylinder axis, two-dimensional diffraction prob- 
lems (independent of the axial coordinate) can like- 
wise be related to equivalent problems in an isotropic 
region [Felsen, 1962c]. 

Simplifications also obtain for certain axially inde- 
pendent problems in a plasma with finite H directed 
parallel to the cylinder axis. Problems of excitation 
by a line source of electric current are completely 
equivalent to the isotropic case; calculation of the 
fields excited by a magnetic line current is more 
involved but can be performed by the separation of 
variables method for the case of the circular cylinder 
[Samaddar, 1962c]. 

E. Summary and Prognosis 

Theoretical activities during the past three years 
have clarified many of the effects to be expected 
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when time-harmonic stationary sources or moving 
charges radiate from or near unbounded or bounded 
isotropic and anisotropic plasmas. To render these 
initial investigations tractable, the plasma has been 
represented by a highly idealized model, and the 
results are of questionable validity when applied to 
certain actual plasma configurations. It may there- 
fore be expected that studies of radiation from 
stationary or moving sources during the next tri- 
ennium will be concerned with more realistic models 
allowing for compressibility and motion of heavy 
particles, and that attention will also be given to 
the calculation of the impedance characteristics of 
variously shaped radiators of finite size in an infinite, 
homogeneous, anisotropic plasma. Because of the 
complicating effects of anisotropy even for the simple 
tensor medium noted above, the influence of inter- 
faces or inhomogeneities on the fields radiated by a 
prescribed source distribution still remains to be 
clarified, especially as regards the excitation of sur- 
face waves or other types of interface waves. Some 
of the many experimental programs may be expected 
to come to fruition during the next three years and 
supply the theoretician with badly needed physical 
data to substantiate his choice of a given plasma 
model. 

New and difficult mathematical boundary value 
problems have also arisen in connection with the 
determination of the radiation or diffraction pattern 
of a source in the presence of variously shaped ob- 
jects embedded in an anisotropic medium character- 
ized by a dielectric tensor as above. To gain an 
understanding of these phenomena for mathemat- 
ically tractable situations, certain idealized shapes 
and medium properties have been considered, 
involving perfectly conducting cylindrical objects 
oriented with their axes parallel or perpendicular 
to the external magnetic field H . Much work still 
remains to be done even in this category, to say 
nothing of the case wherein H is inclined with 
respect to the obstacle axis. It is desirable to rep- 
resent the asymptotic evaluation of the formal 
solution in a manner which highlights and localizes 
the physical scattering contributions from such 
structural effects as surface curvature, edges, etc., 
when the object is large compared to the wave- 
length, or from its volume and orientation when the 
object is small. Much of the effort in diffraction 
theory in anisotropic media is thus likely to be 
devoted to (a) obtaining exact solutions and asymp- 
totic approximations for scattering by simple shapes, 
and (b) developing the extension to anisotropic 
regions of such techniques as Rayleigh scattering 
approximations for long wavelengths or Keller's 
geometrical theory of diffraction for short wave- 
lengths. 
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3.5. Electromagnetic Wave Propagation in Inhomogeneous Plasmas and/or 

Magnetoplasmas 

W.C.Hoffman 1 ' 2 



In addition to the kind of inhomogeneity usually 
associated with electromagnetic wave propagation, 
viz, that of spatially varying permittivity, inhomo- 
geneities of an entirely different nature can arise out 
of the fact that plasmas are characterized by a finite 
temperature distribution. For example, interactions 
resulting from noncentral peaks in the velocity dis- 
tribution function can modulate or grow waves in an 
otherwise uniform plasma. "In coherent scatter" is 
another case in point (of which we shall have more to 
say later). Another such situation is that wherein 
a periodic variation of permittivity, arising from a 
plasma oscillation, causes an incident electromagnetic 
wave to be scattered [Rosen, I960]. Other such 
interactions can occur and generate actual r-f emis- 
sion by the plasma. Since we are restricted to cold 
stationary plasmas, such situations will, however, be 
in the minority. The propagation of electromagnetic 
waves in plasmas whose complex permittivity ex- 
hibits a spatial variation will certainly be the case of 
greatest importance. We first of all consider the 
situation where no ambient magnetostatic field is 
present. 

Many investigations of the transmission character- 
istics of an inhomogeneous plasma have been moti- 
vated by the "plasma sheath" surrounding a reentry 
body in the upper atmosphere. Klein, Greyber, 
King, and Brueckner [1961] have investigated plasma 
sheath interaction with incident microwave radia- 
tion at both normal and oblique incidence over a 
wide range of the parameters of interest (electron 
density, wavelength, and collision frequency). At 
normal incidence absorption dominates reflection. 
For oblique incidence, the absorption was deter- 
mined in the WKB approximation and found to have 
a steady value below 45° but to drop off rapidly for 
larger angles of incidence. The electromagnetics 
of the plasma sheath problem have also been investi- 
gated by Sisco and Fiskin [1959, 1961]. Miller [1962] 
has given a theory, based on a set of pseudo-plane 
wave functions, for electromagnetic wave propaga- 
tion through an exponentially varying two-dimen- 
sional plasma sheath over an infinite plane conductor. 
The case of a plane TE wave in a stratified medium 
whose ionization density varies as exp(z/z ) has been 
considered by Taylor [1961]. Taylor [1962] has also 
investigated the effect of an irregular variation of the 
degree of ionization of a slightly ionized gas, and 
finds, in distinction to Budden, that such irregulari- 
ties give rise to an "anti-Lorentz term" which in 
general would hinder rather than aid whistler propa- 
gation through the lower ionosphere. The reso- 
nance that occurs as a result of the interaction of a 
high frequency electromagnetic signal with electrons 
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trapped in the potential troughs of an ion wave has 
been analyzed by J. E. Drummond [1961b]. Drum- 
mond obtains an approximate expression, based on a 
nonlinear ion wave theory, for the high-frequency 
conduction current density, and is thus able to pre- 
dict the existence of an "ion wave window" for VHF 
electromagnetic waves in a dense plasma (coJ>^>co 2 ), 
the nature and location of the window being con- 
trolled by the frequency and amplitude of an ion 
wave in the plasma. The interaction of low-fre- 
quency electromagnetic waves with a semi-infinite 
plasma has been investigated by Turcotte and 
Schubert [1961] from the standpoint of coupling be- 
tween electromagnetic and acoustic modes. An 
interesting result has been obtained by Ichimaru, 
Pines, and Rostoker [1962] in connection with a 
theoretical analysis of plasma wave instabilities due 
to drift of electrons relative to ions. They find that 
such instabilities are accompanied by a strongly 
enhanced scattering of electromagnetic waves similar 
to critical opalescence in a gas. 

Electromagnetic wave propagation in a magneto- 
plasma with spatially varying conductivity tensor 
unites the analytical difficulties of anisotropy and 
inhomogeneity, and accordingly there has been 
correspondingly relatively less progress in this area. 
A significant advance in the application of the matri- 
zant solution of the first-order vector wave equation 
to propagation through an inhomogeneous magneto- 
plasma has been made by H. B. Keller and J. B. 
Keller [1962]. The propagation of electromagnetic 
waves in an anisotropic stratified medium has also 
been studied by Hougardy and Saxon [1963], who 
showed that the wave equation can be reduced ap- 
proximately to a scalar one-dimensional form, and 
cross coupling taken into account by iteration 
methods. A penetrating study of the absorption 
and reflection of electromagnetic and plasma waves 
through a plasma whose density and/or magnetic 
field are slowly varying functions of position has been 
made by Stix [I960]." Hoffman [1962] has formu- 
lated the case of an axially symmetric magneto- 
plasma as a vector integral equation, which can be 
approximated by a set of scalar integral equations 
having as solutions mean square convergent Neu- 
mann series. The exterior problem for a perfectly 
conducting sphere in an inhomogeneous magneto- 
plasma has also been considered by Hoffman [1963a, 
b] in terms of an asymptotic form of the vector wave 
equation in dipolar coordinates as well as a vector 
integral equation involving the Green's dyadic for a 
sphere in free space. A formal analysis in terms of 
the Green's dyadic for an inhomogeneous aniso- 
tropic medium has also been given by Villeneuve 
[1961]. 

The subject of "in coherent scatter" from ir- 
regularities in electron density has received much 
attention during the last triennium. Since this 
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phenomenon can be regarded as the result of a 
sort of fine-structure in homogeneity in the plasma, 
its theoretical aspects may legitimately be discussed 
here, a complete discussion of 1 he overall phenomenon 
bein^ reserved lor the Commission 4 report. The 
physical model involves radar scattering from local 
random fluctuations in plasma density, the spectral 
width being that corresponding to the ion velocities 
[Bowles, 1961]. American contributions to the 
mathematical theory of such "incoherent scatter" 
have been made bv Bowles, Oclis, and Green [1962], 
Buneman [1961, 1962], W. E. Drummond [1962], 
Renau [I960], Renau, Camnitz, and Flood [1961], 
Renau 1 1962], Rosenbluth and Rostoker [1962], and 
Salpeter [1960a, 1960b, 1961a, 1961b]. The theo- 
retical approach has involved either a generalized 
Nyquist theorem or solution of the Bolt zm aim 
equation, although extensions to nonequilibrium 
plasmas have been carried out by Buneman [1961], 
Rosenbluth and Rostoker [1962], and Renau [1962]. 
Analyses of the case where an ambient magnetic 
field is present have been carried out by Buneman 
[1961 1, extending a previous detailed study by 
Hagiors [1961], and by Salpeter [1961a, 1961b]. 
The effect of an ambient magnetic field is to cause 
discrete spectral peaks to appear at multiples of the 
ion gyrofrequencies whenever the angle between 
the ambient magnetic field and the direction of 
propagation becomes close to 90°. 

Under certain circumstances plasma oscillations 
or other wave motions within a plasma can generate 
electromagnetic radiation. For instance, tin 4 theory 
of excitation of electromagnetic radiation by gradient 
coupling in plasmas, first analyzed by Field [1956], 
has been refilled and extended by Tidman and Boyd 
[1962], Tidman and Weiss [I960', 1961a], and Wetzel 
[1961, 1963]. The amount of electromagnetic radia- 
tion produced by nonlinear coupling between large 
amplitude transverse and longitudinal plasma oscilla- 
tions has been calculated by Tidman and Weiss 
[1961b]. Harris and Simon [I960] have shown that 
incoherent radiation from a plasma cannot be cal- 
culated on the basis of the Vlasov equation alone. 
Kino [1960] has developed a traveling wave para- 
metric amplifier theory employing a plasma (or 
magnetoplasma) as the nonlinear propagating struc- 
ture, and some refinements of the theory for small 
diameter plasma columns have been given by 
Paik [1962]. 

Electromagnetic Wave Propagation in Homogeneous 
Magnetoplasmas 

In this part of the report we consider the inter- 
action of an electromagnetic wave with a uniform 
plasma containing an ambient magnetic field. 
However, phenomena involving energetic particles 
(such as synchrotron radiation) and transport 
processes (such as electron beams and hydromag- 
netic waves) are specifically excluded. 

Generalizations of the classical conductivity tensor 
in the direction of finite temperature effects will be 
discussed first of all. Sen and Wyller [1960] have 
solved the Boltzmann equation for a uniform 



Lorentz gas in an ambient magnetic field by expand- 
ing the velocity distribution function according to 
the ( 1 hapman-Enskog method and thus arrived 
at a generalized Appleton-Hartree formula. The 
effects of velocity dependence upon the collison 
frequency between electrons and neutral air molecules 
are fully incorporated in the Sen and Wyller analysis. 
A similar study for the solar atmosphere (proton- 
electron and electron-electron collisions) has been 
carried out by Wyller [1961]. The wave normal 
surfaces associated with a finite temperature Apple- 
ton-Hartree formula have been considered by Allis 
[1961] and Papa and Allis [1961], who discuss not 
only the two electromagnetic waves but also a third, 
plasma wave whose velocity is of the order of the 
sound velocity (the "magnetoacoustic wave")- The 
conductivity of a partially ionized gas in the presence 
of an ambient magnetic field has also been analyzed 
by Drummond, Gerwin, and Springer [1961] and 
Kelly [I960]. An analysis of the polarization of 
electromagnetic waves propagating in a magneto- 
plasma has been carried out by Haskell and Holt 
[1961]. 

A number of books dealing at least in part with 
the electromagnetic phenomena associated with 
magnetoplasmas have appeared since those ref- 
erenced in the last USA National Committee report 
on plasmas [Smullin, 19(H)]. The American publi- 
cations include the volumes edited bv Clauser 
[1960] and J. E. Drummond [1961] and the excellent 
treatise by Rose and (lark [1961]. 

Electromagnetic wave theory for homogeneous 
magnetoplasmas has received considerable attention 
during the last triennium. Some of this work has 
already been described above under generalizations of 
the Appleton-Hartree formula [Papa and Allis, 
1961]. Finite-temperature effects associated with 
longitudinal propagation in a magnetoplasma have 
also been investigated bv Pradhan [1957], Narasinga 
Rao, Verdeyen, and Goldstein [1961], Scarf [1962], 
and Willett [1962]. Interesting effects are found to 
occur in tin 4 vicinity of electron-cyclotron resonance. 
The effects of strong magnetic fields upon plasma 
propagation have been investigated by Poeverlein 
[1961], who found, for one propagation mode, per- 
fect guidance along the lines of force and, for the 
other, isotropic propagation. Harley and Tyras 
[1961] and Rothman and Morita [1961] have also 
studied the effects of an intense magnetic field. 
Chow [1962] has extended Bunkin's analysis of radi- 
ation in a magnetoplasma to a gyro-electro-magnetic 
medium and shown that the Sommerfeld radiation 
condition carries over. Kogelnik [1960] and Mee- 
cham [1961] have given dyadic Green's theorems 
appropriate to radiation in uniform magnetoplasmas. 
Wait [1961] has obtained exact solutions for some 
boundary value problems involving a two-dimen- 
sional magnetoplasma, and determined the reflection 
coefficients of such a stratified plasma in planar and 
cylindrical geometries. 

Buchsbaum, Mower, and Brown [1960] have 
investigated the caAaty modes of oscillation of a 
confined magnetoplasma. Hodara and Cohn [1962] 
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have obtained dispersion equations for propagation of 
electromagnetic waves within an infinitely long 
plasma slab containing an axial magnetic field. 

The interaction between electromagnetic waves 
and magnetoplasmas is usually analyzed by means 
of Maxwell's equations plus the Boltzmann equa- 
tion. Solutions of the Boltzmann equation as such 
are therefore of interest. American contributions in 
the area have been made by Pradhan [1957], who 
solved the initial value problem (E(zfi) given) for 
longitudinal propagation and Scarf [1962], who 
treated the corresponding boundary value problem 
(E(0,t) given). The latter analysis was actually 
directed toward the Landau damping of transverse 
electromagnetic wave packets and admits the possi- 
bility of determining the temperature of energetic 
electrons in the exosphere through the proper use of 
whistler data [Liemohn and Scarf, 1962]. The elec- 
tron distribution function of a weakly ionized mag- 
netoplasma in an arbitrary time-dependent electric 
field has been determined by Zmindzinas and 
Wu [1961]. 

We now turn to a consideration of mechanisms for 
radiofrequency emission from cool stationary mag- 
netoplasmas. (Bremsstrahlung, Cerenkov radiation, 
multiple stream interaction, and the like are thus 
specifically excluded.) Hirshfield and Brown [1961] 
have studied the thermal power that can be emitted 
according to KirchliofFs radiation law from a mag- 
netoplasma of radiation temperature T r , and give 
curves for absorption and emission via cyclotron 
radiation. Wyld [1960] has shown that electro- 
magnetic energy can be radiated from the surface 
bounding a finite magnetoplasma whenever plasma 
oscillations exist in the latter. He calculates this 
thermal radiation in the limit of small ratio of 
gyrofrequency to plasma frequency, and finds it is 
less than the synchrotron radiation by orders of 
magnitude. 

The nonlinear interactions of an electromagnetic 
wave incident normally upon a plasma slab in a 
transverse magnetic field have been studied by 
Whitmer and Barrett [1961, 1962]. The dependent 
variables in the Boltzmann equation (including the 
nonlinear terms) and the (coupled) Maxwell equa- 
tions are expanded in a Fourier series in time, and the 
resulting differential equations for each harmonic of 
the electric field solved for plane-wave propagation. 
In the first paper, the properties of the coupled field 
at the second harmonic frequency are discussed in 
terms of the following parameters: ambient magnetic 
field, electron density, electron-neutral particle col- 
lision frequency, field strength of the incident wave, 
and the thickness of the plasma layer. In the second 
paper the dependence of the third and fourth har- 
monics on the same set of parameters is analyzed, and 
the conversion loss obtained for the hth harmonic. 
A general equation is given which enables one to take 
account of interactions with all higher and lower 
harmonics in computing the electric field at the hth. 
harmonic. Coherent nonlinearities excited in a 
magnetoplasma by intense radiation, and the energy 
conversion from the intense wave to a harmonic 



wave or the parametric increase of another, weak 
wave present have also been discussed by Eastman 
[1961]. J. E. Drummond [1961a, ch. 11] has ana- 
lyzed the generation of harmonics and subharmonics 
at multiples of the cyclotron frequency by intense 
microwave illumination of a magnetoplasma. 
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Part 4. Classical Diffraction and Scattering 



4.1. Diffraction and Scattering 
L. B. Felsen 1 and V. H. Weston 2 



A. Introduction and Summary 

Work on electromagnetic diffraction and scattering 
since the last General Assembly has proceeded essen- 
tially along' the lines discussed in the previous report, 
with the great majority of contributions concerned 
with high-frequency effects. We distinguish again 
between exact solutions for specially shaped objects 
(canonical problems) and approximate methods ap- 
plicable to arbitrary shapes. In the first category 
a variety of new problems has been solved, and more 
accurate asymptotic solutions of known problems 
have been derived. No new general techniques of 
analysis of high-frequency scattering have appeared; 
Keller's geometrical theory of diffraction and im- 
proved versions of Fock's theory together seem 
capable of interpreting quasi-optic diffraction phe- 
nomena, and additional evidence has been presented 
to support their validity. Interest has been stim- 
ulated to achieve a similar general formulation for 
diffraction at low frequencies and several contribu- 
tions have appeared in this area. Little progress 
has been made in bridging the analytical gap be- 
tween the low-frequency and high-frequency asymp- 
totic expansions for general scatterers although one 
or two ideas have been presented on this important 
topic. 

Though not exhaustive, the references cited in 
this review are believed to be quite representative 
of American activity in electromagnetic diffraction 
and scattering, and are taken from recognized 
scientific journals only. An exception has been 
made in including papers appearing (in summary 
form) in the preprints of the Symposium on Electro- 
magnetic Theory and Antennas, held in Copenhagen, 
Denmark, during June 1962. Since the Proceedings 
of the Symposium will appear in book form it did 
not seem out of place to anticipate this publication. 
Multiple scattering and scattering in anisotropic 
media has been specifically excluded; these topics 
are covered in other reports. 



The authors express their appreciation to J. B. 
Keller who collaborated in the early phases of prep- 
aration of this summary but was unable, due to 
other commitments, to participate further. 



B. High-Frequency Diffraction 
1 . Canonical Problems 

a. Constant impedance wedge: A number of the 
previously reported studies involving wedges with 
constant impedance boundary conditions have been 
further extended. Karal and Karp have given 
additional applications of their formulation of dif- 
fraction problems involving right-angle wedges 
which have reactive surfaces on one or two sides, 
with major interest centering around the excitation 
of surface waves. The cases considered include 
plane wave diffraction (Karal and Karp, 1962b) and 
excitation by a line dipole at the vertex of the wedge 
(Karp and Karal, 1960). Theoretical aspects of 
diffraction of a surface wave incident on the apex 
from one of the reactive wedge faces were investi- 
gated by Karal, Karp, Chu, and Kouyoumjian 
[1961], while numerical and experimental data are 
given by Chu, Kouyoumjian, Karal, and Karp 
[1962]. Radiation pattern computations and meas- 
urements in the last paper are compared with 
approximate calculations based on the Kirchhoff 
method and good agreement was noted. The results 
are of interest in connection with the radiation from 
finite surface wave antennas. 

As regards the now familiar problem of diffraction 
by a lossy half -plane or by the related configuration 
of two coplanar half-planes with different surface 
impedance, Bazer and Karp [1962] have published 
results pertaining to coastal refraction which were 
obtained by them several years ago, while Mar- 
cinkowski [1961] has presented numerical data for 
the scattering from a strongly absorbing surface 
which he has compared with, and found to resemble, 
the scattered fields derived for a Sommerfeld-type 
black screen. 

b. Infinite cylinder: Investigations of high-fre- 
quency diffraction by cylinders and spheres have 
aimed primarily at improving the asymptotic 
representations for perfect conductivity; at assessing 
the influence of imperfect conductivity and of surface 
coatings; and at providing precise numerical results 
computed from the exact, but slowly convergent 
harmonic series, especially in the resonance region 
where neither the quasi-optic nor the low-frequency 
approximations are very accurate. As regards the 
perfectly conducting case, Hasserjian and Ishimaru 



i Department of Electrophysics, Polytechnic Institute'of Brooklyn. 
2 Conductron Corporation, Ann Arbor, Mich. 



490 



[1962] have obtained solutions for the surface cur- 
rents excited by an arbitrarily oriented small slot 
and by axial and transverse half wavelength slots. 
The results are expressed in terms of those for an 
infinite plane sheet times a correction factor account- 
ing for curvature; the latter is given in terms of 
optical ray path coordinates, thereby suggesting the 
applicability of the results to other convex shapes. 
Verification by experiment shows that the formulas 
hold for observation points even very close to the 
slots. Logan et al., [1962] have presented more 
accurate expressions than previously available for 
the zeros in the order v of Hankel functions having a 
fixed argument, and for the pattern function asso- 
ciated with the radiated field from a circumferential 
slot on a circular cylinder. The effect of spherically 
symmetric waves incident on infinite cylinders has 
been treated by Brick [1961] who examines the formal 
and asymptotic connection between the currents 
induced by a normally or obliquely incident plane 
wave and by point source excitation, and by Brysk 
[1960] who has investigated the radar cross section. 
Numerical calculations of the equatorial plane radia- 
tion pattern of a circumferential half -wavelength 
slot in the interval l<ka<30 have been presented 
by Knop and Batista [1961], while the radiation 
pattern of a radial dipole in the low-frequency range 
0.05<(a/X)-\0.5 has been computed by Levis 
[I960]. Here, k = 2w/\ denotes the free-space wave 
number, and "a" the cylinder radius. 

As regards high-frequency scattering from dielec- 
tric coated cylinders, Kodis [1961] has developed ex- 
pressions for the diffracted field produced by a line 
source either inside or outside the coating. The field 
in the region containing the source is expressed as a 
sum of two terms, the first corresponding to the result 
that would be obtained if the outer surface were also 
perfectly conducting. For plane wave scattering, 
perturbation calculations are carried out for the three 
limiting cases of thin dielectric, low refractive index, 
and very small surface curvature. The special case 
of backscattering has also been considered [Kodis, 
1962], and the high-frequency asymptotic results are 
given in terms of the directly reflected rays and the 
rays which suffer internal reflections in the dielectric 
coating before leaving the outer surface in the back- 
scattered direction. Since grazing rays are neglected, 
these results are best for coatings of moderate thick- 
ness. High-frequency formulas for possibly lossy 
coatings have also been presented by Helstrom [1962]. 
Numerical calculations from the harmonic (Mie) 
series expressions have been carried out for plane 
wave scattering by two concentric circular cylinders 
composed of different material [Kerker and Mati- 
jevic, 1961b]; for plane wave backscattering when a 
dielectric shell is spaced a resonant distance from the 
surface of a large, perfectly conducting cylinder 
(Plonus, 1960); for the equatorial plane radiation 
pattern due to a circumferential half-wavelength slot 
in a perfectly conducting cylinder surrounded by 
various adjacent lossy layers [Knop, 1961]; and for 
the radiation pattern of an axial magnetic dipole in 
the presence of a C3dinder having shallow axial 



grooves [Wait and Conda, 1961]. The harmonic 
series expansion for the scattering of a normally 
incident plane wave by a magnetized ferrite cylinder 
was obtained by Eggiman [1960a] for the case where 
the incident electrical vector and the d-c magnetic 
field are parallel to the axis. Numerical calculations 
for thin cylinders were given. 

c. Sphere: The geometrical optics cross section of 
a large perfectly conducting sphere was recomputed, 
but incorrectly, by Vea et al. [1960], and their error 
was pointed out by Kennaugh et al. [1960] and by 
hogan [I960]. In addition, hogan [1960] has given 
terms in the geometric optical scattered field expan- 
sion up to the order of (ka)~ 4 , where "a" is the radius 
of the sphere, thus extend previous knowledge in the 
field. Kazarinoff and Senior [1962] discuss the fail- 
ure of simple creeping wave theory as regards the 
representation of the minor component of the surface 
field diffracted into the shadow region. However, 
as pointed out by Logan (to be published in the IRE 
Transactions PGAP), this phenomenon has previ- 
ously been noted and corrected in the Russian liter- 
ature. A technique proposed by T. T. Wu, in which 
it was possible through the use of boundary layer 
theory to obtain higher order terms in the series 
expansion for the field in the shadow region of a 
circular cylinder, was extended by Tang [1961] to 
the sphere problem. He derived the first two terms 
in (his series for observation points in the shadow 
region, including the vicinity surrounding the axial 
caustic, and the boundary layer near the surface of 
the sphere. Weston [1961] deals with the near zone 
backscattered field produced by a plane wave inci- 
dent on a large perfectly conducting sphere; the 
geometrical optics field is computed by using two 
terms of the order of (ka)~ 2 , thereby making this 
result valid for any location of the observation point. 

Nonperfectly conducting and coated spheres have 
also been treated. Hiatt et al. [1960], discuss fine- 
grained random roughness on the surface of a large 
perfectly conducting sphere via an impedance bound- 
ary condition. The cross section is computed and 
compared with measurements, and it is found that 
the result is within 0.1 db of that for a smooth sphere, 
for roughness depths up to 0.01X. Weston and 
Hemenger [1962] have derived the geometrical optics 
and creeping wave expressions for the high-frequency 
backscattered field produced by a sphere coated with 
a thin layer of material of complex permittivity and 
permeability. Plonus [1961] uses the classical series 
expansion to derive a formal solution and calculate 
numerical values for the diffraction of a plane wave 
by a perfectly conducting sphere with a concentric 
shell, see also Garbaez [1962]. 

High-frequency scalar scattering from a penetrable 
sphere has been treated by Rubinow [1961]. By the 
use of the Poisson sum formula, he expresses the 
exact series for the scattering amplitude in integral 
form and decomposes it into two sets of integrals. 
One set is evaluated in terms of residues and leads to 
the diffracted field contribution, while the other set is 
evaluated by the method of stationary phase and 
gives rise to the geometric optical field, which includes 
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the bow term. The resulting expressions fail in the 
forward and backward directions, but have also been 
modified to include these cases [Rubinow, 1962]. 
Numerical tables of scattering coefficients of light 
have been computed for concentric dielectric spheres 
[Kerker et al., 1962]; for single dielectric spheres 
[Meehan and Hugus, 1961, and Kerker and Matijevic, 
1961]; for partially absorbing spheres [Deirmendjian 
et al., 1961]; and for colored spheres [Chromey, 
I960]. Penndorf [1962b] obtains an approximate 
formula for the forward scattering of light by a 
spherical aerosol, for particles whose radius is greater 
than X. In addition, approximations are obtained 
for the scattering and extinction coefficients for small 
aerosols (r<X) [1962c]. He also discusses the for- 
ward, sideward, and backward scattering patterns 
based on data for spheres with refractive index 
^=1.33 [1962a]. 

Nonhomogeneous transparent bodies have been 
treated via the Born, and associated, approximations 
[Schiff, 1962]. In particular, nonhomogeneous 
spheres or cylinders were studied by Albini [1962] and 
Albini and Nagelberg [1962], respectively, and 
general nonspherical particles by Greenberg [I960]. 

d. Noncircular cylinder and spheroid: The principal 
purpose of studying high-frequency diffraction b}^ 
noncircular cylinders and by spheroids is to obtain 
rigorous asymptotic field representations for objects 
with nonconstant surface curvature. Such results 
may then be compared with predictions derived from 
approximate theories applicable to bodies of arbitrary 
shape, thereby furnishing important tests of the 
validity of the general methods. A radial eigen- 
function solution for the line source Green's function 
for the elliptic cylinder was given by Levy [1960] and 
was evaluated by him in the asymptotic limit where 
all radii of curvature are large relative to the incident 
wavelength. The resulting expression confirmed 
predictions made by Keller's geometrical theory of 
diffraction. A similar confirmation was obtained by 
Clemmow and Weston [1961] for the field scattered 
from a slightly noncircular cylinder having smooth 
shallow periodic corrugations around the periphery; 
this nonseparable boundary value problem was 
analyzed by a perturbation method. Measurements 
of currents on perfectly conducting elliptic cylinders 
[Wetzel and Brick, 1960] verify the validity of the 
Fock approximation near the shadow boundary even 
down to kR =2.2 where R is the radius of curvature 
at the tip; the agreement is not good in the deep 
shadow since the Fock approximation contains as a 
parameter the distance along the tangent ray rather 
than the (correct) distance along the surface. 

Levy and Keller [1960] employed a method of 
analysis similar to that mentioned above for the 
elliptic cylinder to show that the high-frequency 
diffracted field for a spheroid agrees with that calcu- 
lated by the geometrical theory of diffraction if all 
radii of curvature are large compared to the wave- 
length. They then show that for high frequencies 
the electromagnetic backscattered fields are de- 
rivable from those for the two scalar cases, a result 
which is expected to hold generally for nose-on 



incidence on perfectly conducting bodies of revolu- 
tion. The exact solution for a thin oblate spheroid 
illuminated by a spherically symmetric wave cen- 
tered on the axis was treated by Goodrich et al., 
[1962]. They derive a high-frequency expression for 
the field on the surface in the parameter range 
&a>> 1, (H? ) 1/2 < <1, where R Q =b 2 /a is the (smaller) 
radius of curvature at the edge, and u a" and u b" are 
the semimajor and semiminor axes, respectively. 
It is shown that the reflection coefficient at the edge 
is proportional to the square root of kR , the same as 
for the thin elliptic cylinder. The asymptotic treat- 
ment also holds for the limiting case of the disk. 

e. Elliptic cone: One of the most interesting canoni- 
cal bodies, the elliptic cone, has been treated by 
Kraus and Levine [1961]. By the method of sepa- 
ration of variables in a "uniformized" version of the 
sphero-conal coordinate system, they have con- 
structed scalar Green's functions which satisfy either 
the Neumann or the Dirichlet boundary conditions 
on the cone surface. The solution involves simply 
periodic and nonperiodic Lame functions as well as 
Bessel functions. Since a plane angular sector is a 
degenerate surface, the problem of diffraction by such 
a sector is included as a special case. These solu- 
tions can serve as a starting point for the determina- 
tion of tip conditions (this result is to be published 
by Levine) and the diffraction coefficients for a 
vertex needed in the geometrical theory of diffraction. 

/. Slit and circular aperture: By an application 
of the integral equation procedure described in the 
last report, Sheshadri and Wu have calculated the 
high-frequency scattering of electromagnetic plane 
waves by a circular aperture [Seshadri and Wu, 
1960b] and by an infinite slit [Seshadri and Wu, 
1960a] in an infinite, plane, perfectly conducting 
screen. The results for the ratio of the transmis- 
sion cross section to the geometric-optical cross sec- 
tion are given to 0[(ka)~ 5/2 ] for the aperture and to 
0[(ka)~ n/2 ] for the slit, with "a" representing the 
aperture radius and slit width, respectively. In the 
slit problem, emphasis is placed on the case of 
grazing incidence. A novel method of solving the 
problem of diffraction by an infinite, perfectly 
conducting strip was presented by Timman and 
Kleinmann at the 13th UKSI General Assembly. 
Via their procedure, the line source Green's function 
for a strip is represented in terms of a complex 
double integral with variable integration limits — 
a complication compensated for in part by a simple 
kernel involving Bessel functions. Because of the 
involved form of the solution it has been difficult to 
verify that all the required boundary conditions are 
satisfied. Kleinmann [1962] has now proceeded to 
investigate the limiting case of plane wave incidence. 

A comparison between measured field distribu- 
tions induced in a slit by an obliquely incident 
plane wave and between exact values calculated 
from the Mathieu function series has been given by 
Hsu [1960], who has noted satisfactory agreement 
for slit widths ranging from 1.27 to 3.5 wavelengths. 
Macrakis [1960] has carried out a comparison be- 
tween measured values of the plane wave back- 
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scattering cross section of a strip and calculated re- 
sults obtained from the Mathieu function series, 
the high-frequency asymptotic solution of Levine 
[1959], and the variational method. A quasi-optic 
approximation for the slit (or strip) problem was 
employed by Plonsey [1961] who represents the cur- 
rents induced by a cylindrical wave incident on a 
strip in terms of their geometric-optical value 
which is corrected for edge effects by the inclusion 
of equivalent edge line currents. Measurements 
in a parallel plane setup (with the electric field 
parallel to the edge) showed reasonable agreement 
with calculated values for a strip width as small as 
1.6 wavelengths. As is to be expected, the effects 
of the edge currents are most pronounced when the 
wave arrives near glancing incidence. 

g. Unidirectional screen: A relatively new class of 
diffraction problems, involving unidirectionally con- 
ducting surfaces, has received detailed attention. 
Such an anisotropic surface constrains electric 
currents to flow along a specified rectilinear direction 
and is an idealized model of an array of closely spaced 
insulated thin wires. The electric field tangential to 
the preferred direction on the screen is required to 
vanish while the perpendicular electric and the 
tangential magnetic field components pass unchanged 
through the surface. The first problem of this kind 
was considered by Toraldo di Francia [1956] who 
evaluated the low-frequency scattering cross section 
of a unidirectionally conducting disk. Karp then 
showed that the problem of diffraction by a unidirec- 
tionally conducting half plane could be solved by an 
application of the Wiener-Hopf integral equation 
technique, and this recognition has stimulated the 
adaptation of integral equation procedures employed 
previously for smooth screens to the study of a variety 
of discontinuity problems in unidirectionally con- 
ducting surfaces. These surfaces have the additional 
interesting property of being able to support both a 
true surface wave which decays exponentially, and a 
TEM-type (multiwire) mode which propagates along 
the wire direction and decays algebraically, with 
distance from the surface. 

The problem of plane wave diffraction by a slit in a 
plane screen having an arbitrary direction of con- 
ductivity with respect to the slit edges was treated by 
Seshadri who solved the integral equation in the 
high-frequency limit [1961]. He also obtained the 
low-frequency approximation for the plane wave 
scattering from a strip [I960]. In the former case, 
the solution is given up to terms of 0[(ka)~ 5/2 ] while 
in the latter, terms up to 0[(ka) 2 ] are included, 
where (2a) is the slit width and k is the free space 
wavenumber. Seshadri has also considered the 
efficiency of excitation of the surface wave by a line 
source of electric current (1962b), and the diffraction 
of an incident surface wave by the edge of a uni- 
directional semi-infinite screen, a planar junction of 
two different unidirectional semi-infinite screens, and 
a planar junction between a unidirectional and a 
perfectly conducting half plane [1962a]. The ex- 
citation of the TEM-type wave by a dipole source 
was treated by Karal and Karp [1962a]. 



h. Variable impedance surfaces: While the high- 
frequency diffraction by such diverse structural 
effects as slow surface curvature and abrupt (edge) 
discontinuities on a perfectly conducting scatterer is 
reasonably well understood at present, new phenom- 
ena may arise when the surface properties exhibit 
variation. For gradual changes (over a length in- 
terval equal to the local wavelength) , the quasi-optic 
concepts of "local" diffraction may still be expected 
to apply; for rapid changes, however, the rate of 
variation of the surface properties in the vicinity of 
the quasi-optic diffraction point may be expected to 
exert an influence on the scattered field. The effects 
of surface impedance changes near an edge discon- 
tinuity were investigated by Felsen [1961] who 
considered diffraction by a wedge with a linearly 
varying surface impedance and exhibited the de- 
pendence of the edge diffracted field on the rate of 
impedance variation. Variations on a curved surface 
were investigated by Felsen and Marcinkowski [1962] 
and Marcinkowski and Felsen [1962a, 1962b] via the 
model of a circular cylinder with a sinusoidally 
changing surface impedance around the periphery. 
For rapid oscillations they find that the reflected 
field in the illuminated region can be interpreted in 
terms of a spectrum of diffracted rays which are 
analogous to those encountered for a planar grating, 
while slowly changing surface conditions permit the 
combination of diffracted rays into a single specularly 
reflected ray having a reflection coefficient deter- 
mined by the surface impedance at the reflection 
point. Since the latter condition forms the basis of 
the geometric-optical approximation, criteria can be 
derived which delimit the validity of the geometric- 
optical procedure for constructing the reflected field 
from a gently curved object having a variable surface 
impedance. 

i. Finite cylinder: The very difficult problem of 
diffraction by a capped cylinder of finite length has 
also received attention. Kuehl [1961b] employs an 
approximate procedure suitable when the cylinder is 
thin and many wavelengths long, with the source 
not located near the ends. He first obtains the 
induced currents on the infinitely long cylinder and 
then calculates the radiation from these (unchanged) 
currents when the cylinder is truncated, thereby 
neglecting reflections from the ends. Asymptotic 
evaluations of the radiation pattern, carried out for 
a radial electric dipole located near the cylinder 
surface, agree reasonably well with experimental 
measurements. Kieburtz [1962] has used an essen- 
tially rigorous integral equation and function-theo- 
retic method wherein multiple reflections between 
the ends are taken into account, with due regard for 
the edge condition at the ends of the cylinder. He 
has utilized his results to obtain the scattering cross 
section of thin cylindrical wires. Cassedy and 
Fainberg [1960a] have calculated the backscattering 
from thin lossy cylinders of finite length by the 
variational method, using zeroth and first-order trial 
functions. They find, and have confirmed by experi- 
ments on copper and platinum wires, that loss in the 
wires reduces the backscattering cross section but 
leaves the resonant length virtually unaffected. 
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2. Asymptotic Theories 

As anticipated, most of the effort on general aspects 
of high-frequency diffraction has centered around the 
improvement and refinement, and around further 
applications, of the asymptotic theories of Fock and 
Keller which were discussed in the previous Assembly 
Report. Keller [1962a] studied the initial A r alue 
problem for the time-independent wave equation in 
a region exterior to a polygonal cylinder and has 
deduced the asymptotic behavior of the time- 
harmonic solution from the nature of the singularities 
in the time-dependent solution. The result obtained 
for the leading term in the time-harmonic asymptotic 
series in inverse powers of k agrees with that pre- 
dicted by Keller's geometrical theory of diffraction, 
thereby providing a justification of the theory when 
applied to determine the initial term in the high- 
frequency asymptotic expansion for scattering by 
polygonal cylindrical structures. The verification 
for the higher order terms remains to be carried 
out. 3 The geometrical theory of diffraction has also 
been applied by Keller [1960, 1961] to the determina- 
tion of the plane wave backscattering from a finite 
cone, and by Keller and Magiros [1961] to the study 
of the diffracted fields arising from a cylindrically 
tipped half -plane. In the former case, where no 
exact solution is available, agreement was confirmed 
by comparison with experiment. In the latter case, 
comparison was made with values computed from the 
exact series solution and good agreement was found 
for &a>2, where a is the cylinder radius. Another 
application by Keller and Karal [1960] has dealt with 
the excitation of surface waves by a line source above 
a reactive surface. In this instance, complex rays 
are required and the geometrical theory so formulated 
has been checked against certain special configura- 
tions for which exact solutions are available. 

Buchal and Keller [1960] have employed boundary 
layer techniques to deal with corrections to the geo- 
metrical optics field and to the field computed from 
Keller's diffracted ray theory in the vicinity of 
caustics, edges, and shadow boundaries where these 
fields possess either singularities or discontinuities. 
The actual fields are finite and continuous but 
undergo rapid transitions in thin layers termed 
boundary layers about these regions. Specifically 
treated are the problems of the field in the vicinity 
of a two-dimensional smooth, convex caustic and 
also the special case where this caustic degenerates 
into a line; the field in the vicinity of the edge of an 
aperture in a plane screen; and the field in the 
vicinity of a shadow boundary. W. P. Brown 
[1962] also discusses the boundary layer approach 
and estimates the boundary layer thickness (the 
distance from the caustic) to be the order of 
(ka/2)~ 2/2> a, where "a" is the radius of curvature of 
the caustic. 

Refinements of the Fock theory, dealing pri- 
marily with the description of the diffraction field 



3 A recent summary of the geometrical theory of diffraction is contained in 
Keller [1962b]. 



near the shadow boundary on a smooth convex 
object, have likewise been carried out. For convex 
cylindrical objects, Logan and Yee [1962a] have 
performed extensive computations of what has com- 
monly been referred to as the Fock function, and 
they have achieved a uniform formulation which is 
valid in the illuminated, shadow, and transition 
regions. Thus, high-frequency diffraction phenom- 
ena arising from smooth convex objects can be 
characterized in terms of a single function, thereby 
providing a significant unification of the theory. 
Logan and Yee have also shown how this scattering 
function transforms into the solution predicted by 
the geometrical theory of diffraction whence their 
results can be employed to correct that theory in 
transition regions. Related phenomena have been 
investigated by Rubinow and Keller [1961] who 
studied the shift of the geometric-optical shadow 
boundary on an opaque circular cylinder when 
(X/a) is small but finite. Both zero and finite values 
of surface impedance were assumed, and the shadow 
boundary shift was related to the correction in the 
geometric-optical scattering cross section. The re- 
sults are formulated in a manner which suggests 
their applicability also to other convex two-dimen- 
sional and three-dimensional shapes. 

While the preceding considerations have dealt 
with surfaces whose shapes are describable by analytic 
functions, additional diffraction effects arise when 
the surface contours are nonanalytic. The influence 
of a discontinuity in surface curvature was investi- 
gated by Weston [1962] who employed an integral 
equation procedure to evaluate the high-frequency 
currents induced by an incident plane wave on a 
perfectly conducting cylindrical surface comprised 
of two different but smoothly joined parabolic 
half cylinders. He finds that the discontinuity in 
curvature launches diffracted (creeping) waves on 
the surface which are smaller by 0(1 /k) than those 
originating from a discontinuity in slope (i.e., an 
edge). It is again suggestive that the excitation 
coefficient for these creeping waves can also be 
applied to other cylindrical shapes having a dis- 
continuity in surface curvature. 

Various other contributions, listed in section 1, 
have also led to a formulation of asymptotic results 
which may be expected to apply to certain classes of 
scatterers of relatively arbitrary shape. 

From an analysis based on the theorem which 
relates the scattering cross section and the forward 
scattered field amplitude, some general remarks 
have been made concerning the effects of lossy 
coatings on large scatterers [Hiatt et al., 1960a]. 
It was concluded, and verified experimentally from 
measurements on spheres, spheroids and finite 
cones, that the coatings have a negligible effect on 
forward scattering but an appreciable effect on 
backscattering, see also Ryerson [1962] and Seigel 
and Hiatt [1961]. 

3. Experimental Aspects 

While experimental phases associated with various 
investigations are included at other appropriate 
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places in this report, we list here several contribu- 
tions which deal with experimental techniques in the 
measurement of diffraction fields. Plonsey [1962] 
has dealt with the measurement of surface current 
by an electric probe and has paid attention to the 
error introduced by the interaction between the 
probe and the obstacle. Strait and Cheng [1962] 
have used an auxiliary scatterer for the accurate 
experimental determination of RF magnetic fields. 
Greenberg et al., [1961] have discussed a microwave 
technique which has been developed for obtaining 
total cross sections and angular distributions arising 
from the scattering of electromagnetic radiation by 
nonspherical particles whose size is of the order of a 
wavelength. They compare experiments and ap- 
proximate theoretical results for prolate spheroids 
and finite cylinders. C. I. Beard et al., [1962] deal 
with the error inherent in the use of "midfield" 
experimental results to obtain the far zone forward 
scattering by almost transparent spheres with radii 
large compared to the wavelength. They show 
theoretically that by not being in the true "far- 
zone/' a major error occurs in tin 1 phase of the 
forward-scattered field. 

C. Low-Frequency Diffraction 

As in the case of high-frequency diffraction, though 
on a considerably smaller scale, studies of scattering 
of low-frequency waves can be divided into those 
dealing with exact solutions for special shapes, and 
with genera] methods applied to arbitrary shapes. 
In the first category, the celebrated problem of 
diffraction by a circular disk (or aperture) has been 
subjected to increasingly detailed and sophisticated 
study and as been shown to be susceptible to solu- 
tion by a variety of methods, thereby being rem- 
iniscent of the simpler two-dimensional half-plane 
problem. In most of the recent formulations, one 
seeks to specify the field on the disk implicitly by a 
Fredholm integration equation of the second kind 
which can then be solved by power series expansion in 
terms of the small parameter k. 

Heins and MacOamy [1960] have calculated the 
first few terms in the low-frequency asymptotic ex- 
pansion for the scalar field on a disk excited by an 
obliquely incident plane wave, while Bazer and 
Hochstadt [1962] have performed a similar cal- 
culation for the aperture and have also given the 
scattered far fields and the transmission cross section. 
Expressions for the higher-order terms in the induced 
electric and magnetic dipole moments, and in the 
electromagnetic scattered far fields of a small disk 
and aperture have been presented by Eggimann 
[1960b, 1961], and the results have been employed 
by Eggimann and Collin [1962] to deal with the 
scattering from a planar array of circular disks. 

The scalar problem of diffraction of a plane wave 
by a prolate spheroid is considered by Senior [1960a]. 
He expands the exact solution in ascending powers 
of ka up to and including terms of 0(ka) 6 , where 
"a" is the interfocal distance, and also considers the 
particular cases of the sphere and disk. The uni- 



directional conducting infinite strip has been treated 
by Seshadri [1960] who formulates the problem in 
terms of an integral equation and obtains expressions 
for the far zone fields and the first two terms in the 
series expansion of the total scattering cross section. 

While it is suggestive that certain novel integral 
equation methods applied to the disk can also be 
employed to solve for the diffraction by rotationally 
symmetric scatterers of arbitrary shape [Heins, 
1962a, b], no explicit results for other objects have 
as yet been obtained. In general terms, it is desired 
to infer the low-frequency scattering properties of an 
object from a knowledge of the zero-frequency 
(static) solution. An approximate procedure has 
been utilized by Siegel [1962] and by Senior, Darling, 
andHiatt [1962], who seek to determine the unknown 
coefficients in the dynamic scattered field expansion 
by matching to the known static solution on the 
smallest ficticious spherical surface which completely 
encloses the arbitrarily shaped scatterer. Siegel 
has applied this method to the scattering from a 
cone-sphere combination, the static solution for 
which was obtained by Darling. Karp [1962b] 
seeks to bridge the gap between low-frequency and 
high-frequency asymptotic expansions by proceeding 
from a formulation involving a Fredholm integral 
equation of the second kind which is, however, not 
solved by iteration (thereby yielding the customary 
low-frequency expansion) but rather by the method 
of Fredholm determinants which converges for all 
real frequencies. 

Some new electrostatic results have also been 
found. Van Bladel and Mei [1962] calculated 
the capacitance of a metallic cube, using a varia- 
tional principle. The electrostatic dipole moment 
of a dielectric cube introduced into a uniform electro- 
static field is evaluated by Edwards and Van Bladel 
[1961] for a range of dielectric constants between 
1 and 1000 by numerically solving an appropriate 
surface integral equation. 

The effect of absorbing coatings on perfectly con- 
ducting objects of various shapes has been investi- 
gated theoretically and experimentally by Hiatt 
el al. [1960], with the conclusion that the Rayleigh 
scattering cross section is not markedly affected by 
the coating. 

D. Resonance Region 

No new techniques are in evidence for the calcula- 
tion of the scattering by objects having character- 
istic lengths of the order of a wavelength. Available 
results are usually obtained by numerical evaluation 
of series solutions for special shapes, and have been 
referred to at other appropriate places in this report. 

E. Pulse Solutions 

Gardner and Keller [1962] have treated a pulsed 
vertical dipole situated in the plane interface between 
two semi-infinite homogeneous dielectric media. 
They first find explicitly the two-dimensional field 
due to a line dipole by using the conical flow trans- 
formation and then solve the three-dimensional prob- 
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lem via a representation of the field of a point source 
as a superposition of line sources. They show that 
the field on the interface agrees with that obtained 
previously by Pekeris and Alterman. Ting [1960] 
generalizes the results of Keller and Blank for the 
diffraction of a plane Heaviside pulse by a wedge to 
include arbitrary two-dimensional pulses. He also 
considers diffraction of an arbitrary pulse by a cone. 
An approximate method has been used to calculate 
the scattering of a periodic train of short pulses by a 
sphere [Kennaugh, 1961] and by a finite cone capped 
with a spherical segment [Kennaugh and Moffat t, 
1962], and certain properties of the continuous wave 
scattering are deduced from the pulse solution 
[Kennaugh, 1962]. The results for the cone sphere 
are compared with those of Keller [1960], see also 
Siegel [1963a, b]. 

F. Some General Results 

Uniqueness and existence theorems involved in the 
diffraction of electromagnetic waves by a metallic 
obstacle have been considered by Wilcox [1962]. 
C. S. Morawetz [1961] has shown that if a disturb- 
ance which is initially confined to a finite region 
propagates in free space according to the wave 
equation and is reflected from a star-shaped body, 
the rate of decay with time t at a fixed point in space 
is at least like t~ 1/2 . Marcuvitz [1962] has discussed 
the use of abstract operator methods in the formula- 
tion and formal solution of electromagnetic diffrac- 
tion problems. 

In a paper by Felsen and Karp [1960], an exact 
equivalence is established between certain ring- 
source excited three-dimensional, and line source 
excited two-dimensional problems, both in the same 
medium. The former class includes perfectly con- 
ducting bodies of revolution superimposed upon the 
edge of a perfectly conducting half plane, excited by 

a ring source with variation exp ( ± i h \ where <p i 

the azimuthal variable. The second class involves 
line source excitation of perfectly conducting cylin- 
drical objects. 

In connection with two-dimensional diffraction 
problems, Karp [1961] has shown that for a two- 
dimensional outgoing wave function regular in a 
region outside the circle r>a and satisfying the time- 
harmonic wave equation, there exists a convergent 
expansion in the domain r^>a. This expansion can 
be written as a linear combination of the Hankel 
functions of order of zero and one, the coefficients of 
which are expanded in inverse powers of r. Addi- 
tional conditions on the analytic properties of the 
coefficients are proved, and a relationship between 
the leading terms of these coefficients and the com- 
plex radiation pattern is shown. Karp [1962a] has 
also considered the inverse two-dimensional scat- 
tering problem and has given some new prop- 
erties of the plane wave scattering amplitude 
function. If the far-field pattern function for a 
scatterer with Dirichlet boundary conditions has the 
form/(0, o )=/(0— O ), where O is the angle of inci- 
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dence of the plane wave and is the observation angle 
at infinity, then it can be shown that the scatterer is 
a circular cylinder. 

Welch [1960] uses Rumsey's reaction concept to 
derive certain variational expressions for electro- 
magnetic waves scattered from perfect conductors. 
This is later extended to imperfectly conducting 
media through the use of advanced as well as re- 
tarded solutions [Welch, 1961]. Senior has reviewed 
the derivation of the Leontovich impedance bound- 
ary condition [1962], 
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4.2. On Scattering of Waves by Many Bodies* 

J. E. Burke and V. Twersky 1 



1. Purpose 

This report surveys some of the recent analytical 
work on multiple scattering of waves by fixed con- 
figurations and by random distributions of many 
scatterers. Only papers concerned with discrete 
scatterer formalisms for the reduced wave equa- 
tions are discussed explicitly, but representative 
papers of other formalisms are cited in order to 
indicate the relevance of other approaches. 2 Since 
the analogous previous report to URSI [Twersky, 
1960a] contains historical and introductory dis- 
cussions to the several categories of problems of 
interest, the present report will be relatively brief. 
Citations to specific papers and reports will be re- 
stricted more or less to those published during the 
last three years. Work available in research report 
form for the last survey will be cited again if since 
published in Journals. Some previous citations 
will also be mentioned to facilitate discussion of 
recent generalizations; however, these will be in- 
dicated by an asterisk and not listed in the present 
references (e.g., 1959* means look in previous survey 
for particulars). 

2. Fixed Configurations of iV-Scatterers 
2.1. General Procedures 

Most generally, the problems of present interest 
correspond to solutions (SF) of reduced wave equa- 
tions consisting of a source term (<p) plus several 
scattered terms subject to appropriate boundary 
conditions (in a general sense) at the surfaces of the 
scatterers, and governed by different wave equa- 
tions outside and inside the scatterers. The methods 
of primary interest to the multiple scattering prob- 
lem are those which are essentially independent of 
the specific scatterers involved; using such methods, 
we need not apply boundary conditions explicitly 
and need not rework parts of the problem corre- 
sponding to isolated scatterer problems. Thus, 
assuming that the scattered waves (u s ) for the com- 
ponent scatterers in isolation are known, one seeks 
to embody them in functional equations for the mul- 
tiple-scattered fields which make the effects of the 
configuration explicit. 

For concreteness, we consider a source <p= e ikT 
exciting scatterers located at b s (s = l,2, . . . N). 
We write the total scattered wave as 21 (Jcv) = HU s {k\v 
— b s |)6 ikb „ where U s , the multiple-scattered field 
of body s, reduces to its known single-scattered 
value u s as the other bodies recede to infinity. We 



*Prepared for the Commission 6, USA National Committee Report to the 
XlVth General Assembly of URSI (International Scientific Radio Union) 1963. 

1 Sylvania Electronic Defense Laboratories, Mountain View, Calif. 

2 The related URSI survey by W. C. Hoffman (on propagation in random 
media) discusses various papers which analyze distributions as perturbed con- 
tinuous media, as well as other papers which use essentially particle scattering 
methods. 



seek to express 21 and U s in terms of u s , or in terms 
of auxiliary known functions. The most useful 
auxiliaries are the "scattering amplitudes" defined 
through the far-field form V~cH(kr)F(o,i); here V 
stands for 21, U s , or u s ; (H{kr) equals e iklxl , the 
asymptotic form of H^ (kr) , or h^ l) (kr) in one, two, 
or three dimensions; and F (with o and i equaling 
unit vectors in the directions of observation and 
incidence) stands for the corresponding scattering 
amplitudes Q, G s , or g s . (We use a scalar formalism 
for simplicity but the extension to vector problems 
is straightforward.) 

We also consider the associated "scattering co- 
efficients" J%(u), A(i) and a(i) which arise in the 
series representations of V and F in terms of special 
functions, e.g., in the Fourier series form g s (9 0l 00 = 
l>Z a ns(Qi)£ indo - As in the previous report, the systems 
of integral equations G(g) or algebraic equations 
A(a) obtained from a multiple-scattering analysis 
will be labeled "self-consistent" equations, and their 
usual iterated forms in powers of g or of a will be 
called "successive-scattering" series. 

The earlier general multiple scattering formalisms 
were restricted to configurations of nonidentical 
scatterers for which the isolated scatterer problems 
could be solved by separations of variables. Thus 
working in spherical coordinates, Kasterin [1897*] 
developed the required addition theorems to treat 
spheres in terms of Legendre-Bessel series, and 
Zaviska [1913*] and Ignatowsky [1914a, b*] did the 
same for circular cylinders in terms of Fourier-Bessel 
series. Using essentially similar procedures, Row 
[1955*] and Twersky [1952-54*, 1962g] considered 
configurations of circular cylinders, and Urick and 
Anient [1949*], Ament [1959], Morse [1956*], Kor- 
ringa [1947], Segall [1957], Huang and Yang [1957], 
Eyges [1957, 1958], Waterman and Truell [1961], and 
Schick [1961], among others, considered configura- 
tions of spheres. The analogous expansion in terms 
of Mathieu functions was used by Saermark [1959a, 
b, c, 1960] who developed some of the required addi- 
tion theorems in order to treat configurations of 
strips. 

Essentially two kinds of analytical representations 
have been obtained from such special function ex- 
pansion procedures: infinite sets of linear algebraic 
equations A(a) for the multiple-scattered coefficients 
A ns in terms of the known single-scattered coefficients 
A mt ; asymptotic forms (closed or series) for large 
separations (say &6>>1) for the multiple-scattered 
amplitude G s in terms of the single-scattered values 
g t . The equations A(a) have led to iterated series 
expansions, to closed forms for scatterers small com- 
pared to wavelength, and they have also been treated 
variationally and numerically. The asymptotic 
forms are related to Zaviska 's [1913*] finding that if 
two circular cylinders of arbitrary size are in each 
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other's far fields, then the problem reduces essentially 
to that of one cylinder excited by two plane waves. 
This is the basis for the asymptotic expansions de- 
rived by Twersky [1952b*, 1953a*] by retaining only 
the largest terms involving kb in each order of suc- 
cessive scattering of A(a), and for the closed forms 
he obtained by summing resulting geometric 
progressions. 

Today, the separations of variables derivations of 
the algebraic equations A (a) are only of academic 
interest. It has been shown (Twersky, 1962a, b) 
that such systems of equations can be obtained for 
nonseparable as well as separable shapes without 
knowing the addition theorems for the special func- 
tions of the system in which it is convenient to treat 
the problem. Similarly, the separations of variables 
derivation of the asymptotic forms in kb is super- 
fluous. Thus Karp [1953*] obtained the leading 
terms for 0(g) in 1/kb for arbitrary scatterers, and 
Karp and Zitron [1961a, b; 1959*] considered sub- 
sequent terms of the expansion; for two arbitrary 
cylinders, they obtain the asymptotic expansion cor- 
rect to order (kb)~' i/2 , and for two arbitrary scatter- 
ers in three dimensions they obtain the convergent 
expansion correct to order (kb)~ 2 . Twersky [1962a] 
developed a procedure for obtaining the complete 
asymptotic expansion routinely for two-dimensional 
problems, and applied it for N scatterers to exhibit 
the terms correct to order (kb)~ 5/2 ; he also derived a 
closed form in terms of a differential operator for two 
scatterers, which gave the complete asymptotic 
series. Similarly, for three dimensions, Twersky 
[1962b], exhibited the corresponding converging 
series for N scatterers to order (kb)~' A and derived the 
closed form in terms of a differential operator for 
two scatterers; expanding the denominator of the 
closed form for this case generates the exact con- 
verging series in inverse powers of kb for the two 
scatterer problem. 

Although the separations of variables procedure is 
redundant for the above purposes, the representation 
of the fields for nonseparable shapes as series of 
cylindrical functions, or spherical functions, etc., 
may be of interest for other purposes. On the one 
hand, such representations enable us to exploit di- 
rectly for more complex shapes the explicit results 
previously obtained by special function methods for 
the "array part" of a problem involving elementary 
shapes; e.g., one of the procedures used by Burke and 
Twersky [I960*] for the grating of elliptic cylinders 
was based on expanding the G and #'s as Fourier 
series (instead of Mathieu function series), and sim- 
ilarly, Twersky [1962g] used Fourier-Bessel series 
for arbitrary cylinders in order to exploit the 
corresponding computations for circular cylinders. 
In addition, the known characteristics of the special 
function series may obviate convergence questions 
raised by more general representations. (Of course, 
as far as the single-scattered amplitudes g and co- 
efficients a involved in the representations of G and 
A are concerned, one uses the special function forms 
where convenient to treat distributions of the shapes 
in question; however, for such special computations 



these series representations for g and a need only be 
introduced after one has treated the array part of the 
problem.) 

Most of the current work on multiple scattering 
problems is based on integral equations, or on the 
volume integral and surface integral representations 
for the scattered fields obtained by using Green's 
theorem. Thus we may use the symbolic form 

V(r) = { ^ (£| r -r'|), V(r')} =cf( ${ «b n V 

-Vd n &{ Q )dS, 

where c^o is the free-space Green's function; if 

V= 21 (or u s ) , then the integral is over any surface 

that isolates the collection of scatterers (or the 
scatterer) from the observation point; if V=U S7 
the surface also isolates scatterer s from the others. 
The corresponding scattering amplitudes equal 

F(o,i) = [e- iko ' T8 \ V(r s ';i) } . Ignatowsky [1914a*] used 
essentially such forms to obtain integral equations 
for the current density on one cylinder of a grating, 
and related forms have been used by Schwarzschild 
[1902*], Lax [1951*, 1952*], Ekstein [1951*, 1953*], 
Karp [1953*], Storer and Sevick [1954*], Millar 
[I960*, 1961a, b], Schick [1961], and Twersky 
[1956*-1962]. 

Aside from their utility for specific problems, the 
integral representations and Green's function meth- 
ods lead directly to various theorems that the 
functions of interest must satisfy. Thus Q fulfills 

the same theorems as g (e.g., Q (o,i)= Q (— * — o) r 

the energy theorem, etc. (see Twersky [1962a, b]) r 
and this leads to relations between the G s in the set. 
Of course, G s does not in general fulfill relations as 
simple as those for g (i.e., the G s relations take 
account of the other sources which excite scatterer 
s); however, relatively simple relations emerge for 
special distributions [twersky, 1957*-1962]. 

To generalize the previous discussion of the 
separable shapes, Karp (1953*) uses essentially 
U 8 ={H Ql &} at the surface of a scatterer, to obtain 
asymptotic representations for large kb. The pri- 
mary feature of Karp's procedure for scatterers in 
each other's far fields is to use the fact that the con- 
tribution to the excitation of scatterer s arising from 
the scattered field of scatterer t is effectively a plane 
wave on arrival (i.e., U t ^G t H(k\b a -b t )<p); then, 
since the response of s to a plane wave <p is known, the 
response of s to all such far-field multiple-scattering 
excitations can be obtained by superposition. The 
self-consistent set of algebraic equations G(g) so 
obtained can be solved in closed form for special 
configurations, or iterated to obtain the correct lead- 
ing terms of series for G in terms of g and inverse 
powers of kb. To obtain higher order terms of the 
series, Karp and Zitron (1961a, b) derive additional 
terms of the expansion of U t in plane waxes and in- 
verse powers of kb at scatterer s and use a successive- 
scattering procedure to treat the two-dimensional 
case to (kb)~' s/2 and the three-dimensional case to 

(kb)-\ 

To generalize the above, Twersky (1962a, b) uses 
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the complete representation of the fields in terms of 
complex integrals of free-space plane waves, i.e., 
V(r;i) = J*e ik9 ' r F(p,i), where p is a complex unit 
vector, and f represents integration over an appro- 
priate complex range; in two dimensions, p r = 

rcos(0 o — t) and c/t/t? with c equal to the 

general path used by Sommerfeld for H (1) ; in three 
dimensions, //12(p)/2ry is based on WeyPs 

representation for h a) . (The complex integral 
form holds at least for r outside the smallest sphere 
enclosing the scatterer.) The total scattered field 
may thus be written 1/ = 2j t G s (p,i) exp ik[(p — i) 
•6 s + p.r] with > fl f i(o,i)= fl f l (o > i)+2 1 y* fl f l (o,p)(7,(p > i) 
e l}c ($—i)'b ts following directly from the superposition 
principle applied to the plane wave expansion of 

SF— U t ; here 2 1 means s^t, and b ts =b t —b s =b ts b ts . 
Earlier papers on bounded infinite distributions 
[Twersky, 1953-1959*] were based essentially on 
special cases of the above, and a summary of problems 
which can be treated by such representations is 
given elsewhere [Twersky, 1962c]. Millar [1961c, 
1962] independently developed the general two- 
dimensional form G{g) and applied it to the semi- 
infinite grating. 

To obtain the expansions of G t in inverse powers of 
kb, one transforms the complex set of integral equa- 
tions directly to G t (o, i)=g t (o, i) -f-^oF^o, 

b ts )G s (b ts , i), where J ts = &( (kb ts )e- ik ' b t*[Dt s +l] with 
D equal to a known series in inverse powers of 
kb st whose coefficients involve derivatives with re- 
spect to angles. For three-dimensional problems 
the series is convergent; this follows from the fact 
that h ^ can be written as h ( o = cK times a poly- 
nomial in inverse powers of kb — -an expansion used 
previously by Sommerfeld to obtain a series in in- 
verse powers of kr for the single scattering problem. 
For two-dimensional problems the above multiple- 
scattering series based on the asymptotic form of 
H { n ] is asymptotic; however, the exact expansion 
for this case can be obtained in terms of 
e-^ h tsH^{kb ts )L ts and e- i] * h t*H { \>(kb t8 )M ta where 
L and M (which are analogous to the above D) 
follow from the exact Lommel representation for 
H ( n ] (x) as H^ times one polynomial in 1/x plus 
iJ ( i } times another such polynomial. The exact 
series for H£ )f has recently been used by Karp 
(1961) to derive the corresponding exact expansion 
for the single scatterer problem in two dimensions. 
The set of complex integral equations G{g) also 
leads directly to various systems of algebraic equa- 
tions for the scattering coefficients A(a) for separable 
and nonseparable problems without requiring the 
addition theorems for the special functions used to 
represent the a's. Thus Twersky [1962a] sub- 
stitutes Fourier series representations for G and g 
to obtain A{a) in circular cylindrical coordinates 
(for noncircular and circular cylindrical scatterers), 



and substitutes angular Mathieu function series for 
G and g to obtain A(a) in elliptic coordinates. 
Similarly, in three dimensions, Twersky (1962b) 
writes G and g as series of spherical harmonics 
and substitutes into the complex integrals G(g) to 
obtain a system of algebraic equations A (a) for not 
necessarily spherical scatterers. 

2.2. One-Dimensional Problems 

The above has stressed plane wave representations 
which are novel for two and three-dimensional prob- 
lems. However, such representations are the usual 
ones for one-dimensional problems corresponding to 
scattering of plane waves by configurations of parallel 
infinite slabs, or for the analogous transmission line 
problems. Discrete and continuous distributions 
are considered by a variety of methods for various 
applications by Kay and Silverman [1958], Red- 
heffer [1961, 1962], Bellman and Kalaba [1956, 1958, 
1960], Brekhovskikh [1960], Bloembergen and Per- 
shan [1962], Collin [1960], Ferraro [1962], Melzak 
[1962], Richmond [1962], Twersky [1962d], Wing 
[1962], Yen [1962], Baumeister [1962], Young [1962], 
Mendlowitz and Simpson [1962], Catalan [1962], 
Heavans [1960], and others. 

2.3. Several Separable Scatterers in Two Dimensicns 

The results for two arbitrary cylinders and large 
separation to order (kb)~ 3/2 obtained by Karp and 
Zitron [1961a] were also specialized by them to circu- 
lar cylinders. The corresponding specialization for 
two elliptic cylinders was given by Zitron [1962]. 
The closed operational form for two cylinders and 
the expansions g(a) (in terms of Fourier series and 
angular Mathieu functions series) given in Twersky 
[1962a] constitute the complete asymptotic expan- 
sions for large separations for pairs of circular cylin- 
ders and pairs of elliptic cylinders; the explicit forms 
of the terms for N scatterers and the associated series 
g(a) for the separable scatterers give their expansions 
to(kb)- 5/2 . 

Millar [I960*] considers the simultaneous integral 
equations for a row of perfectly conducting cylinders, 
and (using a procedure analogous to Bouwkamp's 
for the single strip) reduces them to a system of 
linear algebraic equations for elliptic cylinders with 
significant dimensions small compared to wavelength 
("narrow cylinders"). He gives the multiple scat- 
tering closed form for two elliptic monopoles, and 
considers the corrections to the single scattered 
values for elliptic monopole plus dipole. Twersky 
[1962a] gives the multiple scattering closed forms 
for arbitrary angles of incidence for two different 
narrow circular cylinders and for two different narrow 
elliptic cylinders for arbitrary boundary conditions; 
these results generalize those obtained originally by 
Zaviska [1913*] and Twersky [1952*] by separation 
of variables. The closed forms for two different 
monopoles, two different dipoles, and two different 
pairs of monopole plus dipole are derived, and shown 
to satisfy the general scattering theorems for the 
configuration. Illustrative graphs for the intensities 
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and phases for a pair of monopoles and for a pair of 
dipoles (e.g., for E parallel and perpendicular 
respectively to a pair of narrow dielectric rods) are 
given in Twersky [1962c] as functions of the separa- 
tion (kb). Computations are made with the closed 
forms using both exact and asymptotic forms of the 
functions of kb; the two sets of computations are 
indistinguishable for kb^>5, and are close enough for 
practical purposes even down to kb=2. 

Hansen [1959] uses integral equations and Bouw- 
kamp's procedure to treat H parallel to several 
narrow slits in a perfectly conducting infinite plane. 
Two slits are considered in some detail, and the 
difference between approximate results for the limit- 
ing case when the two slits coalesce and the cor- 
responding exact result for one slit of double width 
is shown to be small. Numerical results are given 
for several cases, e.g., the transmission coefficients 
for 1, 2, 4, and 6 slits are compared with Miles' 
[1949*] results for the infinite grating. The ratio 
of separation of centers to slit width (d) is 2.5, and 
the range Jcd=0 to 1 is considered. 

Saermark [1959a, b, c, 1960] uses separation of 
variables in elliptic coordinates to treat multiple 
scattering by several coplanar perfectly conduct ing 
strips (H parallel, and angle of incidence arbitrary). 
The corresponding transmission coefficient versus 
scatterer separation for two slits in a plane is com- 
puted and compared with the result for one strip 
[1959a]. Similar results for the transmission coeffi- 
cient for more than two strips [1959b] are compared 
with the results for the infinite grating. 

Related work on antenna arrays includes papers 
by Wu [1961], (lien and Kino- [1961], and Ishimaru 
[1962]. 

2.4. Several Separable Scatterers in Three Dimensions 

Kumagai and Angelakos [I960] obtain approxima- 
tions for back scattering of electromagnetic waves by 
several spheres large compared to wavelength by 
taking into account double reflections. Collins 
[1961] uses integral equations to consider diffraction 
by two rigid disks with a common axis of symmetry. 

Karp and Zitron [1961b] give results to (kb)~ 2 for 
two arbitrary scatterers at large separation. Twer- 
sky [1962b] gives a closed operational form for two 
arbitrary scatterers, and the forms of the terms to 
(kb)~ 3 for N scatterers at large separations; these 
plus the spherical harmonic series for g(a) constitute 
the analogous special results for spheres. He also 
gives closed form approximations for the multiple- 
scattered coefficients for arbitrary angles of incidence 
on two different monopoles, two different dipoles, 
and two different scatterers specified by monopole 
plus dipole terms. 

In the wave mechanics literature, pairs of mono- 
poles and dipoles have been considered by Brueckner 
[1953*], Watson [1953*], Droll and Verlet [1955*], 
and others. Schick [1961] compares the separation 
of variables procedure, Green's function procedure 
based on volume integrals, and the more abstract 
^-matrix formalism (Lippman and Schwinger) for 



scattering by two radially symmetric scattering 
regions. The separation of variables procedure has 
been discussed in some detail for wave mechanical 
scattering and eigenvalue problems by Eyges [1957, 
1958], and related problems are considered by Huang 
and Yang [1957], and Rodberg [1957]. 



2.5. Periodic Structures 

Millar [1961a, b] considers the integral equations 
for an infinite grating of perfectly conducting 
cylinders and rewrites the kernel as an asymptotic 
series plus an additional function in order to facilitate 
consideration of the " grating resonances" or "Wood 
anomalies. " The modified equations this leads to 
are applied to obtain explicit expressions for narrow 
elliptic cylinders. Results are given For E and H 
parallel, both for the transmission grating and for 
the analogous reflection grating of protuberances on 
a ground plane. The behavior for values of the 
parameters near the resonances is emphasized, and 
both qualitative and numerical results are given for 
the two polarizations. Approximations for the finite 
grating are also discussed. 

Twersky [1962g; 1958b*] treats the grating of 
circular cylinders for arbitrary angles of incidence 
and arbitrary boundary conditions by substituting 
Fourier series for the scattering amplitudes into the 
functional equation 6(g) for the general grating 
[Twersky, 1956a*, 1957b*]. The resulting algebraic 
equations A(a) for the scattering coefficients (which 
he also derives by extending the separations of vari- 
ables procedure introduced by Ignatowsky [1914*] 
for normal incidence) are used to construct various 
explicit approximations. For arbitrary angles of 
incidence, closed form approximations in terms of 
known results for one 4 conducting or dielectric cyl- 
inder and certain series of elementary functions 
[Twersky, 1961; 1958c*] are given for monopoles, 
dipoles, and monopoles plus dipoles. Multiple 
scattering or coupling effects are most pronounced 
for closely spaced cylinders and for the Wood 
anomalies. It is shown that the low-frequency 
" packing effects" for E perpendicular to the axis 
merely increases the dipole moment of the isolated 
cylinder; in this range, the circular cylinder within 
the grating is equivalent to an isolated elliptic 
cylinder whose size and shape are independent of 
the angle of incidence. For both polarizations, ex- 
plicit low-frequency closed forms are given which 
include all coupling effects up to poles of order 2 5 . 
For the Wood anomalies, for narrow cylinders, the 
results for G in terms of A (a) are applied for an 
evanescent mode "near grazing" and compared with 
the previous results for 6 ((f). Illustrative plots for 
the resonance situation for the reflection grating of 
protuberances on a ground plane are given in 
Twersky [1962c]; for three propagating modes, the 
intensities and phases versus kb (or angle of inci- 
dence) are compared; for five modes, the phase and 
intensity of one mode versus kb are compared for 
four different values of scatterer radius. 
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A selective survey of grating literature is given by 
Larsen [1962a]. Grids parallel to interfaces are con- 
sidered by Wait [1958*, 1962] and Wait's results are 
applied numerically by Larsen [1962b] . Nondiscrete 
scatterer approaches to grating problems based on 
Rayleigh's [1907a*] Fourier series representation for 
the grating's profile, and other procedures are con- 
sidered by Stroke [1960, 1961], LaCasce [1961], 
Bottema [1959], Namioka [1961], Ingard [1960], 
Marsh [1961], and others. 

Kieburtz and Ishimaru [1961] apply a variational 
method to an infinite planar lattice of square aper- 
tures to consider scattering when several propagating 
modes exist; the equivalent susceptance is obtained 
from an eigenvalue of an energy operator. Numeri- 
cal results for scattering coefficients versus wave- 
length are given for primary and higher order modes. 
Kieburtz and Ishimaru [1962] also consider the 
integral equation for the boundary value problem 
of the periodic lattice of apertures. For square 
apertures of side a and lattice constant d, they give 
numerical results for fa:/ =5 and ka=l, 1.5, and 2. 
Additional analytical results are given by Kieburtz, 
Ishimaru, and Held [1961], and by Kieburtz [1961]. 
Collin and Eggimann [1961] consider scattering by a 
two-dimensional lattice of circular disks in a rec- 
tangular array and obtain closed form approxima- 
tions plus correction terms for the range where the 
single disk can be specified by electric and magnetic 
dipoles. A related waveguide problem is considered 
by Kerns [I960]. 

Problems related to three-dimensional lattices of 
discrete scatterers are considered by Kohn and 
Rostoker [1954], Adler [1961], Korringa [1947], and 
Segall [1957]. 

3. Random Distributions of Discrete Scatters 

A general introduction to this topic was given in 
the previous survey, with emphasis on papers by 
Rayleigh [1899*], Reiche [1916*], Foldy [1945*], 
Lax [1951*, 1952*] and Twersky [1950-1960*]. 
Starting with the multiple-scattering equations for 
the field of a specific configuration, one introduces a 
statistical ensemble of configurations and seeks the 
average coherent field <C^>, the average total 
intensity <^|^| 2 >, and the corresponding average 
Poynting flux. Proceeding heuristically, one ap- 
proximates averages of \I> with two variables held 
fixed by the average for one fixed variable, or uses 
some equivalent approximation to obtain a determi- 
nate set of equations for <^>. These have been 
treated by several different procedures leading to 
explicit approximations for <^> (or equivalently 
for the "bulk parameters" of the medium associated 
with coherent propagation) in terms of single- 
scattered amplitudes g. 

Alternate procedures for the coherent effects in 
"random media' ' are considered by Keller [1960], 
Hoffman [1959*, 1960], Bremmer [1958*, 1962], 
Silver [1962], and others to be cited subsequently. 
For present purposes we note that Keller [1960] 
gives a critical review and comparison of heuristic 
and analytical procedures, Silver [1962] emphasizes 



the probabilistic restrictions implicit in current 
work, Bremmer [1962] gives a comprehensive review 
of the perturbed continuum literature, and Hoff- 
man's book [1960] provides samples of the varied 
statistical problems in this subject. 

The analytical expressions for ^l^l^ — |<^>| 2 r 
and for the corresponding incoherent scattering /, 
are much more complex. These have been con- 
sidered by Foldy [1945*] and by Lax [1951*], but 
relatively little has been done with them to get 
explicit results for specific problems. For special 
distribution, Twersky [1957*-1963] has used <» 
plus the energy theorem to obtain / explicity. 
Most writers ignore the coherent field and use 
essentially particle scattering procedures based on 
the work of Hopf [1934*] and Chandrasekar [1950*] 
to determine the incoherent scattering. Recent 
papers include those by Chu, Churchill, and Pang 
[1962], Grosjean [1962], and Sekera [1962]. 

3.1. Planar Distributions 

Previous results for scattering by random planar 
distributions of arbitrary scatterers [Twersky, 
1957a*] are specialized to elliptic cylinders by Burke 
and Twersky [1963]. Substituting the known scat- 
tering amplitude for an isolated perfectly conducting 
elliptic cylinder into the general results, they obtain 
the corresponding coherent field (intensity and phase) 
and differential scattering cross section (such that 
the total average intensity satisfies the required 
energy theorem). Both "random screens" of elliptic 
cylinders, and distributions of perfectly conducting 
elliptic protuberances on a ground plane, are treated 
for E perpendicular and parallel to the plane of 
incidence. For scatterers small and large compared 
to wavelength, the reflection coefficients, etc., for 
the two polarizations are plotted versus angle of 
incidence for various eccentricities ranging from 
perpendicular strips, through semicircles to the 
other limiting case of "flat strips." For low fre- 
quencies, the real and imaginary parts of the isolated 
scattering amplitudes required in the analysis are 
given to the eighth and sixth powers of frequency, 
respectively [Burke and Twersky, I960]. Analogous 
results for spheres and circular cylinders are applied 
numerically by Twersky [1962c] to compare electro- 
magnetic scattering by a "uniformly rough surface" 
of small perfectly conducting hemispheres on a 
ground plane and by the analogous "striated surface" 
of semicylinders. For both polarizations, for both 
distributions, the following functions are plotted for 
the full range of angle of incidence (normal to 
grazing): the coherent power reflection coefficients, 
the total phase of the coherent reflected field, and the 
"forward" and back scattered differential scattering- 
cross sections per unit area. The graphs show 
"pseudo-Brewster" effects, etc. 

Recent work on random screens and rough sur- 
faces based on other approaches include papers by 
Ament [1960], Beard [1961], Beckmann [1961, 1962a, 
b], Bennett and Porteus [1961], Briggs [1961], Clay 
1960], Davies [1954], Garrison, Murphy, and Potter 
[1960], Hagfors [1961], Hiatt, Senior, and Weston 
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[I960], LaCasce [1961], LaCasce, McCombe, and 
Thomas [1961], Longuet-Higgins [1960a, b, c], 
Marsh [1961], Marsh, Schulkin, and Kneale [1961], 
Senior [1960], Wait and Conda [1961], Winter [1962], 
Siegel and Senior [1962], Eby, Williams, Ryan, and 
Tamarkin [1960], and Mercier [1961]. 

3.2. Random Volume Distributions 

O'Konski [1960] considers the bulk electrical 
properties of suspensions of ellipsoids (ionic polar- 
ization of macromolecules in polyelectrolytes). 
Mandel [1960] discusses the dielectric constant and 
dispersion in suspensions or oriented prolate sphe- 
roids. Amount [1959] used separations of variables 
and other procedures for obtaining the average 
coherent field and the corresponding bulk physical 
parameters of media composed of statistical dis- 
tributions of spheres. He considers the field for a 
plane wave incident on a semi-infinite region of 
scatters, as well as the characteristics of a spherical 
region containing small spheres in suspension 
Approximations for spheres small compared to 
wavelength are given for problems sucn as the 
electromagnetic field in an artificial dielectric of 
dipole scatters, ferromagnetic exchange effects at 
optical frequencies, thermoelectric power of sintered 
mixtures frequency dependence of artificial dielec- 
trics, sound waves in suspensions, and optical 
activity. 

Waterman and Truell [1961] begin with Foldv's 
[1945*] formalism for the case of an electromagnetic 
wave normally incident on a semi-infinite half space 
of spheres and then use separations of variables in 
spherical coordinates to obtain an approximation 
for the propagation coefficient K in terms of the 
forward and hack scattered values of the scattering 
amplitude of an isolated sphere. Essentially the 
same form had been obtained previously for the 
acoustic case of spheres by Urick and Anient [1949*], 
and had been generalized to arbitrary angle of 
incidence by Twersky [1953b*, 1958a*]. To general- 
ize the results to arbitrary scatterers, Twersky 
[1958a*, 1962] starts with the ensemble average of 
the Green s function surface integral representation 
for arbitrary angle of incidence on a slab region 
of distribution. Approximating the multiple-scat- 
tered amplitude G in the kernel by a simple set of 
free-space scattering amplitudes leads to explicit 
results for K and for the bulk parameters (e and /*) 
in terms of g(k' ,k) where k = hi, and k f is either k 
or its image in the slab face. (The corresponding 
approximation for the incoherent scattering is 
based on the energy theorem.) It is pointed out 
that such procedures are "not quite self-consistent/' 
since one assumes initially that a fixed scatterer 
within the distribution is excited by a set of free- 
space fields of the form e ikr (i.e., a set of "micro- 
scopic" waves in &-space), and then finds that the 

form e iK ' r (the 
m /i -space J. Several limi- 
tations of the explicit results are discussed and 
standard procedures ("hole corrections," etc.) for 
continuing what is essentially the start of a pertur- 



excitation consists of waves of the 
'macroscopic" field in i^-space) 



bation procedure around the characteristics of free 
space are mentioned; as they stand, all such results 
are restricted to sparse distributions. 

In order to obtain a more rapidly convergent 
representation for denser distributions than one 
based on g(k' , k), and in order to develop a procedure 
that is fully self consistent on the macroscopic level, 
Twersky [1958d*, 1962e] introduces a generalized 
isolated scatterer amplitude #(k, K) corresponding 
to an obstacle excited in i^-space but radiating into 
t-space. Approximating the field at a scatterer 
fixed with the distribution by the sum of the coherent 
field (waves of the form e' K ' r ) plus the field radiated 
by the fixed scatterer (a wave having propagation 
number k), Twersky [1959a*, 1962f] obtains func- 
tional equations for K, e, and /jl in terms of g(k, K). 
The equations are applied to small spheres of arbi- 
trary e' and // and the "Lorentz-Lorenz form" is 
obtained for each bulk parameter e and /*; it is shown 
that e and fi are independent of the direction of 
incidence and of the polarization of the incident 
field (i.e., the coherent field defines a unique Max- 
wellian medium). As a second approximation, the 
"loss terms" corresponding to incoherent scattering 
are included; these appear in the bulk e and m in the 
roles of electric and magnetic "conductivities." 
Another illustration considered is normal incidence 
on a slab of arbitrarily shaped scattereis large 
compared to wavelength and with parameters close 
to those of free space ("large tenuous scattereis"); 
for this case, the bulk parameters e and \x equal the 
bulk index of refraction. Illustrative computations 
for the coherent phase and attenuation coefficients, 
and for the coherent and incoherent intensities 
versus the fractional volume for a "pure gas" are 
given in Twersky [1962a]. 

The above forms for K in terms of g(k/ k) and 
g(k, K) are compared with that of Lax [1951*, 
1952*] for K in terms of g(K, K) [Twersky, 1960c] 
and the differences in the explicit results for large 
tenuous scatterers are noted. All forms agree at 
low concentrations, but they diverge at higher con- 
centrations: each formalism is essentially the start 
of a different perturbation procedure around the 
characteristics of free space. Twersky [1962h] gives 
an alternative, "more macroscopic," derivation for 
the forms e and ^ in terms of g(k, K) ; the forms 
follow on direct comparison of the integrals for the 
coherent transmitted and reflected fields in terms of 
the multiple-scattered amplitude 67, with the usual 
Green's function representation for the fields trans- 
mitted and reflected by a uniform slab. He also 
compares the approximation for the incoherent 
scattering constructed to satisfy energy conserva- 
tion (say /) with that obtained by inserting an 
exponential "shielding factor" in the volume integral 
for / to account for the "attenuation" that the in- 
coherent scattering undergoes in leaving the slab. 
For large tenuous scatterers, it is shown that the 
"shielded form" is merely the first stage of a series 
of successive scatterings of the incoherent intensity, 
and that / is the corresponding closed form. Thus 
constructing / to satisfy the energy principle auto- 
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matically takes into account all such orders of 
incoherent multiple scattering. Generalizations for 
distributions of nonidentical scatterers are derived, 
and it is shown that the results for such "mixtures" 
reduce consistently to those for identical scatterers 
if the volume available to a scatterer is taken as 
the original volume less the space taken up by the 
total number of scatterers. 

In addition to the papers by Keller [I960], Silver 
J1962], Bremmer [1962], and Hoffman [1960] cited 
previously, recent papers using other methods to 
analyze physical problems in this field include those 
of Parrent, Shore, and Skinner [1962], Skinner 
[1961], Hoffman [1960], Butler [1962], Chernov 
[I960], Tatarski [1961], Yamada [1961], Ford and 
Meecham [1960], Bowhill [1961a,b], Budden [1959], 
Bugnolo [1960a,b; 1961a,b], Yeh [1962], Grosjean 
[1962], Chu, Churchill, and Pang [1962], Sekera 
[1962], Mulliken [1962], Beckmann [1962], Hessemer 
[1961], Kottler [1960], Isihara [1961], Kielich [1962], 
Melzak [1962], and Musgrave [1959]. 

3.3. Generalized Distributions 

The previous sections were limited to either 
periodic distributions or to "gas-like" distributions 
completely specified by the one-particle distribution 
function, say p s . For the periodic case, p s is a series 
of delta functions; for the gas, p s equals the number 
of scatterers divided by the available volume. More 
generally, it is of interest to consider distributions 
which give the results for the case as one limit 
of some parameters, and those for the periodic case 
as the other limit. Consistent generalizations of the 
full range from gas to crystal based on a two-particle 
distribution function pj(\r s —r t \), where / is the 
"radial distribution" function, are at present limited 
to the one-dimensional case for which explicit enough 
forms for/ exist; see discussion in Twersky [1960a,b], 
and work on planar distributions of cylinders by 
Zernicke and Prins [1927*] and by Twersky [1953a*, 
1959b*]. Cole (1961) summarizes recent progress 
and reviews the literature of radial distribution 
functions. 

Results characteristic of " quasi-random 7 ' and 
"quasi-periodic" distributions may also be obtained 
for a " two-phase" system in which some scatterers 
are considered to be in a "gas phase" and some in a 
"crystal phase." The scatterers in both phases con- 
tribute explicitly to the coherent field, but only those 
in the gas phase contribute explicitly to the inco- 
herent scattering. Thus Twersky [1962i] considers 
a two-phase system of iV=2V r +^V" x large tenuous 
scatterers such that N T are in a gas phase and N x 
are in a crystal phase, and obtains explicit results 
for the coherent and incoherent intensities, etc., in 
terms of fundamental scattering parameters. In 
order to investigate the behavior of the fields during 
a "compression process" in which the system starts 
as a pure gas and ends up as a pure crystal, he intro- 
duces a simple model for the change in relative popu- 
lations as a function of total fractional volume (the 
fraction of space occupied by the N scatterers, say 
co ) : as co increases from to some maximum value 



(e.g., 0.74 for a face-centered cubic crystal of 
spheres) , the ratio N x /N (the fraction of crystal-type 
scatters) increases linearly to unity and the system 
changes from a gas to a crystal. Illustrative graphs 
of the fields for the two-phase system are compared 
with those for the pure gas for the full range of co . 
He also applies the results for mixtures [Twersky, 
1962h] to include a "multicomponent" gas phase 
consisting of differently oriented asymmetrical scat- 
terers and different "clumps" (some scatterers stay 
single, some form doublets, triplets, etc.); in general 
the scatterers when clumped give more incoherent 
scattering than otherwise. The final expressions 
are simple enough in form to facilitate inverting 
measurements when seeking such "crystalizing func- 
tions" and "clumping functions" experimentally. 

The theoretical results are compared in detail 
with parallel microwave experiments by Beard, 
Kays, and Twersky [1963] on a large scale dynamical 
model of a "compressible gas" of large styrofoam 
spheres. The theoretical model of a "multicom- 
ponent two-phase" system seems adequate to ac- 
count for many of the observed phenomena. (Other 
large scale dynamical models of "compressible 
gasses" developed to further comparisons of theory 
and experiment under controlled conditions are 
described in Twersky [1962c]. 

3.4. Relations Between Scatterer Statistics and Signal Statistics 

Under appropriate measuring conditions, the en- 
semble averages for configurations of discrete 
scatterers are equivalent to time-averaged measure- 
ments on a collection of scatterers in motion (i.e., 
time constants of the motion large compared to 
periodicity of incident wave and small compared to 
time interval for the measurement). The relevance 
of configuration averages to measured time averages 
can be determined in the course of actual experiments 
[Beard and Twersky, 1958*, 1960a*, 1960b*] and we 
may restrict considerations to situations where these 
are equivalent. Thus we equate time averages of the 
instantaneous signal <9 V T (where c7 = |c7| and r is 
real) to configuration averages of the field \p=<p-{-*LU s 
where <p is the incident wave and U s is the multiple- 
scattered wave of scatterer s. The configuration- 
averaged coherent field <^> = V^ ta (where C 
is the "coherent intensity" and a is the "coherent 
phase") are related to time-averaged measurements 
through <7H<c7e ZT >| 2 and tan a=<c7 sin r>/ 
<c7 cos r>. Equivalently, we may write o7V T = 
-yjCe iaJ r3e iy such that <5V 7 >==0. Similarly, the 
configuration averaged total intensity T=<^\\f/\ 2 ^> 
equals the time-averaged function <c7 2 >. The 
corresponding incoherent intensity is 

i=T-c=<\+\>>-\<*>\*=<&>>-\<&e* 

For many problems, we may use the approximation 
J~S S <|?7 S | 2 >, i.e., 1 is approximately the sum of the 
average multiple scattering intensities of the com- 
ponent scatterers. 
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The papers mentioned in the previous sections 
were concerned with deriving expressions for C, a, 
and 1 in terms of fundamental scattering parameters; 
as indicated above, these can be compared with 
corresponding time-averaged measurements. How- 
ever, additional time-averaged functions can be 
isolated in the course of experiments. Thus for the 
dynamical model of a gas [Beard and Twersky, 
1958*], Beard [1963] has measured the variances of 
the phase quadrature components of the instanta- 
neous field <x 2 > = <^ 2 cos 2 7> and <?/> = < J7 2 
sin 2 7 > (such that <£ 2 > + <?/> = <# 2 > =J) and 
their covariance <^xy^> = < Q 2 cos y sin y^>. For a 
Rayleigh distribution, we have <^xy^>=0 and 
<£ 2 > = <?/ 2 >:=<i7 2 >/2. More generally, how- 
ever, these second moments give the standard devia- 
tion and correlation coefficient for a bivariate descrip- 
tion of the ensemble of signals. Because of this, and 
because the fundamental scattering parameters (size, 
shape, dielectric properties, etc.) were expected to 
enter differently in the second moments than they 
do in C, a, and I, it has been of interest to relate 
these moments to fundamental scattering theory. 

The required relation was recently derived by 
Twersky (1962j). He introduced the "asymmetry 
function" P= |P|e"such that 2P^ 2 = <(^-<^» 2 > 
~2<£7s> 2 . In terms of I and P, the second 
moments of the ensemble of signals are given by 
<£ 2 >=^/+|P| cos p, <?/ 2 > = |J— |P| cos p, and 
<£?/> = |P| sin p. Thus measurements of these 
second moments give 1, \P\, and p, which together 
with C and a provide five independent functions for 
isolating scattering parameters from microwave 
measurements on distributions. Explicit approxi- 
mations of P for the gas, two-phase system, etc., of 
large tenuous scatterers are given by Twersky [1962i]. 

In addition, the scattering function representation 
for P provides explicit relations for the leading terms 
of the moment expansions of other measurables. 

Thus for "small randomness"< c7 2 >-<c7> 2 «K 

+ \P\cos (p-2a), <&>>-C~iI-\P\ cos (p-2a), 
<r>-a~(\P\/C) sin (p-2a), <y>-< r >*« 
[}I—\P\ cos (p-2a)]/C, etc. 

3.5. Series Developments 

The papers mentioned in the previous sections 
were based on more or less standard procedures, e.g., 
approximating an average with two variables held 
fixed by the average for one fixed variable, etc. 
The present section considers more analytical work. 
We begin by indicating the relationship between the 
heuristic and analytical procedures, and then sum- 
marize some recent results. 

Consider the solution \f/ for the scattering of a 
wave by a given configuration of N scatterers, say 
xp(l, 2, . . . N), where 1 stands for all significant 
properties of scatterer "one" (including its loca- 
tions r) , etc. We introduce an ensemble of configura- 
tions in terms of an appropriate probability distribu- 
tion function W(l, 2, . . . N), and write the average 



of i/', i.e., the coherent field as <l^^> = f . . . fty 
(1, 2, . . . N)W(1, 2, . . . N)d/d2 . . . dN. 
Similarly, we average the intensity \\p\ 2 over the 
ensemble, and write the average total intensity as 
<W 2 > = \<t> l *+V, where |<»| 2 is the coherent 
intensity, and where the incoherent intensity V is 
the average absolute squared deviation of yp from its 
mean value <^\p^>. We may also consider the aver- 
age of the corresponding Poynting vector, and the 
average of other scattering functions. 

In general, we proceed by stating conditions on 
\p(l, 2, . . . N) to make the scattering problem for 
one configuration determinate, then introduce an 
ensemble of configurations in terms of W(l, 2, . . . N), 
and then seek the expectation values <C^> 7 "Cl^l 2 ^ 
etc. We have essentially two kinds of representa- 
tions for \{/: "compact representations" in terms of 
the solutions of a set of simultaneous integral equa- 
tions (say \[/ c ) and "expanded" representations as 
series of orders of scattering (say \p c ). If we work 
with \p e , then we obtain <Cip^> = <i\pe^>, etc., by in- 
tegration; the integrated forms of the terms of the 
series \f/ e may be sufficiently simpler than the orig- 
inals and lead to closed or relatively compact forms 
for <C^c^> (see for example Twersky, 1953a*). More 
usually (practically all the papers cited in the pre- 
vious sections) we average the compact representa- 
tion \p c , proceed heuristically to construct a set of 
"average equations" that are simpler than the orig- 
inals, and then solve these to obtain say <C*Ar^> ~ < C X P^ > - 

A more general discussion than the above, one 
that applies both to a random distribution of scat- 
terers and to a randomly perturbed continuum, is 
given by Keller's [1960] parallel description of "hon- 
est" and "dishonest" methods; Keller uses both 
methods to derive results for scintillations of star- 
light arising from inhomogeneities in the atmosphere. 
The conditions for which <C*/v> approaches <^ e > 
for a class of one-dimensional scattering problems 
have been determined by Bazer [1959*]. 

For three-dimensional distributions of discrete 
scatterers, the relationship of the heuristic c-forms 
to their corresponding e-forms is discussed by Twer- 
sky [1963]. The compact heuristic representations 
<^yp c ^> in the integral equation forms obtained by 
Foldy [1945*], Lax [1951*], and Twersky [1962d] 
differ from <^ e > in that <i\pc^> neglects terms 
corresponding to common path differences in a 
given order of scattering: each scatterer is "used 
only once" in a given order of scattering of <^>, 
so that, e.g., the back and forth interactions be- 
tween a pair are neglected. Solving the two-body 
problems implicit in the general distribution, Twer- 
sky [1963] constructs an alternative heuristic com- 
pact form for <^\p^> to restore some of the missing 
processes ignored by <^c>. Similarly, he extends 
the heuristic integral equation forms <C\^c\ 2 ^> 
derived by Foldy and Lax by taking into account 
additional two-body processes. Comparison of 
<C\>Pc\ 2 ^> with the average of the product of the exact 
series <C|^e| 2 > indicates that the heuristic forms 
^I'/'cl 2 ^ are poorer approximations of <i\$e\ 2 ^> 
than the <^\p c ^> are of <^ e >. In order to obtain 
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a form of <|^ c | 2 > consistent with the approximation 
<^ c >, Twersky introduces into <|^ e | 2 > the same 
simplification that converts <^ e > to <iAc>, i.e., 
each scatterer is "used only once" in a given chain 
of successive scatterings of the single-scattered 
value u. This simplification in the series <|^ e | 2 > 
leads directly to a relatively simple integral equa- 
tion whose inhomogeneous term is \<^\p c ^>\ 2 and 
whose kernel involves a "new single scatterer" 
ue; the function ue satisfies the integral equation 
analogous to <jp c ^> whose inhomogeneous term is 
u, i.e., ue is the corresponding "Green's function" 
(in the sense the term is used in physics) for an 
angle dependent source u. (For the case of mono- 
pole scatterers, Foldy [1945*] obtained two series 
representations for the kernel function in his in- 
tegral equation for <|iAd 2 >; the initial term of one 
of these series corresponds to the new integral 
equation specialized to monopoles.) 

In order to exhibit the relationships between the 
heuristic and exact forms explicitly, Twersky 
[1963] restricts consideration to a very simple class 
of scattering problems: the radiation fields of the 
individual scatterers are significant only in the for- 
ward half -space, particularly near the forward 
scattering direction ; the scatterers are located within 
a slab region of space on which a plane wave is 
normally incident; the scatterers are identical (or, 
equivalently, the averages over isolated scatterer 
parameters are identical), and their position statis- 
tics are those of a "rare gas." For the e -forms, he 
works in general with a finite number of scatterers 
(N) in a finite available volume (V); for the corre- 
sponding c -forms, N and V are infinite but their 
ratio is bounded. This illustrative problem is con- 
structed around scatters which do not need to be 
taken into account more than once in each order of 
scattering. Thus, for this case, the heuristic form 
<i£ c > used explicitly by Foldy, Lax, Twersky, and 
used implicitly by practically all papers in the 
earlier physics literature is the correct limit of <iA e > 
for N— >°o, V— >°o [i.e., we have essentially <C\^c> 
= (l+//A^) iV -^e / =<i/' c >]. However, the situation 
is quite different for <|iA| 2 >- the heuristic form 
constructed with the aid of the energy theorem 
[Twersky, 1962b] gives the correct value that 
<|i// e | 2 > approaches in the limit iV-^°° ; the analog 
of Foldy's [1945] heuristic integral equation <|^ c | 2 > 
for which the kernel function is obtained in closed 
form for this simple illustration, is incorrect; the 
"<i// c >-consistent" heuristic integral equation for 
<|i// c | 2 > is identically the limit of <|iAe| 2 > for 
jy-^oo f and this is also the initial term of the analog 
of one of Foldy 's two series. 
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4.3. Passive Communications Satellites, Review — 1960-1962 

I. Russell Burke ! 



In the period from January 1, I960 to January 1, 
1963 Echo I, the first successful communications 
satellite, was launched on August 12, 1960, and 
successful communications experiments were con- 
ducted using the satellite. In addition, two other 
attempts were made to orbit passive communications 
satellites without resulting in successful communica- 
tions experiments. Of these, one was the first 
Echo launch on May 13, 1960 in which orbit was not- 
attained. The other was the West Ford launch of 
1961. 

The objectives of Project Echo were to demon- 
strate the feasibility of using large, inflatable spheres 
as communications reflectors, and to study the 
behavior of large, lightweight erectable struct nits 
in the space environment. Echo I is described in 
detail in Jaffe [I960]. Figure 1 shows the Echo 
structure inflated on the ground. The satellite was 
made of Mylar 0.0005 in. thick coated with vapor- 
deposited aluminum to provide radio reflectivity. 
The 100-ft. diam sphere weighed 136 lb. and was 
folded into a 26-in. diam container for launching. 
Two sublimating powders were sifted into the 
satellite (20 lb. of anthroquinone and 10 lb. of 
benzoic acid) to cause inflation of the satellite in the 
vacuum of space. Five free space ballistic tests of 
the balloon and inflation system were made using a 
solid fuel booster shown in figure 2 prior to orbital 
launch. These tests developed the packaging, 
ejection, and inflation techniques which resulted in 
the successful deployment of Echo 1 in orbit. The 
orbit achieved by Echo I was circular, 900 nautical 
miles in altitude, inclined 47° to the equator. Com- 
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munications experiments were carried out via Echo I 
between ground stations at Holmdel, N.J., Gold- 
stone, Calif., and Stump Neck, Md., and by inde- 
pendent experimenters located elsewhere in the 
United States and also in Europe. These experi- 
ments and their results have been described in 
detail [Project Echo, 1961; Beamer, 1961; Stevens 
and Victor, I960]. In brief, two-way telephone 
transmissions of excellent quality were made via 
the satellite between the east and west coast of the 
United States. Voice and music were transmitted 
from Holmdel to Jodrell Bank, in England. Fac- 
simile transmission via the satellite also took place. 
Many modulation techniques were employed. 

\o deviations from propagation theory were 
observed, and returned signal strength from Echo I 
was within 1 db of theoretical during the pressurized 
Life of the satellite (1.4 days). 

In addition to the communications experiments, 
Echo I has contributed significantly to understanding 
of the space environment and the causes of orbital 
perturbations; because of its large size and small 
relative weight, Echo I is very sensitive to the effect 
of atmospheric density and solar radiation pressure. 
Jaffe [1961]; Jastrow, and Bryant [1960] summarizes 
results of the Echo I experiments in this regard. 



DIAMETER: 100 FT. WEIGHT: SATELLITE 130 LBS. 
INFLATION MATERIAL 30 LBS. 
CONSTRUCTION: 0.0005 INCH THICK PLASTIC FILM COATED 
WITH VAPOR DEPOSITED ALUMINUM 

Figure 1. 




Figure 2. 
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Figures 3 and 4 respectively show the periodic varia- 
tion of Echo I perigee height with time, due largely 
to solar pressure, during its first year in orbit, and 
the variation of period with time from which atmos- 
pheric density data have been derived. 

Echo I is still in orbit, and although radar observa- 
tions indicate a reduction in cross section of 3 to 6 
db, as recently as May 1962, a television signal 
was transmitted from Massachusetts to California 
via the satellite by Lincoln Laboratory of MIT. 

In Project West Ford, the objective is to orbit a 
belt of half-wave dipoles; in this case a volume of the 
belt illuminated by the ground transmitter would 



serve as the signal source for the ground receiver, 
rather than a single large scatterer such as Echo I. 
Details of the project are contained in [Morrow; 
I960]. 

Echo I, although an excellent isotropic scatterer 
when pressurized, began to display scintillations 
after all of the pressurizing gas had escaped. These 
became larger with time and it has been deduced 
that either the structure is insufficiently rigid to 
resist the environmental forces or that the pattern 
reflecting the satellite's previous folded condition is 
gradually becoming enhanced by the plastic memory 
of the mylar. As a result, NASA has entered a 
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Figure 5. 

second flight program for orbiting spherical passive 
satellites. In this one, Echo II, the satellite will 
be 135 ft in diameter (3 db larger in radar cross 
section than Echo I) and will be constructed of an 
aluminum-mylar-aluminum laminate which, al- 
though only about 50 percent thicker than the Echo I 
material will be almost one hundred times stiffer, 
resulting in a sphere sufficiently rigid to resist both 
the environment and the mylar's plastic memory. 
The material is twice as heavy as the Echo I material; 
as a result Echo II will be approximately four times 
as heavy as Echo I. A detailed description of Echo 
II is contained in Jaffe [1961]. The Echo II orbital 
launch is scheduled for 1963; during 1962 two sub- 
orbital launches were made to test the ejection, 
inflation and deployment systems. Because of the 
heavier weight of the structures, and because it was 
decided to photograph and televise the ejection and 
deployment from the booster, a launch vehicle of 
greater capability was required than for the Echo I 
suborbital tests, and a modified Thor DM-21 was 
employed. This gave a maximum altitude of 1500 
km and a zero g, space vacuum condition for 20 min. 
The first test was conducted on January 15 and the 
deployment was unsuccessful due to excessive 
accelerations during the unfolding process. In the 
second test on July 18, a subliming material of re- 
duced vapor pressure was employed (benzoic acid; 
acetamide had been used for the first test). Ejec- 
tion, inflation and deployment were successful; 
however, radar observations showed scintillations 
as large as ±12db for C-band which indicated that 
the final pressure reached was inadequate to produce 
an un wrinkled surface. Figures 5 and 6 show 
respectively the results of the first and second sub- 
orbital tests a few seconds after the inflation process 
commenced. A controlled inflation system is being- 
developed in which the pressure source used for 
deployment is not that used for final stressing of the 
skin and hence deployment can take place slowly, 
the higher pressures being applied only after the 
sphere is fully deployed. 



DEPLOY 




Figure 6. 

Although there has been only one successful com- 
munications satellite experiment during the period 
reviewed, the general area of passive communica- 
tions satellites has been the source of many papers. 
Some of these present methods for comparing passive 
with active satellites [Campbell and Pollack, I960]. 
However the area having the most interest has been 
demonstrated in methods for increasing the radio 
cross section of passive communications satellites 
without significant penalty in weight or size. Tech- 
niques for accomplishing this can be assigned to 
several categories. The first category involves some 
kind of predetermined orientation. An example of 
this is a spherical segment, vertically oriented which 
can be designed to have a modest gain, and still 
accept transmissions from all visible ground stations 
[Gillespie, I960]. The second category involves 
satellites which do not require orientation. In one 
case, a rather modest gain can be obtained from a 
sphere if the surface is diffuse instead of specular, 
provided that the included angle between ground 
transmitter and receiver does not exceed 83.7° 
[Raabe, 1961]. Other cases include spheres studded 
with Luneberg lenses, corner reflectors, lens-reflector 
combinations, etc. [Stabler, and Johnson, 1961; 
Ryerson, 1960] or represent spheres which are them- 
selves large Luneberg lenses which may be defocused 
to control angular coverage [Siegel et al., 1962]. 
Another case is the use of arrays such as the Van 
Atta array which are relatively insensitive to orienta- 
tion, to reradiate received energy [Hansen, 1961]. 
Another category is represented by proposals to 
increase the size to weight ratio of spherical struc- 
tures, and decrease the response to solar pressure 
through the use of openwork or mesh materials rather 
than a continuous web of conducting material 
[Bradley, 1962; Jastrow I960]. 
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4.4. Passive and Active Reflectors 
J. Kaiser and I. Kay 



In recent years considerable interest has been cen- 
tered on reflectors in space for use as long distance 
communication relays between ground stations. 

One of the earliest papers on this subject was by 
A. C. Clarke [1945] who in 1945 described a global 
communications system based on three uniformly 
spaced equatorial 24-hr orbit communication satel- 
lites. 

Since that time the relative merits of passive 
versus active reflectors in space have been frequently 
discussed in the literature. Active reflectors, gen- 
erally use active electronic elements and some form 
of electric power sources to amplify and repeat in- 
coming signals whereas passive reflectors have no 
active elements or power sources. 

As pointed out by a number of writers [Bartow 
et al., 1959; Bond et al., 1960; Jacoby, 1960; Kaiser, 
and Raisbeck, 1961; Pierce, 1959], the ground sta- 
tions for communications systems based on active 
repeaters require generally much less transmitted 
power and antenna gain than those based on passive 
repeaters of equal weight, orbit altitude and com- 
munication capacity. 

Active reflectors have been proposed for orbit alti- 
tudes from 1000 statute miles (s.m.) to 23,000 s.m., 
the choice of the particular altitude being governed 
by booster capabilities, communication system cover- 
age, circuit availability, system economics, and ap- 
plications. Passive reflectors for practical systems 
would generally be placed into the lower altitude 
orbits, unless they can be made extremely large 
[Bradley, 1961]. 

An interesting case is the moon, which is being 
used by the U.S. Navy for a 10-channel teletype com- 
munications system between Hawaii and Washing- 
ton, D.C. The effective scattering cross section of 
the moon is found to be 9X10 n m 2 , corresponding to 
a perfectlv conducting sphere 670 miles in diameter 
[Jakes, 1961]. 

Experimental communication systems based on 
active and passive reflectors in space have been 
tested. Telstar and Relay are examples of active 
reflectors or repeaters. The moon and the Echo I 
balloon are examples of passive reflectors. The ex- 
perimental performance of the Echo I balloon agreed 
with the predictions for a perfectly reflecting sphere 
of the same (100 ft) diameter. 

Of the active reflector designs which have been 
considered, probably the most novel is due to Van 
Atta [Sharp, 1958; Sharpe et al., 1960; Van Atta, 
1959]. The Van Atta array consists of antenna ele- 
ments interconnected by transmission lines whose 
electrical lengths are so chosen that the phase dis- 
tribution of a plane wave incident at an arbitrary 
angle will automatically recombine to generate a 
transmitted plane wave in the opposite direction. 
Modifications of the Van Atta array were considered 
by Bahret [1961; Bauer, 1961; Waneslow, 1962, and 
Rutz and Kramer, 1962]. An evaluation of the Van 



Atta array for communications was made by Hansen 
[1961] who concluded that its usefulness would be 
for long distance point-to-point transmission. 

A number of workers have analyzed the possibility 
of enhancement of the signal return from passive 
reflectors in orbit. J. R. Pierce and C. C. Cutler 
[1959] list several forms of passive reflectors which 
may be useful for satellite radio relay, including bal- 
loons, saucers, corner reflectors, clustered poly- 
hedrons, and other shapes. (See also M. G. Chate- 
lain [I960].) 

The prototype of the corner reflector, whose cross 
section is characterized by a pronounced dependence 
on wavelength, is the flat plate. At t lie other 
extreme, the perfectly reflecting sphere has been 
regarded as a standard for wide-angle reflectors 
since in the optical limit the sphere is an isotropic 
scatterer. Wide-angle reflectors such as the per- 
fectly reflecting sphere are characterized by cross 
sections which are independent of wavelength when 
the wavelength is sufficiently short. 

Passive reflectors having a relatively wide angle 
of return, whose cross section does, however, depend 
on frequency have been discussed by W. E. Morrow 
[1960, 1901J. A large number of dipoles randomly 
oriented in orbit would present a maximum realiza- 
ble average scattering cross section -A 5 =0.16X 2 where 
X=wavelength. 

Although the corner reflector docs return a large 
amount of power in one direction, it would be 
necessary to provide attitude control in order to use 
it as a relay because of its highly directional reflection 
characteristics. On the other hand, while the per- 
fectly reflecting sphere provides wide coverage its 
cross section is considered to be disappoint ingly small 
for its practical use as a satellite relay. 

The attempts to surmount these difficulties also 
fall into two general categories. One attempt in- 
volves the use of a large number of corner reflectors 
either randomly placed or attached symmetrically 
to a large sphere so that the attitude control is no 
longer necessary [Bradley, 1961; Cutler, 1962; Har- 
rington, 1962;Herriott, 1959; Morrow, 1961; Sinclair, 
1959; Stabler and Johnson, I960]. The other is to 
use some form of dielectric lens which orients the 
reflected ravs in suitable directions [Bobev, 1959; 
Graham et ah, 1961; Kay, 1959; Reyes, 1959] order 
to avoid the problem of attitude control in the latter 
type of reflector the lens is usually required to be 
spherically symmetric. It can be shown that there 
is a maximum cross section which can be obtained 
over a given angle for any spherically symmetric 
passive reflect or. Thus the coverage requirement of a 
system must necessarily limit the return from such a 
reflector of a specified diameter. 

Y. E. Stabler and A. L. Johnson [1960] have 
analyzed and experimented with reflectors having 
increased reflectivity. They have concluded that 
considerable improvement in directivity, above 13 
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db, can be obtained by placing reflecting elements or 
lenses combined with reflecting surfaces on spherical 
satellites. In the 1962 spring URSI meeting, R. 
Graham, K. M. Siegel, et al., [1961] discuss the 
optimum parameters for a spherically symmetric 
corner reflector. They obtained the best results 
for a solid spherically symmetric reflector, which is a 
Luneberg lens with a uniform coating of power 
reflection coefficient approximately equal to one- 
third and a power transmission coefficient of two- 
thirds. Such a lens was actually constructed and 
an experiment performed during the meeting to 
demonstrate the result. They estimated that the- 
oretically for a lens the size of the Echo balloon the 
backscattering cross section at X-band is 60 db 
above that from Echo. 

C. C. Cutler [1962] discusses some of the funda- 
mental limitations on the possibility of obtaining 
useful directivity from a nonoriented, passive re- 
flector in orbit. He derives a formula for the net 
advantage of a reflectivity enhanced scatterer such 
as a sphere covered with small reflecting elements 
over the reflectivity of a smooth sphere. The ad- 



vantage is equal to 4 



/sin</>Y 



2 where 20 is the apex 



angle of a cone within the sphere which reflects the 
impinging waves and 6 is the half-power beamwidth 
of the reflected beam. He shows by an example 
that for a satellite height of 4,000 miles and ground- 
station separation of 4 ; 000 miles the gain is less than 
unity and at a height of 22,300 miles the gain is 
8.5 db. He concludes that for the amount of gain 
obtained over that of the smooth sphere, the prob- 
lems encountered with increased directivity scatters 
will probably negate their usefulness. 

J. Kaiser and G. Raisbeck [1961] analyzed and 
evaluated passive communication satellite systems. 
These authors state that an important parameter 
for the evaluation of passive reflectors for com- 
munications is the effective area of the reflector 
which determines the power received from the 
ground transmitting station. It is noted that ex- 
tremely large areas are required to receive sufficient 
power to compete favorably with active communi- 
cation satellites having a radiated power of 1 w. 
W. L. Bradley [1961] discusses orbiting passive 
electromagnetic reflectors constructed of fine metallic 
wire mesh, which, if engineering problems of un- 
folding and erection could be solved, would allow 
large scattering cross section reflectors having 
relatively small weight. 

A general conclusion which can be drawn from 
the various studies made on active and passive 
reflectors is that at the present time active reflectors 
seem to be the preferred relays for most communi- 
cation systems, while at the same time there do 
exist special situations for which passive reflectors 
are probably more efficient. 

Although the most dramatic application of passive 
and active reflectors seems at present to be in the 
field of long-distance communication systems, this 
report would not be complete without some mention 



of the many general purpose studies which have 
been made of radar enhancement devices. Most of 
these studies have been concerned with passive re- 
flectors, and some very interesting results have been 
obtained. For example, D. Atlas has observed that 
a homogeneous dielectric sphere having a diameter 
less than 6.4 times the wavelength will, if it is capped 
at its rear pole by a small piece of metal, have a 
backscattering cross section in excess of that pro- 
duced by a Luneberg lens reflector of the same 
diameter [Atlas, 1962]. 

Methods for designing reflectors having desired 
reflection properties were studied by S. T. Bobey 
[1959] L. Peters [1959a, 1959b] and Kennaugh and 
Bauemler [1958] of the Ohio State University Re- 
search Foundation. A synthesis procedure for 
inhomogeneous spherically symmetric radar en- 
hancement lenses was prescribed by A. F. Kay 
[1959] who made a general analysis of this type of 
reflector. 

Twersky [1961] discussed coherent and incoherent 
scattering by orbital distributions of scatterers. 
He derived mathematical expressions for scatterers 
distributed in regions of space bounded by various 
surfaces and illustrated the results with an example 
from geometrical optics. 

In addition to such efforts many cross-section 
calculations and measurements were made for di- 
electric and metallic reflectors for particular shapes 
by other investigators. Some papers which describe 
work of this type are included in the reference list 
below although no specific mention of them has 
been made elsewhere in this report. 
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